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Preface 


This publication was formerly entitled The Deep Space Network Progress Report. 
Although the practice of reporting progress in the development and operations of the 
Deep Space Network continues, the report has been expanded to include developments in 
Earth-based radio technology as applied to other research programs. These programs are: 

(1) Geodynamics: For several years, the laboratory has been developing radio 
interferometry at microwave frequencies for application to geodetic measure- 
ments. This branch of telecommunications technology is now being applied to the 
study of geodynamics. 

(2) Astrophysics: The deep space stations, individually and in pairs as an inter- 
ferometer, have been used by radio astronomers for astrophysics research by 
direct observations of radio sources. 

(3) An activity closely related to radio astronomy’s use of the deep space stations is 
NASA’s continuing program of radio search for extraterrestrial intelligence in the 
microwave region of the electromagnetic spectrum. 

Each succeeding issue of this report will present material in some, but not all, of the 
following categories: 

Radio Astronomy 

Search for Extraterrestrial Intelligence 
Radio Interferometry at Microwave Frequencies 

Geodetic Techniques Development 
Spacecraft Navigation 

Orbiting Very Long Baseline Interferometry 

Deep Space Network 

Description 
Program Planning 

Planetary and Interplanetary Mission Support 
Advanced Systems 

Network and Facility Engineering and Implementation 
Operations 

Spacecraft Radio Science 
Planetary Radar 
Energy 

In each issue, there will be a report on the current configuration of one of the seven 
DSN systems (Tracking, Telemetry, Command, Monitor and Control, Test Support, 
Radio Science, and Very Long Baseline Interferometry). 

The work described in this report series is either performed or managed by the 
Telecommunications and Data Acquisition organization of JPL. 
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SETI Downconverter 


B. Crow 

Communications Systems Research Section 


The SETI (Search for Extraterrestrial Intelligence) breadboard subsystem was begun in 
order to develop the technology to efficiently implement a SETI instrument capable of 
searching wide bandwidth with high resolution. Due to the recent directive from Congress 
to terminate the program, NASA has tasked JPL to complete the SETI breadboard (since 
it is more than half complete) so that it will be available for use by a NASA-selected user 
by the end of FY 82. The downconverter covered in this report is the interface hardware 
between the receiver IF output and ADC used as the input to the spectrum analyzer. 


The SETI (Search for Extraterrestrial Intelligence) down- 
converter has been developed to interface between an IF input 
and provide complex baseband outputs for use in the SETI 
breadboard system (see Fig. 1). 

The SETI downconverter subassembly (see Fig. 2) has been 
built in a modular fashion with 13 modules in the SETI down- 
converter drawer and the required monitor and control and 
power supplies in a second drawer. The downconverter accepts 
both a 30-MHz IF (for use at an Arecibo receiver) and a 
50-MHz IF (for use with a Block III DSN receiver). The 
50-MHz IF is translated to 30 MHz in the downconverter. Two 
types of signal gain control are available: (1) an AGC loop 
with loop time constants of 0.1,1, and 10 sec, and (2) a com- 


puter control gain control that offers a 35-dB dynamic range. 
An ADC is provided to measure the output of the AGC detec- 
tor and is calibrated to provide an indication of output power. 
The IF signal is basebanded in the complex mixer where care- 
ful attention has been paid to the amplitude and phase balance 
between the two baseband signals (see Fig. 3). An internal test 
signal is available to allow stand-alone testing of the downcon- 
verter (and the SETI breadboard). The downconverter is con- 
trolled and monitored by a PDF 1 1/44 and is currently operat- 
ing as designed in the development laboratory. 

Table 1 is a list of the technical requirements and the mea- 
surements that have been made to confirm the downconverter 
performance. 
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Table 1 . Downconverter requirements vs measured performance 


Parameter 

Requirement 

Measurement 

Input bandwidth 

10 MHz 

16 MHz 

Input frequency 

30 MHz 
50 MHz 
Test signal 

30 MHz 
50 MHz 
28 MHz 

Baseband output 

.1 VRMS 

.1 VRMS 

AGC dynamic range 

>30 dB 

50 dB 

AGC time constant 

0.1, 1, 10 sec 

0.1, 1, 10 sec 

Compute controllable gain 
range 

>30 dB 

35 dB (see 
Fig. 4) 

Input impedance 

50 ±5 ohm 

58 ohm 

Baseband frequency 
response 

4 MHz 

4.3 MHz 

Rejection of unwanted 
harmonic 

>30 dB 

2nd harmonic = 
-46 dB 

3rd harmonic = 
-52 dB 

Phase and amplitude 

90 ±3 

See Fig. 3 

balance 

0 ±0.5 dB 

See Fig. 3 

Linearity 

Pin max <-25 dBm 

See Fig. 7 

Noise figure 

- 

See Fig. 6 

ADC power readout 

- 

See Fig. 5 


30-MHz 


50-MHz 



Fig. 1. SETI breadboard 


50-MHz IN 
30-MHz IN 


5-MHr 

REF 


20-MHz 

REF 



Fig. 2. SETI downconverter block diagram 
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The DSN Radio Science System Mark III - 1980 

B. J. Buckles 
TDA Engineering Office 


The DSN Radio Science System supported the Voyager 2 Saturn encounter radio 
science experiments in August 1981. Support for these experiments was provided by all 
the Deep Space Stations of the DSN. However, the critical support for the Saturn occulta- 
tion and ring scattering experiment was provided at DSS 43 by the medium-band open- 
loop recording system. The implementation of the medium-band open-loop recording 
system at DSS 43 culminated many years of effort by the DSN to provide a DSN multi- 
mission radio science capability, and all 64-m DSS open-loop recording systems are now 
essentially identical This article describes the DSN Radio Science System, emphasizing 
the recent implementation at DSS 43. 


I. Introduction 

Figure 1 is a block diagram representation of the DSN 
Radio Science System (RSS) recently used to successfully 
support the Voyager 2 Saturn encounter radio science experi- 
ments. The DSN Deep Space Stations (DSSs) are the heart of 
this system. In particular, DSS 43 was the prime supporting 
station for the radio science experiments at Voyager 2 Saturn 
encounter. DSS 43 supported the Saturn planetary occultation 
and ring scattering experiments with a newly implemented 
medium-band recording system consisting of a highly phase- 
stable, four-channel, open-loop receiver and a high-rate, high- 
density digital recorder. Also, a wideband recording system 
consisting of a four-channel, wide-bandwidth, open-loop 
receiver and digital recorder was installed at DSS 43 to serve 
as a backup to the medium-band system. 

Other new implementations provided to support the Voy- 
ager radio science experiments were modifications to the 


Precision Power Monitor (PPM) to provide more reliable 
operation of the instrument and new operational software for 
the Occultation Data Assembly (ODA) that provided an 
interface to the PPM, allowing PPM data to be recorded by the 
ODA and, thus, system temperature information to be pro- 
vided with the Original Data Record (ODR) recordings. Also 
the capability to play back data from the DSS over the wide- 
band data lines for system testing purposes was provided. 

The following sections describe the DSS subsystems that 
support the DSN Radio Science System in detail. In addition, 
each of the other facilities supporting the DSN Radio Science 
System, the Ground Communications Facility (GCF) and the 
Network Operations Control Center (NOCC) and their support 
functions are described. The narrative in general will describe 
the flow of data from receipt at the antenna through the sys- 
tem as indicated on the Radio Science System block diagram 
(Fig. 1). 
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II. The Deep Space Stations 

The Deep Space Station (DSS) is the Radio Science System 
instrument. DSS system performance directly determines the 
degree of success of the experiment, and system calibration 
determines the degree of uncertainty in the results of the 
experiment. The following paragraphs describe those functions 
performed by the individual subsystems depicted in Fig. 1. 
Specific configuration and calibration requirements are ad- 
dressed in a separate paragraph after the functional descrip- 
tions of the subsystems. 

A. The Antenna Mechanical Subsystem 

The 64-m antenna mechanical subsystems function as 
large-aperture collectors which, by double reflection, focus 
incoming radio frequency (S- and X-band RCP/LCP) energy 
into the S- and X-band feedhorns (part of the Antenna Micro- 
wave Subsystem). 

The large collecting surface of the antenna focuses the 
incoming RF energy onto the hyperboloid subreflector which 
is adjustable in both axial and tilt positions to permit optimi- 
zation of the focusing of energy into the feedhorns. The sub- 
reflector is locked into a fixed position for many radio science 
events. The subreflector then reflects the received energy to 
the dichroic plate, a device that reflects S-band energy to the 
S-band feedhorn and transmits X-band energy to the X-band 
feedhorn. 

Transmitted S-band RF energy emanating from the feed- 
horn is focused by the same reflectors into a narrow cylindri- 
cal beam. Since the beam is very narrow, the antenna must be 
pointed with high accuracy and precision. This is accomplished 
by a series of drive motors and gear trains which rotate those 
portions of the structure that support the reflectors, position 
sensors, and related electronics. Electronic servo amplifiers 
are used to amplify and condition the axes angles or position 
error signals which are received and are provided to the drive 
motor controls. 

Pointing angles are computed from an ephemeris provided 
by the Project, and the antenna is pointed to these angles. 
Once the receiver has acquired a signal to provide feedback, a 
radio source can be tracked by scanning around it (CONSCAN) 
and computing pointing angles from signal-level information 
supplied by the receiver. During periods when excessive signal- 
level dynamics or low received signal levels are expected, i.e., 
occultations and some other radio science experiments, 
CONSCAN cannot be used and angle pointing is accomplished 
by manually inserting offsets to the computed angle predict 
set. 


B. The Antenna Microwave Subsystem 

The antenna microwave subsystem accepts the received S- 
and X-band RCP/LCP signals at the feedhorn from the antenna 
mechanical subsystem. The received signals are transmitted 
through the polarizer plates to the orthomode transducer. 
The pojarizer plates are adjusted so that RCP signals are 
directed to X-band traveling wave maser 2 and S-band traveling 
wave maser 1 and so that LCP signals are directed to X-band 
traveling wave maser 1 and S-band traveling wave maser 2. 
After amplification by the masers, the signals are routed to 
the receiver-exciter subsystem via the microwave switching 
assembly. 

The S-band uplink signal is transmitted via the diplexer 
assembly through the feedhorn to the antenna where it is 
focused and beamed to the spacecraft. The noise diode assem- 
blies under control of the PPM inject known amounts of noise 
power into the received signal path so that accurate real-time 
system temperature measurements may be made. 

C. The Transmitter Subsystem 

The Transmitter Subsystem accepts the S-band frequency 
exciter signal from the receiver-exciter subsystem and ampli- 
fies it to a transmitted output level of 20 kW. The signal is 
then routed via the diplexer to the antenna and then focused 
and beamed to the spacecraft. 


D. The Receiver-Exciter Subsystem 

The Receiver-Exciter Subsystem receives, amplifies, and 
frequency downconverts spacecraft-radiated S- and X-band 
RCP/LCP signals. The closed-loop receivers provide doppler 
and ranging signals to the tracking subsystem. Dedicated open- 
loop receivers provide baseband signals to the radio science 
subsystem, and fixed-tuned, wide-bandwidth, open-loop 
receivers provide backup systems to the narrowband and 
medium-band open-loop receivers. The exciter generates the 
S-band drive signal provided to the transmitter subsystem 
that provides the spacecraft uplink signal. 

The Spectral Signal Indicator provides local displays of 
received signal spectrums and routes spectral data to the radio 
science subsystem. These displays are used to validate radio 
science system data at the DSS, NOCC, and project areas. 

The Precision Power Monitor (PPM) measures system 
temperatures by injecting known amounts of noise power 
into the signal path and comparing the total power within a 
given bandwidth before and during periods of noise injection. 
System temperature measurements are made for each of the 
four masers by utilizing the closed-loop receivers as monitoring 
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devices. This use of the receivers imposes a configuration 
restraint on them when they are used to monitor all four 
masers. That is, the Block IV receivers must monitor X-band 
RCP/LCP masers and the Block III receivers are used to 
monitor S-band RCP/LCP masers. 

E. The Tracking Subsystem 

The Tracking Subsystem receives the ranging spectrum and 
doppler signals from the Receiver-Exciter Subsystem. It gener- 
ates a range code that is routed to the exciter and modulates 
the S-band uplink carrier. The demodulated range spectrum is 
compared to a model of the transmitted range code, and the 
round trip signal delay to the spacecraft is computed, for- 
matted and transmitted to the GCF Data Records Subsystem, 
which produces an Intermediate Data Record (IDR) tape upon 
which the data are delivered to the Project. Similarly the 
doppler phase is counted, formatted, and delivered to the 
Project. The radio metric data supports radio science celestial 
mechanics and general relativity experiments in addition to 
providing a partial backup to the open-loop system. 

In addition, the Tracking Subsystem computes data resid- 
uals and noise estimates, receives and stores predicts, and 
provides partial status information to the monitor and control 
subsystem including receiver AGC levels and lock status. 

F. The Radio Science Subsystem 

The Radio Science Subsystem digitizes, bandwidth reduces, 
and records (l)very narrow, (2) narrow, and (3) medium- 
bandwidth radio science data and digitizes and records wide- 
bandwidth radio science data. It receives radio science fre- 
quency predicts from NOCC, configuration and control data 
from the monitor and control subsystem, and S- and X-band 
RCP/LCP signals from the Receiver-Exciter Subsystem. It 
transmits spectral information from the SSI to NOCC via the 
GCF wideband data lines. It controls the narrow- and 
medium-band open-loop receiver local oscillator (LO) by 
sending frequency control information to the Receiver- 
Exciter Subsystem. 

G. Monitor and Control Subsystem 

The Monitor and Control Subsystem provides control 
messages to the centrally controlled DSS subsystems including 
the Radio Science and Tracking Subsystems. It provides 
partial status information for these same subsystems. It gen- 
erates the angle drive tape needed by the Antenna Pointing 
Subsystem and, in addition, controls the SSI and PPM in the 
Receiver-Exciter Subsystem. It receives the system tempera- 
ture information from the PPM and displays it on the Data 
System Terminal (DST). It relays the system temperature 


information to the NOCC Monitor Subsystem for display in 
NOCC and the project mission support area. 

III. The Ground Communications Faciiity 

The Ground Communications Facility provides the com- 
munication networks required to support the communication 
requirements of the Radio Science System. These facilities 
exist at the DSS and JPL and are briefly described in the 
following paragraphs. 

A. GCF High-Speed Data Subsystem 

The High-Speed Data Subsystem transmits radio science 
predictions from the NOCC to the DSS and CTA21, and 
Radio Science System performance validation data from the 
DSS to NOCC. 

B. GCF Wideband Subsystem 

In real-time, the Wideband Subsystem transmits SSI data 
from the DSS to the NOCC. 

C. GCF Data Records Subsystem 

The GCF Data Records Subsystem formats and provides 
radio science data on computer-compatible tape to the flight 
projects. 

IV. Network Operations Control Center 

A. NOCC Radio Science Subsystem 

The NOCC Radio Science Subsystem (NRS) generates 
open- and’ closed-loop radio science DTV graphics displays, 
and DRS status and configuration displays. In addition, the 
NRS provides the planetary atmosphere refracted trajectory 
input to the NOCC Support Subsystem for the purpose of 
generating radio science predicts. 

B. NOCC Monitor Subsystem 

The NOCC Monitor Subsystem displays system temperature 
information and supports the DSS Monitor and Control Sub- 
system. 

C. NOCC Support Subsystem 

The NOCC Support Subsystem generates DSS frequency 
and tracking predicts using a polynomial coefficient tape 
produced by the Planetary Orbital Ephemiris Atmospheric 
Study (POE AS) software as input. In addition, predicts can be 
generated using manual inputs. 
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D. NOCC Display Subsystem 

The NOCC Display Subsystem provides the NOCC Radio 
Science Subsystem generated radio science displays to the 
Network Operations Control Area and to project radio science 
areas, and provides control data to the NOCC Radio Science 
Subsystem. 

E. NOCC Tracking Subsystem 

The NOCC Tracking Subsystem generates tracking system 
status displays and transmits them to the Display Subsystem 
to be routed and displayed. 

V. CTA21 

A unique part of the Radio Science System exists at 
CTA21. This equipment bandwidth-reduces wide-bandwidth 
radio science data recorded at the DSS. It receives digital wide- 
bandwidth and medium-bandwidth radio science data from the 
DSS Radio Science Subsystem, generates bandwidth-reduced 
radio science data, and provides medium- and wide-bandwidth 
radio science data on computer-compatible tape to the flight 
project. 

VI. Open-Loop Recording System Pre-Pass 
Calibration 

The recorded open-loop receiver signals were the prime 
deliverable data type for the Voyager 2 Saturn occultation 
and ring scattering radio science experiments. For that reason 
it was extremely important that the open-loop system be 
properly set up prior to the recording pass, and that a calibra- 
tion tape be made. This calibration data serves as a bench- 
mark of system performance and is later used to determine 
data quality and establishes the basic uncertainty in the 
results of the experiment. 

VII. Purpose of the Calibration 

The open-loop recording calibration tapes are made to 
accomplish the following: 

(1) Establish the output of the open-loop receivers at a 
level that will not saturate the input to the Occultation 
Data Assembly (ODA) or Digital Recording Assembly 
(DRA) A/D converters. 


(2) Provide data that will establish baseline values for the 
phase, frequency, and amplitude stability of the open- 
loop system. 

To achieve goal 1 , the calibration recording is made using a 
test signal generated by the exciter-translator that is set to 
provide’^the maximum predicted received signal level for the 
recording period. Then the output level of the receivers is 
adjusted to the level determined by the equation 


_(2SNR + 2k^/l^ + k^)\ 

where 

S = receiver output levels (volts rms) 

SNR = expected received signal-to-noise ratio 

k = number of sigma margin desired for saturation 
L - A/D converter saturation level (volts) 

In order to achieve goal 2 for medium-band recording, it is 
necessary to step the test signal across the bandpass of the 
receiver filter. This is done in ten 5-khz (at S-band) steps. 
Stepping the signal in frequency in this manner allows the 
necessary phase, frequency, and amplitude calibrations to be 
made. Frequency stepping is not necessary for narrow band 
recording. 

So that a X-band RCP/LCP test signal may be generated a 
cap is provided that is placed over the X-band feedhorn from 
which a linearly polarized signal is injected into the horn. 
This cap is oriented such that the linearly polarized signal 
is transmitted equally by the orthomode coupler to its RCP 
and LCP output ports, thus providing equal RCP/LCP calibra- 
tion signals. 

A consequence of using a cap over the feedhorn in this 
manner to generate a test signal is that the X-band system 
temperature has been increased approximately 10 dB. It has, 
in effect, been terminated into an ambient load. To counter 
the 10-dB increase in receiver output power due to the pre- 
sence of the cap it is necessary to insert 10-dB pads into the 
signal path at appropriate points. These must be removed after 
the calibration recording is made. 
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Fig. 1. Deep Space Network Radio Science System block diagram depicting the Radio Science System as configured to support the Voyager 2 Saturn encounter 
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This article is a summary of all spaceflight missions that the Deep Space Network has 
supported, beginning with the Ranger Program and ending with the continuing coverage 
of the Voyager Program, 


I. Introduction 

The Deep Space Network (DSN) is managed and operated 
by the Jet Propulsion Laboratory under a NASA contract and 
has been since NASA was formed in 1958. The DSN is tasked 
to support all deep space planetary probes, including the lunar 
exploration. 

Since the Ranger series in the 1960s the DSN has evolved to 
a highly sophisticated telecommunications and data acquisi- 
tion network. The 26-meter antennas have been converted to 
34 meters, and 64-meter antennas have been built at all com- 
plexes. The L-band systems were changed to S-band, with 
X-band now also being used on Voyager and Pioneer 12. 
Ranging systems were developed and are now used along with 
other techniques to determine exact spacecraft locations. 
Many other developments occurred during this period to sup- 
port the requirements of future space flight missions. 

II. The Programs 

A. Ranger Program 

The Ranger missions were a series of lunar probes to obtain 
pictures from different areas of the Moon, study space physics 


in transit, and analyze the human environment from a landed 
instrument package. The primary objective was the collection 
of preliminary information for scientific studies of possible 
landing sites for the NASA manned lunar program. 

Several major technical advances were first used on Ranger. 
First, the spacecraft was fully stabilized in attitude, keeping a 
high-gain antenna pointed at Earth and solar cell panels aimed 
at the Sun, which was a technical challenge in the early 1960s. 
The second technical advance attempted, but not completed, 
was the landing of an instrument capsule on the lunar surface 
by means of a retrorocket. 

B. Surveyor Program 

The objective of these missions was the softlanding on the 
lunar surface of instrumented spacecraft capable of performing 
operations that would contribute scientific knowledge about 
the lunar surface. The Surveyors also returned data to explore 
landing sites in support of the Apollo Project. 

C. Lunar Orbiter Program 

These were photographic missions to obtain high-resolution 
pictures of the lunar surface while in orbit around the Moon. 
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The primary purpose was to search for and survey landing sites 
for the Apollo Project. 

D. Mariner Program 

This series of flight missions consisted of 10 different 
launches, of which three failed, that flew by Venus, Mars, 
Mercury, and orbited Mars, The objectives of these missions 
were to make interplanetary and planetary measurements con- 
sisting of the different science phenomena, and except for the 
first 2 Mariners, to return TV pictures of the different planets. 
The Mars Orbiter mapped the planet and returned data for site 
analysis for the Viking Landers. 

E. Pioneer Program 

The primary objective of the Pioneer 6 through 9 missions 
is to collect scientific data relative to interplanetary phenom- 
ena within a region approximately 0.8 to 1.2 astronomical 
units (AUs) from the sun. Phenomena of particular interest 
include the characteristics of electric and magnetic fields; 
electron density along the earth-spacecraft path; and temporal 
and spatial distribution of plasma, cosmic rays, high-energy 
particles, and cosmic dust. 

The primary objectives of the Pioneer 10 and 11 missions 
are to conduct investigations of the nature of the Asteroid 
Belt, the physical characteristics of Jupiter and its imme- 
diate environment, and the interplanetary medium beyond 
the orbit of Mars and to the extreme of the spacecrafts’ 
communications capabilities. Additionally, Pioneer 11 has 
a primary objective to investigate the physical character- 
istics of Saturn and its immediate environment. 

The Pioneer Venus multiprobe spacecraft investigated the 
characteristics of the Venusian atmosphere down to the lowest 
scale height above the surface. 

The Bus was retarded after the probes were released so that 
it did not enter the atmosphere until after completion of the 
probes’ descent to the surface. In this way, the Bus acted as a 
frequency reference for an interferometric determination of 
the atmospheric wind effects on the probes during their 
descents. 

The objectives of the Orbiter spacecraft were to conduct 
investigations of the characteristics of the upper atmosphere 
of Venus and the region about the planet. Additionally, the 
spacecraft performed radar altimetry and imaging of the 
planetary surface. 

F. Helios Program 

The objective of this mission was to place two solar orbiting 
spacecraft in highly elliptical orbits to achieve a perihelion of 


0. 3 AU and aphelion of 1 .0 AU. With these trajectories the 
spacecraft were able to study the properties of the Sun at 
close range, and enable better understanding of the Sun’s 
influence upon the Earth. One of the spacecraft is still active 
after 7 years in a changing thermal and radiation environment. 
This Program is a joint space venture between the Federal 
Republic of West Germany and the United States of America. 

G. Apollo Program 

This manned spacecraft mission was able to place the first 
men on the lunar surface. The DSN had a secondary role to 
the prime NASA tracking network, and was able to help sup- 
port the Apollo 8 through Apollo 1 7 Missions. 

H, Viking Program 

The primary objectives of the Viking spacecraft were to 
study the planet Mars from orbit and the surface by the use 
of two Orbiters and two Landers. This was done while the 
Vikings took more than 10,000 pictures from orbit and the 
surface; conducted a search for microbial life; performed 
organic and mineralogical studies of the Martian soil; studied 
physical and magnetic properties of the soil; made daily 
weather reports and listened for seismic activity. An extended 
mission known as the Lander Monitor Mission (LMM) is still in 
operation. The primary objectives of the LMM are to continue 
radio science experimentation using Lander 1 as a radio source 
on Mars, to continue to monitor Mars’ meteorology and to 
continue periodic picture taking of Martian terrain around 
Lander 1 . 

1. Voyager Program 

The primary objectives of the Voyager Mission were to 
conduct exploratory investigations of the Jupiter and Saturn 
planetary systems and the interplanetary medium between 
Earth and Saturn. This was accomplished by two launches in 
1977 on flyby trajectories that employ a Jupiter’s gravitational 
assist to reach Saturn.* The primary science objectives were to 
conduct comparative studies of the planetary systems of 
Jupiter and Saturn, including their environment, atmosphere, 
surface, and body characteristics. Additional objectives in- 
cluded the investigation of one or more satellites of each 
planet, the nature of Saturn’s rings, and the interplanetary and 
interstellar media throughout the cruise phase of the mission. 
The Voyager Mission designs included one spacecraft, 
Voyager 2, which is using Saturn’s gravitational assist to 
go on to a Uranus and a Neptune flyby. 

III. Summaty 

All space flight missions previously described in this article 
are summarized in Table 1 . 
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Table 1. Space flights supported by the Deep Space Network 


Mission 

Launch date 

Objective 

Results 

Ranger 1 

23 Aug 6 1 

Lunar flyby 

Decayed in Earth orbit 

Ranger 2 

18 Nov 61 

Lunar flyby 

Decayed in Earth orbit 

Ranger 3 

26 Jan 62 

Soft land seismometer Passed within 22,862 miles 

Ranger 4 

23 Apr 62 

Soft land seismometer Impacted on backside 

Ranger 5 

18 Oct 62 

Soft land seismometer Passed within 450 miles 

Ranger 6 

30 Jan 64 

Lunar photo-impact 

TV system failed 

Ranger 7 

28 Jul 64 

Lunar photo-impact 

Returned 4316 photos 

Ranger 8 

17 Feb 65 

Lunar photo-impact 

Returned 7137 photos 

Ranger 9 

21 Mar 65 

Lunar photo-impact 

Returned 5814 photos 

Surveyor 1 

30 May 66 

Lunar soft landing 

Successful engineering flight; 
over 11,000 photos 

Surveyor 2 

20 Sep 66 

Lunar soft landing 

Failed due to vernier engine 

Surveyor 3 

17 Apr 67 

Lunar soft landing 

Made 3-bounce landing, measured 
surface bearing; 6326 photos 

Surveyor 4 

14 Jul 67 

Lunar soft landing 

Lost during retro maneuver 

Surveyor 5 

8 Sep 67 

Lunar soft landing 

Provided soil analysis; over 
19,000 photos 

Surveyor 6 

7 Nov 67 

Lunar soft landing 

Performed first lateral hop ; 
over 30,000 photos 

Surveyor 7 

7 Jan 68 

Lunar soft landing 

Scientific exploration of a 
highland region; 21,038 photos 

Lunar Orbiter 1 
Lunar Orbiter 2 

10 Aug 66 
6 Nov 66 

Equatorial orbit j 
Equatorial orbit 

1 Photographed complete equatorial 

r belt of Moon for Apollo. First 

Lunar Orbiter 3 

5 Feb 67 

Equitorial orbit j 

1 Earth pictures from deep space 

Lunar Orbiter 4 

4 May 67 

Polar orbit ] 

1 Photographed approximately 98% 

Lunar Orbiter 5 

1 Aug 67 

Polar orbit J 

I of lunar surface 

Mariner 1 

22 Jul 62 

Venus flyby 

Destroyed by range safety officer 

Mariner 2 

27 Aug 62 

Venus flyby 

Passed within 21,594 miles 

Mariner 3 

5 Nov 64 

Mars photo flyby 

Shroud failed to separate 

Mariner 4 

28 Nov 64 

Mars photo flyby 

Returned 22 photos 

Mariner 5 

14 Jun 67 

Venus flyby 

Obtained atmospheric, 
ionospheric, other data 

Mariner 6 

24 Feb 69 

Mars flyby 

Studied surface and atmosphere: 
returned TV images 

Mariner 7 

27 Mar 69 

Mars flyby 

Studied surface and atmosphere: 
returned TV images 

Mariner 8 

8 May 71 

Mars orbiter 

Failed to achieve Earth orbit 

Mariner 9 

30 May 71 

Mars orbiter 

Mapped nearly all of planet’s 
surface; determined landing sites 
for Viking 

Mariner 10 

3 Nov 73 

Venus/Mercury flyby 

Studied clouds of Venus, returned 
TV images of Mercury: first to 
use gravity of one planet as aid in 
reaching another 


Table 1 (contd) 


Mission 

Launch date 

Objective 

Results 

Pioneer 6 

16 Dec 65 

Investigation of inter- 
planetary environment 

Obtained data on solar flares, 
solar wind, interplanetary parti- 
cles and fields 

Pioneer 7 

17 Aug 66 

Investigation of inter- 
planetary environment 

Obtained data on solar flares, 
solar wind, interplanetary 
particles and fields 

Pioneer 8 

13 Dec 67 

Investigation of inter- 
planetary environment 

Obtained data on solar flares, 
solar wind, interplanetary parti- 
cles and fields 

Pioneer 9 

8 Nov 68 

Investigation of inter- 
planetary environment 

Obtained data on solar flares, 
solar wind, interplanetary parti- 
cles and fields 

Pioneer 10 

2 Mar 72 

Jupiter flyby 

Studied Jovian system, asteroid 
belt, and heliosphere 

Pioneer 11 

5 Apr 73 

Jupiter-Saturn flyby 

Made additional studies of Jovian 
system; reached Saturn with assist 
from Jupiter’s gravity 

Pioneer 12 

20 May 78 

Venus Orbiter 

Studied Venusian atmosphere 
from orbit 

Pioneer 13 

8 Aug 78 

Venus atmosphere 
probes 

Studied Venusian atmosphere via 
4 decending probes and instru- 
mented carrier 

Helios 1 

10 Dec 74 

Heliocentric orbiter 

Studied Sun from an orbit near 
center of solar system 

Helios 2 

16 Jan 76 

Heliocentric orbiter 

Studied Sun from an orbit near 
center of solar system 

Apollo 8 

21 Dec 68 

Manned circumlunar 
flight 

Command-service module orbit of 
Moon, Returned TV pictures of 
surface 

Apollo 9 

3 Mar 69 

Manned Earth orbiter 

First flight of lunar module. 
First Apollo EVA, live TV 

Apollo 10 

18 May 69 

Manned circumlunar 
flight 

First lunar module orbit of Moon. 
First color TV from space 

Apollo 11 

16 Jul 69 

Manned lunar lander 

First landing. Obtained surface 
samples, performed EVA and 
surface experiments, live TV 

Apollo 12 

14 Nov 69 

Manned lunar lander 

Second landing. Placed experi- 
ments package on surface, 
returned parts from Surveyor 3, 
obtained surface samples 

Apollo 13 

11 Apr 70 

Manned lunar lander 

Service module oxygen tank 
rupture. Crew returned safely to 
Earth after flight around Moon 

Apollo 14 

31Jan71 

Manned lunar lander 

Third landing. Returned surface 
samples, performed experiments 

Apollo 15 

26 Jul 71 

Manned lunar lander 

Fourth landing. First use of Rover 
vehicle, obtained surface samples 

Apollo 16 

16 Apr 72 

Manned lunar lander 

Fifth landing. Obtained surface 
samples, performed experiments, 
explored with Rover vehicle 


Table 1 (contd) 


Mission 

Launch date 

Objective 

Results 

ApoUo 17 

7 Dec 72 

Manned lunar lander 

Sixth landing. Obtained surface 
samples, performed experiments, 
explored with Rover vehicle 

Viking 1 Orbiter 

20 Aug 75 

Mars Orbiter 

Studied Martian surface and 
atmosphere from orbit; relayed 
Viking Lander data to Earth 

Viking 2 Orbiter 

9 Sep 75 

Mars Orbiter 

Studied Martian surface and 
atmosphere from orbit; relayed 
Viking Lander data to Earth 

Viking 1 Lander 

20 Aug 75 

Mars Lander 

Studied surface conditions; 
analyzed soil; returned TV 
images of landscape 

Viking 2 Lander 

9 Sep 75 

Mars Lander 

Studied surface conditions; 
analyzed soil; returned TV 
images of landscape 

Voyager 1 

20 Aug 77 

Jupiter/Satum flyby 

Studied two planetary systems; 
returned TV images including 
Jupiter’s Red Spot, Saturn’s 
rings, and many satellites 

Voyager 2 

5 Sep 77 

Jupiter-Saturn flyby 

Studied two planetary systems; 
returned TV images including 
Jupiter’s Red Spot, Saturn’s 
rings, and many satellites 
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This article describes ongoing Deep Space Network support of the Helios 1 spacecraft. 
In addition, planning for a Solar Corona Faraday Rotation Experiment during the Decern’ 
ber 1981 14th Helios 1 perihelion and subsequent solar conjunction phase is detailed. 


I. Introduction 

The Helios Mission, a cooperative effort of the Federal 
Republic of West Germany and the United States of America, 
was initiated to study the properties of the Sun via spacecraft 
measurements at close range. Two Helios spacecraft were 
launched; Helios 1 on December 10, 1974 (Ref. 1) and Helios 2 
on January 15, 1976 (Ref. 2). Helios Mission support through 
1978 was regularly documented in the TDA Progress Report 
(Ref. 3). On March 21, 1980, the Helios 2 spacecraft suffered 
a failure of the traveling wave tube amplifier, and (essentially) 
no useful data could be salvaged, although periodic attempts to 
communicate with and acquire data from Helios 2 were made 
until January 8, 1981 , at which time the spacecraft transmitter 
was commanded off. On the other hand, the Helios 1 spacecraft 
continues in reasonably good health, and has been regularly in 
communication with the Deep Space Network (DSN) until the 
present time. This article describes recent DSN support of 
Helios 1 , and, in particular, describes planning for a December 
1981 Helios 1 Faraday Rotation Experiment. 

II. Helios 1 Support 

During 1981, Helios 1 was supported by the German Weil- 
heim tracking station, and by (primarily) the 26-m DSN 


stations at Goldstone, California (DSSll) and Honeysuckle 
Creek, Australia (DSS 44). During this period, important inner 
solar system science data were generated by onboard Helios 1 
science instruments. In turn, these data are being correlated 
with science measurements from other spacecraft, such as the 
Voyagers, Pioneers, and ISEE-3, to further study such topics 
as the modulation of cosmic radiation with solar cycle, the 
observation of coronal mass ejection events from the solar 
limb, and the triangulation of Types III and IV solar events 
(Ref. 4). During 1981, which is marked by the beginning phase 
of the decline of solar activity during the 21st solar cycle, 
interesting phenomena observed in Helios 1 science data 
included magnetic bubbles, short periods of rapid decline in 
solar wind intensity, and the appearance of simple ionized 
helium (Ref. 4). 

The 12th perihelion of Helios 1 occurred on Decem- 
ber 12, 1980, and the 13th on June 13, 1981. The 14th 
perihelion passage will occur on December 20, 1981, with a 
solar conjunction phase immediately - thereafter, and this 
opportunity will be utilized to perform a Faraday rotation 
radio science experiment with the Helios 1 signal and the 
DSN 64-m subnetwork. 
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III. Faraday Rotation Measurements in the 
Solar Corona 

In November 1968, the Pioneer 6 spacecraft was occulted 
by the Sun. The Pioneer 6 spacecraft transmitted a linear 
polarized carrier, which allowed the first experiment to be 
performed using a man-made signal source in determining 
Faraday rotation through the solar corona (Ref. 5). The 
quasilongitudinal approximation for Faraday rotation is given 
by (Ref. 6): 

n = QT^ J N^B^R 

where 

12 = Faraday rotation, deg 
/ = signal carrier frequency, Hz 
Q = 13548X 1Q6 
R = signal path, m 

= electron density, m“^ 

= solar longitudinal component of magnetic field, 
tesla 

Therefore, measurement of coronal-induced Faraday rotation 
provides information about both electron density and the solar 
magnetic field (Ref. 7). To support this class of experiment, 
the 64-m subnetwork was equipped with remotely rotatable, 
microwave linear feeds (Ref. 8). Closed-loop polarimeters 
(Ref. 9) were implemented and are used to automatically mea- 
sure the orientation of the received signal polarization. This 
system maximizes the received signal strength and yields precise 
signal polarization angle data. 


Research using the polarization measurement capability at 
DSS 14, 43, and 63 continues with the Helios spacecraft; 
Volland et al. (Ref. 10) and Bird et al. (Ref. 11) have reported 
Faraday rotation observations from the solar occultations of 
Helios in 1975. More recently, Dennison et al. (Ref. 12) have 
expanded on these results with a measurement of the gravita- 
tional deflection of polarized signal radiation propagating 
through the solar corona. 

IV. Helios 1 Faraday Rotation Experiment 
Planning 

In previous Helios Faraday rotation experiments (Refs. 10, 
11, 12), heavy coverage requirements of other flight projects 
prevented the Helios Project from obtaining continuous cover- 
age during the experiment period, causing sizeable gaps in the 
data and thereby somewhat hampering scientific analysis and 
interpretation of the data because of possible 180° ambiguities. 
In the upcoming 14th perihelion passage and subsequent solar 
conjunction, the geometry is such as to cause the solar conjunc- 
tion phase to be compressed into a rather short period; Figs. 1 , 
2, and 3 detail the perihelion passage and subsequent solar 
conjunction geometry. In addition, total flight project cover- 
age requirements on the 64-m subnetwork during December 
1981 are light by historical standards. The fortuitous combi- 
nation of these circumstances allowed Helios 1 to obtain 
nearly continuous scheduled coverage on the 64-m subnet- 
work from day of year (DOY) 353 at an ingoing Sun-Earth- 
Probe Angle (SEP) of 4.1 degrees, to DOY 359, at an out- 
going SEP angle of 3.8 degrees. The schedule of tracks planned 
at the 64-m subnetwork in support of this experiment are 
listed in Table 1. Demonstration passes for this experiment 
were conducted as follows: 

DSS 14: DOY 344 and 348 
DSS 43: DOY 338 and 339 
DSS 63: DOY 335 and 337 
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Table 1. Planned Helios 1 64-m subnetwork coverage 


Day 
of year, 
1981 

DSS 

Pre- 

calibration, 

Z 

Post- 

calibration, 

Z 

Acquisition 
of signal, 
Z 

Loss 

of signal, 
Z 

352 

43 

16:55 

08:40 

19:05 

08:10 

353 

63 

08:00 

16:20 

09:30 

15:50 

353 

14 

13:15 

00:00 

15:15 

23:25 

353 

43 

22:15 

08:50 

23:45 

08:20 

354 

63 

06:45 

16:25 

08:15 

15:55 

354 

14 

13:25 

00:15 

15:25 

23:45 

354 

43 

17:20 

08:50 

19:20 

08:20 

355 

63 

06:50 

16:30 

08:20 

16:00 

355 

14 

13:35 

18:30 

15:35 

18:00 

356 

63 

10:00 

16:40 

12:00 

16:10 

356 

14 

13:40 

00:35 

15:40 

00:05 

356 

43 

17:30 

09:10 

19:30 

08:40 

357 

63 

06:35 

16:30 

08:35 

16:15 

357 

14 

13:45 

00:50 

15:45 

00:20 

357 

43 

17:40 

09:20 

19:40 

08:50 

358 

63 

06:40 

16:50 

08:40 

16:20 

358 

14 

13:50 

00:50 

15:50 

00:20 

358 

43 

22:20 

09:25 

23:50 

08:55 
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Fig. 1. 1981 -1982 Helios 1 position In the ecliptic plane assuming a fixed Earth-Sun line 
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Fig. 2. 14th perihelion and subsequent solar conjunction phase of 
Helios 1 in the ecliptic plane assuming a fixed Earth-Sun line 
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Application of a dc magnetic field on an rf hydrogen gas dissociator has led to the 
development of a more efficient source of atoms for hydrogen masers. We discuss the 
influence of the applied magnetic field on the plasma and some details of the resulting 
enhancement of the dissociation efficiency. We also consider the ways that this improved 
efficiency can result in the enhancement of the operation and field life of hydrogen 
masers, as well as aid the development of spaceborne masers. 


I. Introduction 

The development of an rf hydrogen gas dissociator, effi- 
cient in both atom production and power consumption, is 
of great interest in hydrogen maser research and develop- 
ment. The significance of a low-powered source for spaceborne 
masers is self evident: it reduces the power requirement for 
the operation of the maser and aids the power economy in 
the spacecraft. In ground-based masers, an efficient source 
can contribute to a more reliable operation as well as ex- 
tended operational life. 

The details of the influence of the rf dissociator on the 
operation of the hydrogen maser have been previously dis- 
cussed (Ref. 1). Briefly, the dissociator produces atoms that, 
upon their exit through a multichannel collimator, form an 
atomic beam of hydrogen. The beam enters the inhomogen- 
eous field of a state selecting magnet which focuses atoms in a 
particular hyperfme sublevel into the entrance of a storage 
bulb within the maser high-Q cavity. The stored atoms interact 
with the microwave field coupled to the cavity, and undergo 
stimulated emission. 


The efficiency of the dissociator in atom production can 
then affect the operation of the maser in three important 
ways. First, since the state selecting magnet cannot alter the 
path of undissociated molecules which exit the collimator, a 
significant number may enter the storage bulb in the cavity. 
The presence of the molecules in the storage bulb, the number 
of which depends on dissociation efficiency, contributes to 
collisional line broadening of the hydrogen line. Second, the 
required back pressure of the hydrogen gas for production of 
a useful flux of atoms, as determined by the source efficiency, 
constitutes the load on the vacuum pumps. The vacuum pump 
saturation and the resulting decrease in their hydrogen capacity 
is a frequently encountered mode of maser failure. Finally, 
parameters of pressure and input power directly influence 
the operational life of the dissociator itself. The failure of the 
hydrogen gas dissociator is also an often encountered mode 
of maser failure. 

Motivated by the above con^sideration, we undertook to 
improve the operational efficiency of the rf hydrogen gas 
dissociator (source), based on an initial study of the rf plasma 
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processes (Ref. 1). Here we report on the development of an 
rf hydrogen source which is capable of producing atoms with 
an efficiency comparable to conventional sources, but which 
consumes a significantly lower input rf power. This source 
can also operate in a low pressure mode, where the efficiency 
in atom production is improved, while the power is kept 
constant. The improvement in the source performance is 
obtained by the application of a dc magnetic field on a con- 
ventional rf dissociator. 

In Section II of this paper we discuss theoretical back- 
ground and the relevant concepts for this atom source. Sec- 
tion III presents some details of the source and its operation. 
Finally, a summary is given in Section IV, together with an 
outline of recommended future work. 


II. Theory 

Molecular dissociation in an rf plasma is due to the colli- 
sion of molecules with electrons that move under the influence 
of the applied rf field. The rf field is coupled to the source in 
a capacitive or an inductive mode. In inductively coupled 
schemes, as employed in JPL masers, the motion of electrons 
may be regarded as circular motion about the axis of the 
oscillating magnetic field. These electrons transfer their energy 
to particles in the plasma through inelastic collisions and are 
eventually lost to the walls of the source or to recombination 
processes with positive ions. The efficiency of atom produc- 
tion in the rf source is therefore related to the energy of the 
electrons, as well as their collisional frequency with hydrogen 
molecules. The energy of the electrons is determined by the 
power input, while their collision frequency is related to 
molecular density via the pressure in the source. It is, however, 
undesirable to increase the power input or the source pressure 
in order to improve the efficiency of the atom production. An 
increase in the power, which is undesirable for spaceborne 
maser application, is responsible for generation of ions with 
higher energies. As the energetic ions collide with the glass 
wall, they cause sputtering, which then results in an increase 
in the recombination rate of hydrogen atoms on the walls. 
With a gradual degradation of the glass walls, the operational 
life of the source is degraded accordingly. As for the increase 
in the number density of molecules by increasing the pressure, 
we have already mentioned the desirability of a reduced gas 
load to the vacuum pumps, as well as a reduced molecular 
density in the storage bulb. 

Consider, however, the motion of an electron in an rf 
plasma which is superimposed by an axial dc magnetic field. 
The electron will experience a force F = e\ X B, where v is 
the electron velocity and B is the magnetic induction. This 
implies that the path of electrons would be modified through 


the action of the field on the component of velocity perpen- 
dicular to it. The result is that the electron will circulate the 
field with Larmor frequency / = eHjlTimc, while its parallel 
motion is unaffected by the field. Here H is the magnetic 
field, m is the mass of the electron, and c the velocity of 
light. The helical motion of the electron thus modified con- 
stitutes a lengthening of its path in the plasma, as compared 
to its path in the absence of the magnetic field. Furthermore, 
the circulating electron in the dc field can now absorb energy 
in resonance from the field as it moves in the direction of the 
field, when its cyclotron frequency - eH/mc is equal to 
the frequency of the oscillating field co. 

In this manner, it is possible to increase the efficiency of 
atom production in the rf dissociator by the application of 
the magnetic field. It should be mentioned that dc magnetic 
fields have been previously applied to rf plasmas used for 
production of positive ions (Ref. 2). While small dc fields do 
not directly affect the path of ions due to ions’ relatively 
large mass, they significantly modify the drift of ions in the 
plasma (Ref. 3). It is the modification of ion drift that pri- 
marily contributes to improved performance of rf ion sources 
with axial dc magnetic fields. Furthermore, because of the 
general dependence of the ionization cross-section on energy 
(Ref. 4), increased electron energy due to the applied dc field 
directly improves ion production in rf ion sources, an effect 
that is not necessarily desirable in molecular dissociation. It 
is therefore quite important to keep in mind that the influence 
of a dc magnetic field on the efficiency of atom production 
in rf dissociators is significantly different from its influence 
on ion current production in rf ion sources. 


III. The Source and Its Operation 

Because the hydrogen gas dissociator assembly of JPL 
masers is inside the vacuum package and therefore not acces- 
sible without breaking of vacuum, a dissociator test bed was 
constructed. The test bed consists of a vacuum canister for 
connecting an ion pump to a source mating flange. A pyrex 
glass bulb was designed and constructed with dimensions 
identical to those of conventional sources used in JPL masers, 
but suitably modified for use in the test bed. The modification 
consists of a glass ring with an 0-ring groove externally attached 
to the wall of the source containing the exit collimating holes. 
The source may then be positioned on the flange with the 
0-ring providing the necessary vacuum seal. The rf exciter 
coil assembly and hydrogen gas feed lines provide the neces- 
sary support for the glass bulb. A picture of the test bed and 
the source is given in Fig. 1. 

For the generation of the dc magnetic field a solenoid was 
constructed from an insulated copper wire wound on a copper 
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cylinder with an inside diameter one inch larger than the 
diameter of the dissociator bulb. The solenoid was placed on 
the mating flange in a coaxial configuration with the source. 
This particular solenoid provided a maximum field of 145 gauss 
on its axis with a 2.0-A current input. The solenoid pro- 
duced a field which was linearly dependent on current within 
the range of zero to 2.0-A, the maximum current which the 
power supply produced. 

The dissociation efficiency of the source was determined 
by measuring the changes in the pressure observed when the 
plasma was on or off. The variation of the pressure was ob- 
tained by noting the current in the ion pump power supply. 
While this method of efficiency determination for the dissocia- 
tion lacks a high absolute accuracy, it does provide a satisfac- 
tory means for a relative determination of the dissociator 
efficiency. 

The source was initially operated under the condition of 
power and pressure similar to those employed in the maser. 
A back pressure of hydrogen gas produced a pressure of 
200 m torr in the dissociator; an rf power input of 3 W pro- 
duced an efficiency measured to be 5%. As the input power 
was decreased, the efficiency of the source deteriorated, and 
the plasma was not sustained if the power was reduced below 
2 W. The source was left to operate under normal operating 
conditions for a period of a few days to allow its performance 
to stabilize. After this period, the solenoid power supply was 
turned on and the solenoid current was varied between zero 
and 2 A. As the magnetic field was increased, the light inten- 
sity of the plasma increased significantly. Furthermore, for 
specific values of the field a change in the plasma mode was 
observed. A measurement of the efficiency of the source at 
this level of pressure and power was performed with the dc 
magnet on. The source exhibited an increase and a decrease 
in efficiency depending on the magnitude of the applied mag- 
netic field. The amount of variation in the efficiency ranged 
from a 10% increase to a 10% decrease of the efficiency 
measured without the applied magnetic field. 

The efficiency measurements were continued as the input 
rf power was decreased. First it was observed that the applica- 
tion of the magnetic field decreased the value of the mini- 
mum power required for maintaining a plasma, depending on 
the current in the solenoid. For a current of 0.4 A corre- 
sponding to a field of approximately 40 gauss, it was possible 
to reduce the power to about 0.3 W. With these parameters, 
the efficiency of the source was determined to be 5%, the 
same value as the efficiency of a conventional source operat- 
ing with 3 W of input rf power. The dc magnetic field then 
resulted in a tenfold decrease in the power consumption. 


This result could be predicted in the following way. At a 
given value of pressure, a specific level of input power yields 
the optimum efficiency. A decrease in this power level would 
decrease the efficiency, while an increase in the level does not 
necessarily lead to increased efficiency because of the differ- 
ence in the dissociation energy of the hydrogen molecule 
(4.8 eV) and its ionization potential (15.8 eV). While the 
dissociation and ionization of the molecule constitute com- 
peting processes in the rf plasma, an increase in the average 
energy of electrons, with a Maxwellian energy distribution 
in the plasma, contributes preferentially to the ionization 
process, rather than the desired dissociation of the molecules. 
Therefore increasing power from the optimal level may in fact 
decrease the dissociation efficiency by contributing to the 
ionization process, via the increase in the energy of elec- 
trons. Conversely, when the power is decreased from the 
optimum level, an appropriate magnetic field strength can 
enhance the dissociation by allowing the electrons circulating 
with cyclotron frequency to absorb energy in resonance with 
the field and regain the condition of optimal energy for 
maximum dissociation efficiency. 

Based on this consideration, we decided to investigate the 
influence of an inhomogeneous field on the efficiency of the 
source. We reasoned that if the energetic electrons in the high 
energy tail of the Maxwellian distribution could be swept 
away from the region near the exit of the source, a larger 
number of electronic collisions with the molecules will result 
in dissociation in that region, rather than the competing ioni- 
zation. A test of this reasoning was demonstrated by distorting 
the homogeneity of the magnetic field through positioning the 
axis of the solenoid away from the axis of the source. When 
the solenoid was positioned with its axis about 1/2 inch away 
from the source axis, an increase in the dissociation was 
observed for a 200-m torr pressure of hydrogen. With a de- 
crease in the pressure to 80 m torr, the measured dissociation 
efficiency of the source increased from a 7% level to 12.8% 
when a current of 0.75 A was applied to the solenoid; this 
is a 70% increase of the dissociation efficiency- The power 
level in both cases was 2.5 W. 


While the configuration of the inhomogeneous field was 
somewhat arbitrary, determined by the constraint of the 
existing setup, it clearly demonstrated the potential value of 
this approach in designing a highly efficient source. It should 
be mentioned, however, that such a design should include 
consideration of influence of the inhomogeneous magnetic 
field on the production of secondary electrons released by 
particle impact with the walls. These electrons participate in 
plasma processes and influence the parameters of the rf 
plasma. . 
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IV. Summary and Conclusion 

The development of a dc magnetic field to the rf dissocia- 
tion has yielded the following results: 

(1) For a specific value of the magnetic field, it is possible 
to maintain the efficiency of the source and yet 
decrease the rf power input by an order of magnitude. 
This conclusion is of significance to both spaceborne 
maser development and improved reliability of ground- 
based masers through the expected improvement in 
the source life. 

(2) It is possible to obtain good efficiency for the source 
and, at the same time, reduce the hydrogen pressure 
with an applied dc magnetic field. This finding implies 
a means for the reduction of the gas load to the maser 
vacuum pumps, thereby increasing their operational 
life, and reducing the frequency of maser breakdowns 
resulting from pump failures. 

(3) The results with the inhomogeneous field indicate that 
it is possible to design a magnetic field configuration 
to sweep energetic electrons away from the region near 
the source exit. This would increase the source effi- 
ciency by preferential dissociation of the molecules, 
as compared to the competing process of molecular 
ionization. 

(4) The application of a magnetic field on the rf plasma 
results in a significant increase of the light output 
from the plasma. While this effect is not of interest 


to maser application (except for reaffirming that there 
is no direct relationship between the efficiency of the 
source and the intensity of the emitted plasma light) 
it may be of significance in other applications. In 
particular, certain light sources such as mercury ion 
light source, and rubidium light source, employ rf 
power for excitation. An application of a dc magnetic 
field can improve the performance of those sources 
and extend their operational life. 

The investigation of the performance of the source was 
determined in a test bed and by obtaining the dissociation 
efficiency through the measurement of pressure variations. It 
is obvious that it would be highly desirable, and necessary, 
to perform the tests using a hydrogen maser. In this way, the 
extent of improvement obtained with the new source in 
relation to the maser vacuum system, power output, and line 
Q may be directly determined. Efforts are presently underway 
to design and construct a source unit which may be mounted 
externally to our existing maser test assembly. We intend to 
optimize the configuration and the strength of the magnetic 
field for a desirable mode of source and maser operation. 
Once this is accomplished, we will substitute equivalent 
permanent magnets for the solenoid, thereby eliminating the 
need for a solenoid power source. Finally, we intend to test 
the influence of the magnetic field on the light output of 
mercury lamps, which will be utilized in a trapped mercury ion 
frequency source. The results of this investigation will be 
reported in a later paper. 
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The development of a multifrequency upconverter-maser system for the DSN has 
created the need to develop a closed-cycle refrigerator (CCR) capable of providing more 
than 3 watts of refrigeration capacity at 4,5 K, In addition, operating concerns such as 
the high cost of electrical power consumption and the loss of maser operation due to CCR 
failures require that improvements be made to increase the efficiency and reliability of 
the CCR. One refrigeration method under consideration is the replacement of the Joule- 
Thomson expansion circuit with a magnetic refrigerator. Magnetic refrigerators can pro- 
vide potentially reliable and highly efficient refrigeration at a variety of temperature 
ranges and cooling powers. This paper summarizes the concept of magnetic refrigeration 
and provides a literature review of existing magnetic refrigerator designs which have been 
built and tested and that may also be considered as possibilities as a 4 K to 15 K magnetic 
refrigeration stage for the DSN closed-cycle refrigerator. 


I. Introduction 

Efficient long-life refrigeration techniques are being devel- 
oped for the DSN masers to improve the present closed-cycle 
refrigerator technology. One area of investigation is to find an 
alternative to the Joul e-Thomson refrigeration process in the 
DSN closed cycle refrigerators (CCRs). In the CCR, as in 
virtually all mechanical refrigerators which are capable of the 
liquefaction of helium, a Joule-Thomson (J-T) circuit is 
added to an expansion engine to serve as the final stage of 
cooling to produce the desired refrigeration at the liquid 
helium (LHe) temperatures. The J-T circuit, while having no 
moving parts or seals, depends critically on extreme high gas 
purity to obtain long-term continuous operation. In fact, a 
major reason for present CCR failures is due to gas impurities 
which plug the J-T circuit. Because the J-T circuit and the 
precooling expansion engine circuit are essentially independent 


of each other, it is important to determine whether other 
refrigeration cycles can be used in place of the J-T circuit 
which will match or surpass the J-T circuit in terms of effi- 
ciency, reliability, temperature span achievable, and ease of 
construction and operation. One such refrigeration mode 
under investigation is the cooling by adiabatic demagnetiza- 
tion. 

Adiabatic demagnetization of a paramagnetic salt (or mag- 
netic refrigeration) was chosen as the possible replacement 
for the J-T circuit due to the following potential advantages 
over the J-T circuit. 

(1) High efficiencies are attainable (more than 70% of 
Carnot efficiency in comparison to 20-25% for the J-T 
expansion stage as measured for a DSN 1-watt CCR 
(Ref. 1)). 
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(2) The refrigeration capacity per unit volume of working 
material is large as high density solids rather than the 
low density gases are used. 

(3) Potential high reliability can be expected since only a 
few slow moving parts are required. 

(4) The paramagnetic material is relatively inexpensive 
(e.g., gadolinium gallium garnet costs about $225/lb). 

There are several disadvantages which must also be considered. 

(1) The need for a high field superconducting magnet 
that must be maintained at liquid helium temperatures. 

(2) Stray magnetic fields which will require shielding to 
protect the maser. 

(3) The use of low temperature seals for moving parts. 

The development work on cyclic magnetic refrigeration to 
date has shown great promise for its application as a viable 
means of providing efficient and reliable refrigeration. Their 
proposed uses include: 

(1) Laboratory material studies in the /ukelvin range 
through adiabatic demagnetization of the salt pill’s 
nuclear spins (Ref. 2.). 

(2) Spacecraft operations to cool detectors and instru- 
ments (Refs. 3, 4). 

(3) Large-scale refrigeration for superconductivity appli- 
cation in power generation and transmission (Ref. 5). 

(4) Complete magnetic refrigeration systems operating 
between 300 K and 20 K (or 4 K) for the production 
of liquid H 2 (or liquid He) (Ref. 6). 

It is interesting to note that none of the material refrigera- 
tors developed to date have been used to produce 4 K refriger- 
ation. The major development work thus far has centered 
around room temperature devices where the use of super- 
conducting magnets could be avoided, or below 4 K where a 
liquid helium bath could be used to submerse the supercon- 
ducting magnet and also be used for the high-temperature heat 
reservoir. 

Van Geuns (Ref. 7) first proposed the idea of a 4-15 K 
magnetic refrigeration stage as early as 1966 to find a more 
efficient alternative to the J-T circuit of a closed-cycle refrig- 
erator. Although the results of the analysis of the proposed 
refrigerator indicated that the prospects for the magnetic 
refrigerator were promising, no follow-up experimental work 
has been reported in the literature to our knowledge. Since 
that time there have been several studies of paramagnetic 
materials for use in a 2-20 K temperature range (Refs. 8, 9), 


and there have been more proposed designs for a 4-15 K mag- 
netic refrigerator (Refs. 10-12), but to date no experimental 
devices have been developed. 

In the next section, the concept of magnetic refrigeration 
will be described and possible thermodynamic cycles available 
for a magnetic refrigeration stage will be introduced. A brief 
discussion will be made of several magnetic refrigerators that 
have been reported on in the literature and how they may be 
used as a 4-20 K refrigerator stage for use in a DSN CCR. We 
will also introduce some of the design constraints that must be 
considered if a magnetic refrigeration stage is to be used. 


II. Principles of Magnetic Refrigeration 

Magnetic cooling is a consequence of the fact that at a 
fixed temperature the entropy of a system of magnetic dipoles 
is lowered by the application of a magnetic field. In the 
absence of an external magnetic field the magnetic moments 
of the electron spin system are randomly oriented, with all 
orientations equally probable. When a magnetic field is ap- 
plied, the dipole moments tend to align parallel to the field, 
thereby increasing order in the system, or equivalently, de- 
creasing its entropy. In doing so the moments go into lower 
energy levels, giving off excess energy to the crystal lattice 
of the paramagnetic material in the form of a heat of mag- 
netization. This heat energy is removed from the material by 
holding it in contact to a constant-temperature reservoir. Now 
by thermally isolating the material and removing the field 
isentropically, the dipoles again become randomly oriented. 
But in doing so, energy must be drawn out of the material’s 
lattice, thereby lowering the temperature of the material. 

The entropy-temperature diagram (or S, T-diagram) of a 
paramagnetic material is shown in Fig. 1. The entropy curves 
at several constant magnetic fields have been shown. The 
temperature Tq is defined as the temperature at which the 
zero magnetic field entropy per mole of material is equal to 
In (2/ + 1) (see Fig. 1). At temperatures above 7^, the 7^ 
contribution to the entropy by the lattice vibrations becomes 
appreciable. Below 7 q the molar entropy depends little on 
the lattice entropy and thus remains fairly temperature- 
independent until the interaction energy £* of the ions becomes 
an appreciable fraction of the average thermal energy kT. 
Spontaneous ordering of die dipoles then occurs and this 
lowers the entropy. This can be characterized by a magnetic 
ordering temperature Here the specific heat of 

the material also shows a pronounced maximum (since it 
varies as T • dS/dT), Between ^ and 7 q the entropy of the 
spin system depends most strongly on the magnetic field, and 
in this region the adiabatic demagnetization method is most 
effective. 
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III. Thermodynamic Cycles 

The concept of cyclic magnetic refrigeration can best be 
described with the aid of Fig. 2. In phase 1 of the cycle, 
thermal switch S2 is open and thermal switch SI is closed so 
that the working material M is in good thermal contact with 
the heat reservoir R during the application of a magnetic 
field. This permits the heat of magnetization Q- TAS to be 
removed from the material while the temperature of the 
material Tm remains nearly constant at the reservoir tem- 
perature Tf^. In phase 2, both thermal switches are opened to 
thermally isolate the working material and the field is reduced 
until Tnj is lowered to the heat source temperature T^. In 
phase 3 switch S2 is closed to permit thermal contact between 
the working material and the heat sources while the field is 
reduced to zero. Heat energy is extracted from the source 
while the final demagnetization is lowering This process 
establishes a somewhat lower source temperature T^. Phase 4 
represents the last step in the magnetic cycle. During this 
phase, SI and S2 are opened to isolate the material and the 
field is increased until is increased to Tj^. The magnetic 
cycle can be repeated by returning to phase 1 of the cycle 
and again applying the full field to the working material. 
The temperature Tg will continue to fall with each cycle of 
the refrigerator until the net refrigeration of the working 
material is balanced out by the heat leak to the source. 

The magnetic refrigerator in the above example operated 
on the basis of the Carnot cycle (two isothermal steps and two 
adiabatic steps). It is an easy cycle to execute but has a limited 
temperature span while requiring a large field change. Other 
magnetic analogs to the gas thermodynamic cycles exist which 
are more difficult to carry out, but they are capable oflarger 
temperature spans while requiring smaller magnetic field 
changes. They are: the Stirling cycle (Ref. 6), which consists 
of two isothermal and two isomagnetization steps; the Ericsson 
cycle (Ref. 6), which consists of two isothermal and two 
isofield steps; and the Brayton cycle (Ref. 6), which consists 
of two adiabatic and two isofield steps. The four cycles are 
shown in Fig. 3, but they do not represent the only refriger- 
ation cycles possible. The operation of a magnetic refrigerator 
based on one of these thermodynamic cycles will depend 
largely on the design of the refrigerator, the ease of implemen- 
tation and the environment in which the refrigerator will 
work. 


IV. Approaches to Actual Magnetic 
Refrigeration Devices 

A. Two-Switch Devices 

It was the Nobel prize winning work of Giauque and 
MacDougall (Ref. 13) in 1933 which demonstrated the idea 


of magnetic refrigeration to cool from 3.5 to 0.5 K. Subse- 
quent magnetic refrigeration work was able to produce tem- 
peratures to below 1 mK. Until recently, the method of 
adiabatic demagnetization was limited to so-called “one-shot 
operations” where the paramagnetic material was isothermally 
magnetized and followed by an adiabatic demagnetization to 
obtain temperatures below 1 mK for limited periods of time. 

In 1954 Heer et aL (Ref. 14) developed the first magnetic 
refrigerator to operate in a continuous mode to provide a 
fairly constant low temperature. The device used supercon- 
ducting thermal valves (as diagrammed in Fig. 1) and a high 
heat capacity material for the low-temperature reservoir to 
smooth out the low-temperature fluctuations. This refrigerator 
was capable of producing 7 juW of refrigeration at 0.26 K. 
Arthur D. Little, Inc. commercially marketed a few magnetic 
refrigerators modeled after this design (Ref. 15) until ^He 
became commercially available in sufficient quantities to use 
as a refrigerant in this temperature range. 

The next refinement was by Zimmerman et al. (Ref. 16) 
in 1962, by using a superconducting solenoid. This refrigera- 
tor was capable of producing 100 juW of refrigeration at 
0.26 K. More recently Rosenblum et al. (Ref. 17) used a 
similar device to produce 0.1 /uW of refrigeration at the much 
lower temperature of 10 mK. All of these refrigerators were 
low power devices and limited to temperatures below 1 K 
because of the superconducting switches used. Steyert (Ref. 
12) recently proposed a 2-10 K refrigerator which uses 
magnetoresistive switches (Ref. 18) with large switching 
ratios to produce several watts of refrigeration at 2 K using a 
Carnot cycle. No such refrigerator has been built to date. 


B. Reciprocating Devices 

The reciprocating magnetic refrigerator requires the com- 
plex movement of two of the following components: the 
magnet, the magnetic material, or the regenerator fluid. 
Described simply, the working material passes through the 
regenerative column of fluid in which a temperature gradient 
has been established after numerous cycles. The regenerative 
fluid does the work of passing the heat from the heat source to 
the heat sink. The mixing of the fluid across the temperature 
gradient by the movement of the material is the key area of 
concern for the design of the device. An illustration of one 
such refrigerator (Ref. 19) is shown in Fig. 4. 

Three reciprocating magnetic refrigerators have been built 
and tested since Van Geuns (Ref. 7) first proposed the idea of 
a 4-15 K reciprocating magnetic refrigeration stage to replace 
the J-T circuit of his closed cycle refrigerator. Brown (Ref. 19) 
utilized a magnetic Stirling cycle to demonstrate a 47 K tem- 
perature span between 272 and 319 K. He used gadolinium 


31 


metal as the working material which passed through a regen- 
erative column of water and ethyl alcohol. Barclay et al 
(Ref. 20) used a salt pill of compressed Gd^ ( 504)3 * ^^ 2 ^ 
pump heat from 2 to 4 K. Using a magnetic Ericsson cycle at 
a 1/60 Hz rate, they were able to develop 52 mW of refrigera- 
tion capacity at 2 K in their initial tests. Delpuech et al. 
(Ref. 21) built a 2 to 4 K double-acting reciprocating device 
in which two identical Gd 2 ( 504)3 niagnetic elements on a 
reciprocating drive shaft alternately entered and left the 
magnetic field. Incorporating this arrangement with a magnetic 
Ericsson cycle, they were able to attain 570 mW of refrigera- 
tion capacity at 2.1 K at a 0.3-Hz cycle rate. In later experi- 
ments using Gd 2 Ga 5 0 j 2 magnetic elements, they were able to 
achieve 1 .2 W at 1.8 K with a 0.95 Hz rate (Ref. 22), 

C. Rotating Wheel Devices 

The wheel-concept magnetic refrigeration device uses the 
magnetic material to form the rim of the wheel. It has a simple 
axial drive but requires a more complicated magnetic field 
profile since one portion of the wheel must be in a large field 
while the opposite side must be in a low magnetic field. How- 
ever, because one portion of the rim is constantly rotating into 
or out of the magnetic field, the wheel design does provide 
continuous refrigeration. Figure 5 illustrates the concepts of 
both a Carnot-cycle and a Stirling-cycle wheel design. 

Steyert and co workers (Refs. 5, 24) have developed several 
rotating wheel designs which have operated at either room 
temperature or between 2 and 4 K. The prototype low- 
temperature Carnot-cycle wheel design (Ref. 24) used Gd 2 - 
(SO^), • 8 H 2 O as the working material and was capable of 
pumping 0.5 W of heat from 2.75 K to the 4.2 K heat sink. 
The wheel rotated at 0.77 Hz. The capacity and source tem- 
perature varied with the rotation rate of the wheel. 

As a feasibility study for the use of low-temperature mag- 
netic refrigeration technology to provide refrigeration for 
superconducting applications in power generation and trans- 
mission, Barclay and Steyert (Ref. 5) built a room-temperature 
wheel with an objective of understanding the device in an 
experimentally convenient temperature range. Their proto- 
type gadolinium metal wheel operated through a 1.2 tesla 
field change and provided 0.5 kW of continuous refrigeration 
over a 7 K temperature span while rotating at a 0.5 Hz rate. 
This Stirling-cycle wheel device used water as the regenerative 
fluid. 

D. Rotating Anisotropic Crystal 

Barclay (Ref. 11) proposed this magnetic refrigerator 
design based on the anisotropic nature of certain single-crystal 
paramagnetic materials. The . anisotropic material, such as 
DyP 04 , is rotated in a constant magnetic field causing a 


drastic change in the material’s magnetic entropy. The temper- 
ature changes associated with the changes in the paramagnetic 
ion polarization are used directly to heat or cool a heat ex- 
change fluid. Each 360° rotation of the crystal represents 
a double execution of the thermodynamic cycle. Barclay 
describes a Carnot-cycle 4 to 20 K rotational refrigeration 
device that is capable of producing 1 W of refrigeration. A 
conceptual drawing of the device is shown in Fig. 6 . 

E. Active Magnetic Regenerator Device 

The active magnetic regenerator is a device composed of 
layers of one or more magnetic materials having different 
Curie temperatures (depending on the temperature span 
desired) that are thermodynamically cycled like the ordinary 
lead regenerator of the CCR, except in this case the tempera- 
ture of the materials can be changed by the application or 
removal of a magnetic field. Each of the different materials 
executes a small Brayton cycle near its Curie temperature. 
Depending on the number and types of different materials 
used, a temperature span of 10-220 K may be achieved. 
This concept was proposed by Barclay and Steyert (Ref. 10) 
and is still in the development stage. A prototype room- 
temperature device has been built and tested (Ref. 5) in which 
initial results show an 18 K temperature span using only 
gadolinium metal as the magnetic material. 


V. A Magnetic Refrigeration System 

The present and future requirements of the CCR system 
for the DSN are listed in Table 1 . With the development of a 
multifrequency upconverter maser system there exists the 
need to develop a new CCR capable of providing 3-4 watts 
of refrigeration power at 4.5 K, A scaled-up version of the 
present DSN IW CCR, capable of more than 3 W of refrigera- 
tion power at 4.5 K, is presently being tested. Also under 
investigation is the development of a 4-15 K magnetic refriger- 
ation stage which will replace the J-T expansion circuit to 
provide the final stage of cooling in the CCR. 

The reasons for considering a 4 to 1 5 K magnetic refrigera- 
tion stage are the advantages it has over the J-T circuit. This 
replacement eliminates the CCR failures that are a direct 
result of J-T circuit blockage due to gas impurities. The 
refrigeration power per unit volume of paramagnetic material 
is large (since we are dealing with a solid as the working sub- 
stance rather than a gas, and the density of the solid is more 
than an order of magnitude greater than that of the com- 
pressed gas). This should facilitate the use of a small and 
compact design. The magnetic cycle itself can operate at 
nearly 100 % efficiency; when the losses associated with the 
magnetic stage are included, the Carnot efficiency of the 
magnetic stage should still be on the order of 70% or more. 
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The magnetic refrigerator stage will require some moving parts, 
but as they will be moving at slow speeds the overall CCR 
reliability should not be impaired. 

In order to incorporate a 4-15 K magnetic refrigeration 
stage into a closed-cycle refrigerator, the refrigerator compo- 
nents and their constraints must first be identified and fully 
understood. The major system components of the magnetic 
refrigerator are listed in Table 2. While the exact magnetic 
thermodynamic cycle to be used will play a part as to what the 
final design configuration of each component will be, it is 
nevertheless important to identify the main components and 
describe their general details. In what follows, we will discuss 
the pertinent points of the first four items listed in Table 2. 


A. Magnetic Working Material 

It is necessary to select the working material and the mag- 
netic cycle which will best operate over the temperature span 
desired. Through the temperature span of interest, both the 
entropy and the heat capacity of the spin system of the 
material must be large compared to those of its lattice system 
so that the low spin temperature attained during a demag- 
netization is not appreciably curtailed by the internal heat 
load of the lattice. The magnetic ordering temperature of the 
material should be a degree or two below the lowest refrigera- 
tion temperatures desired because the magnetic entropy 
changes with external field become very small below this 
point. Other important characteristics to be considered must 
be its thermal, chemical and physical properties. Finally, 
cost, availability and machinability are also important con- 
siderations, 

A study of paramagentic materials suitable for magnetic 
refrigeration in the 2-20 K range has been conducted by Bar- 
clay and Steyert (Ref. 9) for JPL as part of the overall R&,D 
program on magnetic refrigeration at the Los Alamos Scienti- 
fic Laboratory. The study was restricted to gadolinium com- 
pounds. In particular gadolinium gallium garnet (Gd 2 Ga 5 0|2 
or GGG) was shown to be an excellent choice as the working 
material in the 4-15 K range because of its high thermal 
conductivity, its very low lattice heat capacity and its low 
ordering temperature ('^O.S K). The material is the standard 
substrate material used in magnetic bubble domain tech- 
nology and therefore is readily available. Grown as a single 
crystal boule it is also very machinable. This material is pre- 
sently being examined at JPL and results will be forthcoming. 


B. Heat Exchange Mechanism 

During the magnetic cycle, heat is periodically exchanged 
between the working material and the heat sink or heat 


source. There are two distinct ways in which this heat ex- 
change can take place: through a thermal switch such as is 
needed in the isothermal stage of a Carnot cycle or through a 
regenerative fluid flow as would be used in the Stirling cycle, 
for example. Thermal switches require large switching ratios 
if they are to work efficiently. Two likely candidates would be 
magnetoresistive switches (Ref. 18), such as single crystals of 
beryllium, with switching ratios of 1000 at temperatures as 
high at 15K, or helium gas gap switches which are actuated 
quickly by electrically heating charcoal to change the gas 
pressures. Of main importance is that no regenerative fluid is 
used since the field is changed by varying the magnet. How- 
ever, this may cause serious problems for a nearby maser 
unless adequate shielding is provided. As another heat exchange 
mechanism, a heat exchange fluid such as high pressure helium 
gas can also be used to couple the working material to the 
hot or cold sources and to effect the regeneration. Fluid 
pumps will be required to operate at low temperatures to 
force the circulation of the fluid. It is perhaps better still to 
use the working material as the regenerator as in the active 
magnetic regenerator device. The magnetic material could be 
shaped in small spheres and have the helium forced through 
the regenerator bed to attain the necessary heat transfer. 


C. Magnetic Material Drive Mechanism 

The work for the magnetic refrigeration stage is provided 
by the drive motor. For a reciprocating design a force compen- 
sation mechanism must be utilized to pull the material in 
and out of the field because of the large forces involved, A 
more complicated cam and drive shaft are required if both 
the magnet and the magnetic material are to move. For both 
the rotating wheel and the rotating anisotropic crystal a 
simple axial motion of the drive shaft is used. The wheel 
provides its own force compensation by its continuous geom- 
etry; therefore the force from the drive shaft’s torque will be 
smaller than the force required for the reciprocating drive 
units. 


D. Superconducting Magnet 

The biggest disadvantage in a magnetic refrigeration stage 
for the CCR is the need for a high field superconducting 
magnet. Stray fields from this magnet would pose a serious 
problem for maser operation unless effective shielding was 
provided. The shielding should be adequately accomplished 
through the use of high permeability or superconducting 
materials. 

Another problem to be overcome is how to cool the mag- 
net initially. The magnet can be cooled by an initial LHe 
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transfer but this method is not very feasible in light of the 
CCR being installed on the antennas. The magnet will most 
likely be precooled to 12-15 K by the two-stage expansion 
engine as the J-T circuit is now. From there, two options 
are possible: to use a Simon liquefier (expansion of stored 
high pressure helium gas that has been precooled to under 
15 K) to produce sufficient LHe to cool the magnet; or to 
start the magnetic cycle at low fields at 12-15 K and have 
the magnetic stage cool itself down. The latter choice re- 
quires a high superconducting wire such as Nb 3 Sn. In 
either situation, once the magnetic refrigerator is operating, 
the small parasitic heat load required by the magnet operating 
in the persistent mode will require only a small fraction of 
the total large refrigeration capacity of the magnetic refrigera- 
tion stage. 


VI. Conclusions 

The idea of using magnetic refrigeration to provide the 4 K 
refrigeration for the DSN maser has been introduced. The 
main feature of the magnetic refrigerator is its potential high 
efficiency (more than 70% of Carnot efficiency can be ex- 
pected). Additionally, because it uses solid materials it is 
inherently orientation-independent. This, and other design 
constraints laid on the refrigerator by the operating maser 
system on the antenna have been listed. Several different 
refrigerator devices have been presented as possible designs 
for a 4-15 K magnetic refrigeration stage. The question to be 
answered will be which type of magnetic refrigeration device 
using which thermodynamic cycles will prove to be optimal 
for our present and future refrigeration requirements. 
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Table 1 . CCR System requirements for maser cooling 


Reliable 

Efficient 

Multiyear lifetime 
Unattended operation 
Rapid cooldowns 
14 watt cooling capacity 
Compact 

Magnetic field isolation of maser package 

Low microphonics 

mK temperature stability 

Orientation independence 

Continued operation during power failures 


Table 2. Major components of a magnetic refrigerator system 


Solid magnetic working material 
Heat exchange mechanism 
Magnetic material drive mechanism 
Superconducting magnet 
Dewar 

15 K refrigerator precooler 




<I> Tq 


Fig. 1 . S, T-diagram for a typical paramagnetic material showing the 
entropy vs temperature curves for several applied magnetic fields. 
Process 1-^2 shows an isentropic (adiabatic) demagnetization 
which decreases the temperature of the material. Process 2 -> 3 
shows an isothermal demagnetization whereby the material 
absorbs heat from low-temperature heat source. Temperatures <l> 
and Tq are defined in the text 



HEAT SOURCE (S) 


Fig. 2. Schematic diagram of a conventional magnetic refrigerator. 
The cyclic operation of this refrigerator Is discussed in the text 


CARNOT 


BRAYTON 



T T 


Fig. 3. Four magnetic cycles that can be used to make a cyclic 
magnetic refrigerator. The processes of each cycle are mentioned in 
the text (taken from Ref. 6) 


COOLING 



Fig. 4. Schematic diagram of magnetic Stirling cycle. The pro- 
cesses are (a) Isothermal magnetization, (b) isofield cooling In 
regenerator, (c) Isothermal demagnetization, and (d) zero (or low) 
field heating in regenerator (taken from Ref. 19) 
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HIGH FIELD 
REGION — 


HIGH FIELD 
REGION — 



Fig. 5. Schematics of (a) Stirling-cycle wheel and (b) Carnot-cycle wheel refrigerators. The paramagnetic material forming the rim of the 
wheels rotates in a counterclockwise direction. The fluid flow Is forced by a pump and its direction of motion Is indicated by the arrows. 
The high magnetic field region is Indicated by the shaded area. Radially across the wheel Is the low field region. ^ represents the 
inherent temperature change of the working material upon entering and leaving the field (taken from Ref. 23) 



Fig. 6. Side and end views of a wheel of DyP 04 single crystals 
inside a superconducting magnet. The drive shaft and ducts are 
also shown (taken from Ref. 11) 
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Recently, there has been a substantial effort to increase the amount of data that can 
be received from outer planet missions by coherently combining signals from ground 
antennas in such a way as to increase the total effective aperture of the receiving system. 
However, as these signals become weaker, the baseband arraying technique in current use 
degrades somewhat due to carrier jitter. One solution to this problem is Sideband-Aided 
Receiver Arraying (SARA). In SARA, sidebands demodulated to baseband in a master 
receiver at the largest antenna are used to allow slave receivers in the other antennas to 
track the sideband power in the signal rather than the carrier power. The already existing 
receivers can be used in the slaves to track and demodulate the signals in either a residual 
carrier or a suppressed carrier environment. The resultant baseband signals from all the 
antennas can then be combined using existing baseband combining equipment. Computer 
simulations of SARA show increases in throughput (measured in data bits per second) 
over baseband-only combining of 17% at Voyager 2 Uranus encounter and 31%o at 
Neptune for a four-element antenna array and (7, 1(2) convolutional coding. 


I. Introduction 

Many different antenna arraying systems have been pro- 
posed for receiving increased telemetry rates from outer planet 
missions. The performance of two of these schemes, “base- 
band-only combining” (BBO) and “baseband combining with 
combined carrier referencing” (CCR), have been studied by 
the authors (Ref. 1). BBO (Refs. 2, 3) is the scheme that is 
currently used by the DSN for Voyager. In BBO each receiver 
demodulates its own carrier using a conventional phase locked 
loop (PPL). CCR schemes involve summing the carrier signals 
from the arrayed antennas in phase so that the array behaves 
like a single antenna with aperture equal to the sum of the 
apertures of the individual array elements. Two different 
implementations of CCR have been studied in detail. These are 
the “master-slave system” (Refs. 4, 5) and “virtual center 


arraying” (Ref. 6). All of the above-mentioned arraying 
schemes require the presence of a residual carrier for tracking. 

Sideband-Aided Receiver Arraying (SARA) is an arraying 
system that offers increased performance over BBO and CCR 
at high modulation indices (such as those proposed for Voy- 
ager 2 Uranus and Neptune encounters) and the ability to 
operate in a suppressed carrier environment. Suppressed carrier 
transmission, in which the modulation index is equal to 90°, 
offers a secant^(0) improvement in bit energy-to-noise ratio 
over residual carrier transmission with a modulation index of 
Q. For Voyager 2 Uranus and Neptune encounters, suppressing 
the carrier will yield a 0.25-dB improvement in EJN^. 

SARA is a master-slave type system. Either a conventional 
tracking loop or suppressed carrier tracking loop may be used 
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in the master channel to demodulate the carrier and produce 
good data and subcarrier estimates. A suppressed carrier loop 
will yield better performance in most cases, and the concep- 
tual block diagrams described in this article are for this type of 
loop. The master channel is associated with the largest aptenna 
in the array where the signal-to-noise ratio is largest. The slave 
receivers contain conventional tracking loops that are aided by 
the subcarrier and data estimates derived in the master. These 
sidebands are relatively low frequency and low bandwidth 
signals that may be transmitted digitally over existing micro- 
wave links. Hardware must be provided in each slave for 
proper alignment of its modulation with the estimates it 
receives from the master. This aiding allows the slave receivers 
to track the sideband power in the incoming signal (rather 
than the carrier power), with losses that are lower than those 
that would be present in systems with separate suppressed 
carrier loops in each slave. The demodulated baseband signals 
from the master and all slaves are then combined in a baseband 
combiner in the same way as in BBO. 

In the following sections, SARA is described in more detail. 
Loop signal-to-noise ratios are derived for the master and slave 
receivers. Degradations due to subcarrier demodulation, sub- 
carrier alignment, and data alignment are discussed. Sec- 
tion VI discusses the performance of SARA at Voyager 2 
Uranus and Neptune encounters for two-, three-, and four- 
element arrays as determined from computer simulations. 
These results show that SARA, operating in a suppressed 
carrier environment, has less radio loss than either BBO or 
CCR. Graphs of data rate vs time for the two encounters are 
also presented, and these clearly show the increased telemetry 
capability that SARA allows. 


II. The Data-Alded Loop 

In the preferred implementation of SARA, the master 
receiver should be capable of operating with both residual and 
suppressed carrier signals. Three well-known regenerative loops 
can be employed for obtaining a carrier reference from a 
suppressed carrier signal. These are the squaring loop (Ref. 7), 
the Costas loop (Ref. 7), and the “data-aided” loop (Ref. 8). 
Moreover, these implementations provide very good carrier 
reference estimates from residual carrier signals whose carrier 
power is suppressed by more than 12 dB below the data 
sideband power (Ref. 9). In the body of this article, only the 
data-aided implementation will be considered. Comparative 
analyses of data-aided and Costas or squaring loops are given in 
Refs. 7 and 8. A block diagram of a data-aided loop appears in 
Fig. 1. The input to this loop is assumed to be a phase- 
modulated sinusoid of frequency cOq with no subcarrier. Since 
the data-aided loop tracks only the data power in the signal, it 
may be assumed thatx(r) has the form 


x(r) = sin + 7 t/ 2 c/(0 + 0] +«i(0 

= \/2^ dif) cos + 0) + « j (0 

where Pj^ is the data power of the signal, d{t) is the data in the 
form of plus and minus ones that are constant for a symbol 
time Tj, and 0 is the phase of the signal. In this analysis it is 
assumed that 0 changes slowly with time compared to c?(r). 
The noise, is assumed to be a narrowband Gaussian 

process with one-sided spectral density IV/Hz in the region 
of interest. 

The output of the VCO is 

r(0 = VT/Tj cos(cOjjr + $) 

where /Tj is a gain associated with the VCO and 0 is the loop’s 
estimate of the phase error 0. It follows that, after filtering out 
the double frequency terms, 

y{t) = s/P^ d{t) cos (0 - $) + n^{i) 

where is the gain of the upper mixer and is a white 
Gaussian process with one-sided spectral density K 2 Nq 
W/U z over the closed-loop bandwidth. The data detector 
integrates j^(f) over a symbol time T^. The sign of this integral 
defines the data estimate 

e(t) = d(t - T^). 

The probability that a correct estimate is made is 
pr{d = d) = 1/2 + 1/2 erf (y/EjN^ cos [0 - 0 ] ) 

where 

erf(x) = f exp(-r^)c?r. 

Jo 

Therefore, the mean value of e(t) is given by 

7(r) = d(t - r ) pr(d= d) - d(t - T^) pr(d i- d) 

= d(t - rp erf cos [0 - ?] ). 

The quadrature carrier estimate s{t) is given by 
s(r) = - V ^3 sin (coj,/ + $) 
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where is a gain associated with the VCO and the quadra- 
ture generator box (labeled “7 t/ 2” in the figure). After the 
lower mixer, an associated low-pass filter to attenuate the 
double frequency terms, and a one symbol time delay, 

w<0 = d{t - sin ((^ - ?) + 

where is the gain in the mixer and n^(r) is a white Gaussian 
noise process with one-sided spectral density IV/Hz. 

Notice that «2 and are uncorrelated due to the fact that 
they are in quadrature with one another. 

Under the previously stated assumption that 0 varies slowly 
with respect toJ(r), 

z{t) «»7(r) w(r) 

= n/^^ 3^4 sin (0 - ^ erf bJEjN^ cos ((^ - $)] + n^{t). 

For small phase errors <j>- $ this linearizes to 

z(0 = sfP^K^K^ [0- ^ erf(Vp^)+«/0 

from which the standard PLL theory (see Appendix A) shows 
the loop signal-to-noise ratio (loop SNR) of the data-aided 
loop to be 

where is the one-sided closed-loop bandwidth of the 
receiver. 


III. The Slaves 

In SARA, the modulation estimates obtained in the master 
are sent to each slave to aid in tracking. Figure 2 shows a 
diagram of a slave that could be used in conjunction with the 
data-aided loop described in the last section. In this implemen- 
tation, a hard-limited estimate of the data (e(r) in Fig. 2) is 
passed to the slave to derive a carrier reference. In Appen- 
dix B, the rationale for this procedure is developed. The input 
signal is again assumed to have the form 

^ (1) 

where Pds is the data power in the slave and 0^ is the phase of 
the signal. It is assumed that 0^ changes slowly with respect to 
d{t) although it is not necessarily the same as 0. The carrier 
estimate from the VCO is 


sin 

where /T is a gain associated with the VCO and is the phase 
estimate of the slave loop. After filtering out the double 
frequency components and delaying, it follows that 

r,)sin(^,- . (2) 

If e(t) is perfectly aligned, then 

"" [V£‘,/Vo ~ 

+A^3(r) 

where is a white Gaussian process independent of the 
noise in the master. If both the phase errors in the master 
and the slave are small, then this can be linearized to 

erf(^/£^) . 

Standard PLL theory (see Appendix A) then shows that the 
loop SNR in the slave is given by 

Oif Ls 

where is the one-sided noise spectral density in the slave, 
and is the one-sided closed-loop bandwidth of the slave. 

IV. Subcarrier Demodulation 

In Sections II and III, only signals without subcarriers were 
considered. If the incoming signal has a subcarrier, then a 
subcarrier tracking loop must be added to the master receiver 
in order to demodulate the subcarrier before the data can be 
detected. The subcarrier estimate must also be sent to the slaves 
along with the data estimate. Figure 3 shows one possible 
configuration for a SARA master and slave with such a sub- 
carrier loop added. The subcarrier is assumed to be a square 
wave in the following analysis. 

The subcarrier tracking loop in the master is, . itself, a 
data-aided loop similar to that described in Section II. It shares 
the data detection arm with the carrier tracking loop. The 
principal difference is that the input to the subcarrier loop is 
an amplitude-modulated square wave, 

i{t) = d(t) Sin cos (^ - ^ + n^(t) 
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where <t>, and ri 2 are as defined in Section II. 

The subcarrier frequency is assumed to be high compared 
to the channel symbol rate =1/7’^. For Voyager 2 Uranus 
and Neptune encounters, = 360 kHz while < 40 kHz. 
Recall that 0 and $ also vary slowly compared to d{t). Under 
these assumptions, the analysis of the subcarrier tracking loop 
is identical to that of the carrier tracking loop with the 
exception that only the first harmonic of the subcafrier is 
retained after mixing with the subcarrier estimate j(t). Since 
the fraction of the total square wave power in the first har- 
monic is S/tt^ , the subcarrier tracking loop SNR is given by 


P 


sc 


^O^Lsc 


erf2 (Vp^) 


where is the closed-loop bandwidth of the subcarrier 
tracking loop. 

In practice, once subcarrier and carrier lock have occurred 
in the master, the bandwidth of the subcarrier loop can be 
narrowed, since most of the doppler error in the incoming 
signal is removed by carrier demodulation. Hence, in the 
steady state, subcarrier demodulation losses should be small 
compared to carrier demodulation losses in the master 
receiver. 

In order to aid the slave receivers, the subcarrier and data 
estimates are recombined. The subcarrier estimatey(r) must be 
delayed by one symbol time before being mixed with the data 
estimate. The slave is then aided just as in Section III, with 
both data and subcarrier demodulation occurring simulta- 
neously. 


V. Alignment Losses 

Thus far it has been assumed that the modulation estimates 
from the master are perfectly aligned with the modulation in 
the slaves. In practice, there will be losses resulting from an 
imperfect alignment of these signals. This misalignment will 
result in a degradation in the loop SNR of the slaves and 
therefore will result in poorer performance in those receivers. 


If the transition probability in the symbol stream d{t) is p,, 
then the average value of the product of the data signals is 


d{t)d(t^T) = (r>,-2|r|)/(r>,) 


= 1 - 2p^R^\r\ . 

These lead to a degraded slave loop SNR of 

P , = ( V ^)0 - (1 - 

Graphs of the losses due to improper subcarrier and data 
alignment in the slaves appear in Figs. 4 and 5 respectively. A 
subcarrier frequency of 360 kHz and a data frequency of 
38.4 ksps were assumed in these plots. These are typical values 
for the Voyager 2 Uranus encounter. Notice that the align- 
ment loss is much greater in the subcarrier than in the data. 

The primary cause of misalignment in SARA will be inter- 
receiver timing problems due to link jitter. A study of the 
existing GCFIO-DSS 14 microwave link (Ref. 10) shows a 
two-way link jitter on the order of 20 nanoseconds. With 
additional alignment hardware in the slaves to take care of 
gross interreceiver time delays, the alignment degradations due 
to the link can probably be kept to a tolerable level (less than 
a dB of loop SNR). 


VI. Numerical Results 

The performance of SARA on coded telemetry was deter- 
mined by software simulation. The code used was the DSN 
standard (7, 1/2) convolutional code with Viterbi (maximum 
likelihood) decoding. Perfect subcarrier demodulation, symbol 
tracking, and modulation alignment were assumed. The only 
losses were the space loss, a 90% weather confidence attenua- 
tion, and the independent carrier phase jitter in each receiver. 

Carrier loop simulation software was used to generate the 
phase errors for each of n receivers according to the distribu- 
tion 


Suppose that the modulations of the master and slave are 
misaligned by r seconds (where |r| < l/4co^^). Then the aver- 
age value of the product of the subcarriers will not be one as in 
the case of perfect alignment. Instead it will be 


Sin Sin (co^^r + t) = (l/(2co^^) - 2 |t|)/[1/(2o^^)] 


pM-) = 


exp (p , cos 
27 t/^(pP 


(/ = 1,2,3, ... ,n) 


where 0^- is the carrier phase error in the receiver, is the 
loop SNR of that receiver and /q is the zero-order Bessel 
function. The loop SNRs for SARA master and slave receivers 
are those given in Sections II and III respectively. If a signal 
with total bit signal-to-noise ratio is received by this 
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array, then the bit SNR seen at the output of the baseband 
combiner is 

TOTSNR = EJN^ ^^w.cos^.j 

where w. is the fraction of signal amplitude seen by the 
antenna. A software Viterbi decoder was used to monitor the 
decoded bit errors made by the system. 

The performances of BBOand CCR systems were also deter- 
mined by the software. BBO was simulated in exactly the same 
way as SARA except that the were the loop SNRs of 
conventional receivers. CCR was modelled as a single antenna 
whose aperture would receive the total power incident on the 
array. 

The power-to-noise ratios for Voyager 2 Uranus encounter 
were taken from design control tables (Ref. 11). The power- 
to-noise ratios for Neptune encounter were assumed to be 
3.5 dB less than those for Uranus encounter due to the 
increased distance of the spacecraft from Earth. The loop 
SNRs of the various receivers were calculated from these 
power-to-noise ratios according to the formulas given in Ref. 1 
and in Sections II and III. A threshold loop bandwidth of 
25^0 = 30 Hz was assumed in all receivers for comparison 
purposes. 

Tables 1 and 2 show the relative performance of BBO, 
CCR, and SARA at Voyager 2 Uranus and Neptune encounters 
for two-, three-, and four-element arrays. The modulation 
indices were optimized to within 2 degrees. The numbers 
labeled “maximum data rate” represent the highest possible 
information bit rate that would allow transmission at a bit 
error rate of 5 X 10~^ . Since, in practice, Voyager has only a 
small number of available data rates, these numbers should be 
considered as a convenient measure of system performance 
rather than as actual data rates. From the tables it is clear that 
SARA performs substantially better than CCR and BBO for 
these Voyager scenarios. Data rates for SARA are about 1 1% 
higher than either of the other schemes. 

The performance of SARA was also determined at a modu- 
lation index of 76°. The maximum data rates for Uranus and 
Neptune encounters were found to be 26.5 kbps and 
1 1.7 kbps respectively. This means that data rate increases of 
2.3% for Uranus encounter and 5.4% for Neptune encounter 
over CCR are attainable with SARA without suppressing the 
carrier. 

The relative performance of the various arraying schemes is 
also shown in plots of maximum data rate vs time for the two 


Voyager planetary encounters. These plots appear in Fig. 6 for 
Uranus encounter and in Fig. 7 for Neptune encounter. (In the 
case of Neptune encounter, the time axis should be read as 
hours relative to some unspecified zero hour.) Horizontal 
dotted lines represent the data rates that are actually available 
on Voyager. One may determine the number of hours of 
transmission at an average bit error rate of 5 X 10'^ from 
these graphs for each of the arraying schemes, data rates, and 
array configurations. This is done by observing how much of 
the particular data rate curve is above the data rate line in 
question. For example, with a four-element array at Uranus 
encounter, SARA provides about 13.8 hours of viewing at a 
data rate of 19.2 kbps, while CCR allows only about 12.3 
hours and BBO only 11.8 hours. This represents an improve- 
ment in throughput of 12% over CCR and 17% over BBO. At 
Neptune, the corresponding improvements at 8.4 kbps are 12% 
over CCR and 31% over BBO. 

VII. Summary and Conclusions 

The numerical results presented in Section VI show that 
SARA offers improvements in throughput for Voyager 2 
Uranus and Neptune encounters. The improvements are due 
primarily to two effects. The first is the ability of SARA to 
work in a suppressed carrier environment. For example, while 
the optimum modulation index for CCR at Uranus encounter 
is about 76°, SARA can operate at a modulation index of 90°. 
Consequently, SARA can see 0.26 dB (6%) more data power 
than CCR. The second effect is the smaller radio loss of SARA 
as compared to the other two systems. The radio loss is lower 
because the sideband aiding in both the master and the slaves 
greatly increases the receiver loop SNRs over conventional 
phase lock loops. A greater loop SNR means smaller phase 
jitter and better tracking. 

Another advantage to SARA is that it can be made com- 
patible with existing DSN hardware. The slave receivers can be 
existing Block III or Block IV receivers with a mixer added for 
the sideband aiding. The master can be built as described in 
the text, or it can be made out of an existing receiver by 
taking the modulation estimates out of the Subcarrier Demod- 
ulation Assembly (SDA) and the Symbol Synchronizer Assem- 
bly (SSA). The modulation can be sent to the slave receivers 
over existing microwave links, and they can be aligned with 
the modulations in the slaves using reprogrammed Real Time 
Combiners (RTCs). One possible DSN implementation of 
SARA is shown in Fig. 8. 

A number of issues associated with the SARA concept 
remain to be addressed. Optimal loop filter parameters for 
specific situations such as the Voyager planetary encounters 
should be derived. The question of tracking loop acquisition is 
also very important. The sideband aiding in the slaves should 
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help them acquire in a low SNR environment where a conven- 
tional PLL might not lock at all. The master has three coupled 
loops: the carrier loop, subcarrier loop, and the data loop 
(included in the data detection function). An efficient lockup 


sequence should be established. It may be necessary to change 
filter parameters for the acquisition process. Also, a device for 
detecting when lock has occurred in each loop might be 
needed. 
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Table 1. Performance of various arraying schemes at Voyager 2 
Uranus encounter (0 hour GMT, Day 34, 19sd; BER = 0.005, 2B, « = 
30 Hz) 


Array configuration 


Optimal 
modulation 
index, deg 

Radio 

loss, 

dB 

Maximum 
data rate, 
kbps 

64m-34m(T/R): 


BBO 

74 

0.32 

16.1 


CCR 

74 

0.25 

16.4 


SARA 

90 

0.07 

18.3 

64tn-34m(LO)-34in(T/R): 


BBO 

74 

0.29 

20.9 


CCR 

74 

0-22 

21.3 


SARA 

90 

0.10 

23.4 

64m-34m(LO)-34m(LO)-34m(T/R) : 

BBO 

74 

0.30 

25.1 


CCR 

76 

0.22 

25.9 


SARA 

90 

0.08 

28.3 


Table 2. Performance of various arraying schemes at Voyager 2 
Neptune encounter (Peak of day curve; BER = 0.005, 2 B,q = 
30 Hz) 


Array configuration 


Optimal 
modulation 
index, deg 

Radio 

loss, 

dB 

Maximum 
data rate, 
kbps 

64m-34m(T/R): 


BBO 

72 

0.51 

6.75 


CCR 

72 

0.45 

6.84 


SARA 

90 

0.25 

8.04 

64m-34m(LO)-34m(T/R): 


BBO 

72 

0.59 

8.53 


CCR 

72 

0.34 

9.04 


SARA 

90 

0.20 

10.2 

64m-34m(LO)-34m(LO)-34m(T/R): 

BBO 

70 

0.43 

10.4 


CCR 

74 

0.36 

11.1 


SARA 

90 

0.23 

12.2 













TO OTHER 
SLAVES 


Fig. 3. SARA with subcarrier demodulation 




Fig. 4. Degradation in loop SNR of SARA slave due to misalignment 
of subcaniers (a>sc = 360 kHz) 


Fig. 5. Degradation in loop SNR of SARA slave due to misalignment 
of data stream (Rj = 38.4 ksps; pf = .2, .4, .6, .8, 1.0) 
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Fig. 6. Data rate vs time: Voyager 2 Uranus encounter 
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MAXIMUM DATA RATE, kbps 


11 


ARRAY = 64m-34m(T/R) = 30 Hz BER = .005 

UPPER CURVE = SIDEBAND AIDED RECEIVER ARRAY (SARA) 0= 90® 
MIDDLE CURVE = COMBINED CARRIER REFERENCING (CCR) 6 = 72® 
LOWER CURVE = BASEBAND ONLY COMBINING (BBO) 6 = 72° 


AVAILABLE DATA RATES 



0 2 4 6 8 10 12 14 16 18 20 22 24 


TIME OF DAY - NEPTUNE ENCOUNTER 

Fig. 7. Data rate vs time: Voyager 2 Neptune encounter 
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MASTER 



TO BASEBAND 
COMBINER 


Fig. 8. A possible DSN implementation of SARA 















Appendix A 
Standard PLL Theory 


In a standard phase-locked loop, as shown in Fig. A-1, the 
input to the loop filter has the form 

z{t) - y4 sin (0 - ^ + n{t) 

where 0 is the phase of the input signal, ? is the loop’s 
estimate of that phase, and n{t) is noise. In our case, n{t) is 
white Gaussian noise with two-sided spectral density A^q/ 2. If 
the impulse response of the filter is/(r), then the output of the 
filter is 


The function 


I 


a{t) = I z(jc)/(r - x)dx +A^(0 


where N{t) is the filtered noise. It follows that 


^vco 


where Ky^Q is the frequency gain associated with the VCO. 
Upon taking Laplace transforms of the equation, 


H{s) = 




s-^KycoAF(s) 

is called the closed-loop transfer function of the loop. 

The phase jitter in the loop is given by 

(p’ = E 

where E is the expected value operator. If the phase 0 of the 
incoming signal is varying slowly compared to 0 (as is the case 
in SARA), then 

= E[E(^--$]'^ 

= E{H{s) <Kn) - m [<Kn) +N(s)lA]f 
= E[H(s) N(s)IA]'^ 




^io 
— loo 


ds 


s^s) = Ky^^AFQt) {sin [0(s)- $(s)] +N(s)/A} 




where F is the Laplace transform of the function/. For small 
phase errors, this may be linearized to 

s ${s) = Ky^^ A F(s) [0(j) - ^s) +{V(s)M] 


whence 


where 


2B^ = 


f /oo 

|//(s)P 

/oo 


ds 


is the closed-loop bandwidth of the system. The loop signal- 
to-noise ratio is just the reciprocal of the phase jitter. 


Kyco A F(s) 
VCO' 


X VCO — ^ ^ 


P = 
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Appendix B 

Sideband Aiding Without Hard Limiting 


In the body of the text, the slave receivers were analyzed 
assuming a hard-limited estimate of the data from the master 
was utilized to derive a carrier reference. Since some informa- 
tion is lost in the hard limiting of the output of the data 
integrator, the question arises as to whether performance 
could be enhanced if this output itself were used in deriving 
the slave carrier reference. In the following analysis, all signal 
references refer to Fig. 2. The incoming signal to the slave has 
the form 

d(t) cos (o3^t + +N^(t) 

where the various quantities have the same meanings as in 
Eq. (1). The carrier estimate from the VCO is 

r^(0 = - y/TK sin 

Neglecting the double frequency components, it follows that 

d(t) sin (4, - +N^it) 

where the noise is a white Gaussian process of one-sided 
spectral density K K* IV/Hz over the closed-loop band- 
width. If the date detector box of Fig. 1 does not perform the 
previously described hard-limiting function, then the signal 
that is sent to the slave for aiding has the form 

e{t) = c/(r-rp + « 3(0 

where d(t) can be plus or minus one and n^if) is a zero mean 
Gaussian random variable with variance N^j{2E^ over the 
symbol time 

With proper alignment, 

where N ^{i) = d{t) N ^{t) has approximately the same spectral 
density as over the closed-loop bandwidth. The other 


noise term, N 2 {t)n^{t), has zero mean and spectral density 
KK* ^V^q[A^q/(2E’j)] over the closed-loop bandwidth since 
N^it) and n^{t) are independent processes. From standard 
PLL theory (Appendix A), the loop SNR in the slave for this 
“soft-quantized” implementation is 



The corresponding loop SNR for the “hard-quantized” imple- 
mentation of SARA (Eq. 3) is 

The corresponding suppression factors are defined by 



and 

h, =erf2(N/FA^). 

These are shown plotted as a function of EJN^ in Fig. A-2. 
The figure shows that the breakpoint occurs at EJNq = 
-1.55 dB, with the SARA implementation performing better 
above this value and the soft-quantized better below this value. 
The range of values expected for a 64-m anterina arrayed with 
three 34-m antennas during Voyager Uranus encounter is 
-3 dB< EJNq < 0 dB. This brackets the breakpoint. If the 
SARA implementation is used, the slave loop SNR is degraded 
by at most 0.3 dB. Since slave loop SNRs at Voyager planetary 
encounters are expected to be above 14 dB with SARA, this 
degradation is virtually insignificant. Consequently, the rela- 
tive difficulties in passing a one-bit (hard-quantized) vs a 
many-bit (soft-quantized) data estimate to the slave receiver 
will determine which scheme should be implemented. 
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We examine the power-efficient capacity of the noiseless optical PPM channel It is 
shown that even though the capacity per photon can be made to increase without bound, 
the capacity per channel use (for best power efficiency) is always less than 2 nats per 
symbol Furthermore, it approaches 2 nats per symbol as the bandwidth expansion factor 
goes to infinity. 


I. Introduction 

In Ref. 1 a method was described for maximizing the 
energy efficiency of a noiseless PPM optical channel subject to 
simultaneous constraints on the throughput capacity and 
bandwidth. Specifically, we considered a Q-ary erasure channel 
often used to model a Q-ary PPM optical communication 
system. This channel is known to have a capacity per 
channel use of 

= (l - e ^ )ln2 nats/channel use (1) 

where is the intensity of the optical source as seen by the 
receiver during the PPM pulse (measured in photons/second) 
and AT is the PPM slot width in seconds. Since there are QAT 
seconds required for each channel use and each PPM symbol 
contains, on the average, N^AT photons, then one can specify 
the channel capacity in terms of 

= ^ = ^ nats/second (2) 


and 

= vvir " ^ nats/photon. (3) 

S 

In Ref, 1 we simultaneously fixed the throughput capacity 
Cj, and the slot width (or equivalently the system bandwidth) 
AT. Then, was plotted as a function of Q. From such 
plots one could determine Q*,the value of Q for which the 
capacity/photon, Cp^ , was maximized. In Ref. 1 we also 
observed that the product C^ATQ^ was approximately con- 
stant. In this article we will examine this latter topic in more 
detail. This will lead to what can effectively be interpreted as a 
Shannon limit for the noiseless, power-efficient optical PPM 
channel. 

II. Analysis 

We must first determine the requirements on the optimizing 
value of Q. A necessary condition for Q* can be obtained by 
setting the derivative of ^Ph with respect to Q equal to zero. 
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However, before this can be done note that if Cj> and LT are 
fixed, any variation of Q must be offset by compensating 
changes in (see Eq. 2). Thus, we must explicitly show the 
dependence of on Q, From (2) we have that 



Now, substituting (4) into (3) and differentiating we obtain 
the necessary condition 

where 

C^ATQ 

^®) = iiia- . (6) 

Thus, once we specify the product a = C>pAT (which is equiva- 
lently the capacity per slot), Q* can be determined by solving 
Eq. (5) numerically. 

Table 1 shows the results of such calculations for a wide 
range of a’s. Both 2* and C>pATQ* are shown for each a. 
Clearly, as a decreases there is a compensating increase in Q* 
such that for most values of a, Cj^ATQ* is slightly less than 2. 

We will now interpret these results in an even more interest- 
ing way. From Eqs.(l) and (2) we recognize that CrpATQ is 
simply C^, the channel capacity in nats/channel use. Thus, 
Cj^ATQ* is nothing more than evaluated at the most 
energy-efficient operating point. We shall denote this quantity 
CJ. Additionally, we note that if a. decreases, either AT 
decreases (bandwidth used increases) for fixed Cj, or the 
required throughput capacity decreases for fixed bandwidth. 
Both statements are equivalent to saying that the ratio of 
available bandwidth to information bandwidth is increasing. In 
Fig. 1 is plotted as a function of this bandwidth expansion. 
In the appendix we prove that <2 and approaches 2 in 
the limit. Thus, we have established the following fundamental 
property: 


For the noiseless optical PPM channel, the most 
energy-efficient use of the channel results in an 
information per channel use rate of less than 2 

nats per channel use. Furthermore, approaches 
2 nats per channel use as the bandwidth expansion 
factor approaches 


III. Discussion 

It is well known that the capacity of the noiseless PPM 
channel measured in nats/photon can be made infinite by 
allowing the word size (and hence the bandwidth expansion) 
to go to infinity. It is therefore surprising that the capacity in 
nats/channel use limits out at 2. This is because, as a decreases, 
Q* increases, Z* goes to zero and the average energy per pulse 
(i.e., per channel use) also goes to zero as seen by Eq. (4). 
Thus 


r* - ^ 

Pt N*AT 


goes to infinity because goes to zero. 

The formulation of bears a striking resemblance to the 
Shannon limit for the additive white gaussian noise channel. 
For the AWGN Shannon showed that the energy required to 
transmit reliably one bit of information, normalized by 
the one-sided noise power spectral density Nq is lower 
bounded by 



and that E^INq can be made to approach this bound as the 
bandwidth expansion of the signal approaches «>. For the 
noiseless optical PPM Channel, vanishes (quantum noise 
does not) so the energy required per bit can be made arbi- 
trarily small. However, the capacity per power efficient use of 
the channel is still limited and approaches 2 nats/channel use 
as the bandwidth expands. As in the gaussian channel case, the 
optical channel limiting behavior is obtained by a sequence of 
increasing complexity orthogonal (PPM in this case) modula- 
tion schemes. 
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Table 1. Optimized values of 0 and CjATQ 


a 

Q* 

Cj^TQ* 

10° 

0 

0 

O 

1 

p 

7.13 

1.42 

10“‘ 

15.7 

1.57 

10-2 

176 

1.76 

10-2 

1.83 X 10® 

1.83 

10-4 

1.87 X 104 

1.87 

10-2 

1.89 X 10* 

1.89 

10-4 

1.91 X 10* 

1.91 


1.92 X 102 

1.92 

10-® 

1.93 X 10° 

1.93 

10-9 

1.94 X 10° 

1.94 

10-10 

1.94 X 10*° 

1.94 

10-" 

1.95 X 10* 1 

1.95 

10-»2 

1.95 X 10*2 

1.95 

10-14 

1.96 X 10i4 

1.96 

10-20 

1.97 X 102° 

1.97 

10-30 

1.98 X 102° 

1.98 




Fig. 1. Variation of energy-efficient channel capacity /channel use 
with bandwidth expansion 
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Appendix 


We wish to show that < 2 and asymptotically Now, by expanding 1/1 - Z and In (1 - Z) in their correspond- 
approaches 2 as a->0. From (5) we have that the optimizing ing power series representations on the open interval of Z 
Q for any a satisfies the relation yields 

l-Z* " InQ*] ® 

where 


' ' n=0 n=l 


II 

n 


C^^TQ* q 
" Ing* " • 

Using this last equality we can rewrite the optimality condi- 
tion as 


Expanding and cancelling common terms gives 






+ • • • 


C* = f{Z*) 


where 


Each term on the left dominates the corresponding term on 
the right, so the assertion is clearly true whenever 


" Z + (1 - Z) In (1 - Z) 


From (2) we see that 


Z = 


-N„AT 
l-e ^ 


Z + (l -Z)ln(l -Z) > 0. 


By the same procedure, this is true if 


so that 0<Z<1. (Z=l corresponds to infinite energy per ' ^^P^^ding both sides gives 
pulse whereas Z = 0 corresponds to no energy per pulse.) We 

now claim that/(Z) < 2 for all Z e (0, 1). To prove this we ask Z + + • • ■> Z +— +— +— 


Z^ 

Z + (l-Z) In (l-Z) 


9 

<2 


(A-1) 


We will prove in a moment that Z + (1 - Z) In (1 - Z) > 0 for 
Z e (0, 1). Thus, (A-1) can be rewritten as 


Ml 


l-Z 


<ln(l-Z) 


which is clearly true. 

Finally, we note that as the bandwidth expansion factor 
increases, a decreases and Z 0. This is true whether or not 
Q* increases (which it does) by virtue of the fact that 
Cj^ATQ* is bounded. For this reason we are interested in /(O). 
By applying L’Hospital’s rule twice to/(Z) we see that/(0) 
does indeed equal 2. Thus, CJ is less than 2 and approaches 2 
in the limit as a -> 0. 
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In this article a method of determining the group delay effect of DSN microwave 
components on multimegabit telemetry detection is described. That is, when the gain and 
phase shift characteristics of the overall system between the modulator and demodulator 
are given, the group delay effect for the intended data rate can be assessed by examining 
the equation for the recovered signal (equation *) for several rectangular pulse inputs 
with different duty cycles. For example, the group delay loss of 30 Mbitfsec telemetry 
with only a Block II X-band maser is less than 0.2 dB even with a lO-MHz center 
frequency offset. 


I. Introduction 

A simplified model for the transmission and reception of multimegabit telemetry is shown in Fig. 1. The modulating signal 
x(0 is a unit amplitude random binary data signal with bit duration Tj,. The phase-modulated signal >>(r) is 


>>(r) = K sin [iTtf^t + (px(t)] 


where 4(=l/r^) is the carrier frequency (T^ « T^), and <p is the modulation index. The /sT in^O is the signal amplitude, and will 
be taken as unity for convenience. If the channel is ideal, thenr(/) =^^(0> 


y(t) cos (2nf^t - ) = sin (<px(t) + |3j ) + sin (47r/^f + 4>x(t) - ) 


and 


x(r) = sin(0x(O + /3j) 
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where |3j is the phase difference between the two reference oscillators at the modulator and demodulator and will also be ignored 
for later calculations. The modulation index (j> is usually chosen to be less than tt/ 2 for both maximizing energy of recovered signal 
x{t) and efficient carrier referencing. Since we are only interested in the modulation component, we assume 0 = ;r/2. 

The overall system between the modulator and demodulator (any combinations of amplifiers, mixers, diplexers, filters, 
antennas, waveguides, transmission lines, etc.) can be considered as an effective channel bandpass filter. Typical gain and phase 
shift characteristics of the resulting channel are shown in Fig. 2. Let /q be the center of frequency of B?F ;g(f) and|3(/) are the 
relative gain and phase shift at frequency / to those at /q. Usually they have symmetry such that 

g(f^ - A/) = g(f^^Af)mdp(f^ - A/) = . 

It is well known that when a signal containing components throughout the frequency band of the BPF is divided into a 
set of narrowband signals, each of these will be subject to a different time delay. Hence the output signal must have a different 
waveform from the input, and phase distortion appears along with the amplitude distortion. 

It is very difficult to analyze the effect of group delay for a random data signal. Hence for analysis we suggest the use of 
periodic rectangular pulse modulating signals with period Ts{=llfs = nT),) and duty cycle l/n, for several different «’s, instead of 
the actual random data signal. Since in most cases the bit time-bandwidth product BT^ is not smaller than one, the group delay 
together with amplitude distortion of one bit affects only itself and its adjacent neighbors. 

In the next section, we describe the way of finding the recovered signal x(r) including the case where the center frequency of 
the effective bandpass filter is not matched to the carrier frequency of the incoming signal. Then we show that the group delay 
effect can be determined by the losses of the integrated values of5c(r) with group delay compared to those ofx(r) without group 
delay. 


II. The Recovered Signal x(t) 

The Fourier series representation of a periodic signal z(t) with period r(=l//) is 

z(t) = Z(kf) z(kf) =jrf dt . 

fc=-~ (W * 

Let /j, = mnf^ = mnT^) where m represents the number of carrier cycles per bit time and n represents the number of bit times 

over which the bit stream is periodic. The m and n are assumed to be integers so that_y(r) is also a periodic signal with period 
(see Fig. 3) and the Fourier series representation can be used. Then with = mnf^ and its /rth Fourier coefficient 
are 

>>(r) = sin |27r/«n;^t + y 

where 


= 


' i,o<k- crsi<7’^/2« 

£ = 0,±l,+2, 

.-l,TJ2n<\t-HT^\<Tj2 


and 


mn^ "'s' 


sin [(m/I - k)Tr/n] ^ sin [(mn +k)Tr/n] _ sin [(mn ~ k)7i] _ sin [(mn + k)n] 
{mn - /:)7T {mn + k)it {mn - k)2n {mn + k)lTi 


60 


Note that 


Also, when mn » £ > 0 




r„J(m« * «)/,] « y„„ [(m. - 8)/.l = 


Ctt 


Therefore 


( fTj' 

^ (ptnlnf^t) + ^ sin ^ 2;r//] + cos [(mn - £) 2nf,t] } 

' c=i 


where m is some number less than mn/2. Also 


r(t) -^cos(mn2irf/) 


sin (9.7rln) 


c=i 


• ^g[(mn + • cos {(mn + fi) 2nf^t - P [(mn + fi)/^] } 


^g[(mn - £)/^] • cos {(mn - Si)2nf^t - P [(mn - 2)fJ }) 


where w" is determined by the bandwidth of BPF which is usually much less than m\ 


Then finally 


where 


x(t) = const. + ^ sin (2n/n) g)y^] cos {£ 27 r// - ^ [(m« + £)/J } 


fi=i 


+ g [(mn - cos {2271// + p [(mn - }) 

= £ 1MM») cos 


C=1 


• cos (P [(mn + £)41 + p [(mn - 2)// 


(*) 


61 



and 


^ ' /g[(mn + g)4] sin (j3[(mw + fi)/J } - g[(m« - g)/J sin {|3 [(m/i - g)4] } \ 

\ g[{mn + g)/^] cos {P [(m/I + g)/J } + g[(mn - g)/J cos {p [(mrt - g)/J } / ’ 

The additive constant is just the dc bias term. If the BPF has symmetry, i.e.,^Q - Af) = g(f +A/) andj3(/*Q - A/) = -|3 (/q + A/) 

and if /<> = 4 , then = 2^((m« + g) 4 ) and 0 ^ = p [(mn + g)/J . 


III. Example and Conclusion 

To see the group delay effect on 30 Mbit/sec telemetry with a Block II X-band maser, the shapes of x(/) are shown in Figs, 4 
and 5 for n = 2 and 10 respectively, using the data from Fig. 6 of Ref. 1 and equation (*) mihf^ ~ fo 0^ ~ 281), The same 
signal conditions are repeated in Figs. 6 and 7 with =/q ~ 15 MHz. The dotted lines in Fig. 4 to Fig. 7 are the shapes of 3c(0’s 
without group delay; i.e., only amplitude distortions are considered. When we use an integrator for detection of x(t) with 
integrating duration TJrty the losses of the integrated values of x(t) with group delay compared to those of x(t) without group 
delay are shown in Fig. 8 for various center frequency offsets. From these observations, we can conclude that the group delay 
effect on 30 Mbit/sec telemetry with Block II X-band maser is not significantly deleterious even with a substantial center 
frequency offset. 
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This article describes a new technique for displaying the information needed by DSN 
operators for monitoring station performance. Known as a ^^Qualitative Monitor, ” it pre- 
sents continuous variables in quasi-analog form on a digital cathode ray tube (CRT) 
display. Color changes, field reverses and blinking symbols assist the operator in identify- 
ing variables that are within acceptable limits and in performing corrective action when 
needed. 


I. Motivation 

A significant part of a Deep Space Station (DSS) or Net- 
work Operations and Control Center (NOCC) operator’s task is 
to determine if the station equipment is performing as it 
should, and to take corrective action if it is not. In many cases 
the assessment of performance requires the operator to study 
numbers on a CRT display and try to determine if they are 
within specified limits and if they are rising, falling or station- 
ary. If he decides that a change is needed and starts to move a 
control, there is the additional question as to which way to 
move the control to get the desired effect. When the displayed 
number is a multidigit quantity, and perhaps displayed in 
“floating-point” format because this made the programmer’s 
task easier, it is sometimes hard to decide just what the vari- 
able is doing, especially if the quantity is jittering around a 
threshold for change of scaling. 

The Qualitative Monitor Display presents a key set of vari- 
ables to the operator in a uniform format which permits him 
to tell at a glance the state of the information presented and 
whether any further study is need. Additionally, the display 
aids the operator in taking any necessary corrective action. 


II. History 

Although the DSN has made significant progress in the pro- 
cessing of data and in the available throughput of information, 
it has made significant but unsatisfactory choices in the man- 
ner in which information is presented to people operating or 
monitoring equipment and system performance. In this section, 
an analog is drawn by examining data presentation in aviation. 
It seems clear that DSN programmers took technically sound 
display approaches which did not yield a good computer- 
human interface. 

Through the early history of aviation until about the mid- 
dle 1950’s, instruments were all analog devices: fuel gauges, oil 
temperature indicators, artificial horizon gyroscopes, etc. Safe 
operating limits were usually shown by colored striping on the 
meter face, usually green for “good,” and yellow and red for 
“marginal” and “bad.” The pilot could tell with a quick glance 
at the instruments if all the needles were “in the green” leaving 
him free to concentrate on maneuvering and navigation. 

As the science of digital instrumentation developed, engi- 
neers were attracted by digital numeric displays that could 
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show numbers that represented the variables and could do this 
more “accurately,” i.e., with more resolution. However, it 
developed that this was not the sort of information the pilot 
needed for control. When landing or taking off, he wanted to 
know if the engine was performing within required margins, 
not the exact temperature at a cyhnder head or the precise oil 
temperature. For the actual control of the plane, a digital 
presentation of the heading added another complication: if 
your heading is now 303 degrees, and you want to turn to a 
heading of 129 degrees, should you turn right or left? In par- 
ticular, the pilot must be able to immediately tell the attitude 
of the aircraft in instrument (within clouds) conditions. He 
could not afford the luxury of a complex interpretation of 
complex but accurate digital information. Hence the aircraft 
is instrumented with analog artificial horizons and a circular 
heading compass. 

The numeric displays in aircraft are usually accompanied 
by bar-graph displays, colored regions on the face of an instru- 
ment or “good/bad” signal lights. Although CRT’s are now in 
use on the more sophisticated aircraft, flight control informa- 
tion is presented as the old-fashioned artificial horizon, flight 
director, or digital gyroscope on the face of the screen. The 
analog display now is being built of electrons instead of cams, 
levers, gears, etc. The reason is quite simple: humans can 
interpret analog displays of complex information much more 
easily than digital displays. 

ill. Approach 

The Qualitative Monitor Display is made up of a series of 
horizontal lines. A portion of each line represents the scale of 
an analog meter, with a rectangular marker moving back and 
forth like the pointer on an edge-view meter. The “meter” 
excursions are limited to 21-character positions, corresponding 
to scale values of -10 to +10 at the edges and 0 at the center. 
At the left of each line is a label which identifies the variable, 
for example, NOISE or S-DOPPLER. At the right of the scale, 
the numeric value is shown with limited scaling, and at the 
extreme right there is a three-character label that identifies a 
more detailed display where the operator can find further 
information on the variable presented in that line. An example 
of the qualitative display is shown in Fig. 1. 

On any particular line, as the marker moves along its scale, 
it changes color; between +/- 5, it is green for “good” values, 
between 5 and 7 and between -5 and -7 it is blue for “mar- 
ginal,” and at the outer extremes it is red for “danger” or 
“bad.” The variable being plotted is hard-limited so that off- 
scale values all show in the most extreme position overlaid 
with a > or < symbol that blinks to attract the eye. If a vari- 
able is so far out of limits that it cannot be reasonably scaled, 
then it is presented in its log base 10 form with the cursor 


overlaid with a blinking “?”. Thus, the variable will always be 
visible on the line and the variable’s state is always apparent to 
the operator. 

Note that the display format of one quality bar above the 
other also facilitates rapid interpretation of the state of the 
entire set of data on the display. If everything is performing 
within specification, the operator’s eye sees the pattern of a 
vertical column of green cursors. If any single cursor or set of 
cursors is “out-of-line” it is readily visible. Hence the operator 
can quickly assess the condition of many diverse data types 
without exact interpretation of any of them. 

IV. An Example 

In order to assess the usefulness of the display, an example 
was programmed into the software supporting the Metric Data 
Display (MDA), the DSN subsystem which acquires and pro- 
cesses tracking data. Although the display is general and easily 
implemented on any computer system, the specific example to 
be examined was programmed on the Modcomp II used as the 
processor for the MDA. 

The variables to be displayed usually have predicted values 
and tolerances associated with them. To create the display, the 
variables are first biased so the predicted value will appear at 
the 0 position on the line and multiplied by a scale factor 
derived from the tolerance so that the values that are within 
the tolerance will appear in the middle 50% of the line. The 
scaled value determines the cursor position. If the value is in 
the range -10 to +10, the cursor is placed in one of the 21 
available character positions on the line. If the value is greater 
than 10 or less than -10, the cursor is placed at the extreme of 
the line with a flashing “>” or “<” symbol as appropriate. 

If the absolute value exceeds a certain value, recall that its 
logarithmic representation is presented. For the Modcomp, the 
maximum integer value is 32767; hence that is the selected 
limit. If the variable exceeds that value, the line cursor is over- 
laid with a flashing question mark and the integer part of the 
exponent is used to position the cursor. 

A general-purpose Modcomp II qualitative display sub- 
routine was written and has been tested in the MDA. The 
subroutine, named QMD, first resets the current DSN display 
terminal, a Megadata, to put it into a known state. Next, the 
characters which create the lines are formatted and sent to the 
terminal. Any line that represents a variable is sent in four 
parts: the label field with appropriate scaling information; the 
marker in its appropriate position and with the flashing “?”, 
“>”, or “<” overlays; the numerical value of the variable; and 
the key to a more detailed display. On successive displays, the 
new marker position is computed and compared against the 


67 


old position, and if the marker need not be moved, only the 
numerical value is updated. 

The Fortran program for the Qualitative Monitor Display 
was developed in two main phases. The first was to establish 
the exact control sequences needed to command the various 
Megadata functions, and then to test several alternative 
schemes for making the line and marker. Various strategies for 
handling labels and scaling were tried. This phase was per- 
formed in the Section 338 laboratory using any of the several 
Modcomp II machines. 

The next phase was to incorporate the features of the QMD 
in an overlay to be called by the MDA operational program. 
This necessitated adding QMD to the list of allowable direc- 
tives, modifying the executive routine so that it would call the 
appropriate overlay and cataloging the QMD overlay into the 
load module. 

The specific display used to test the QMD concept, shown 
in Fig. 2, combined variables from two frequently used MDA 
displays; doppler residuals and pseudo-DRVID were extracted 
from the appropriate data base. All values were scaled appro- 
priate to their expected values. During tuning operations, how- 
ever, the doppler residual could be very large. The logarithmic 
display feature, however, allows the observer to make use of 


the cursor position, even though it is far off scale, to adjust the 
received frequency appropriately to “zero” the residual. 

The developmental version of the MDA program has been 
extensively tested in the DSN Data Systems Section lab, on 
the MDA at CTA 21, and it is expected to be taken to Gold- 
stone for real-life testing. 


V. Results and Conclusion 

The display has been demonstrated to a broad spectrum of 
potential DSN users. The comments received have been very 
encouraging and the authors expect to see the display used for 
various subsystems in the DSN’s Mark IVA (Networks Consol- 
idation) upgrade. DSN consultants in the area of the human 
interface have reviewed the display and were very pleased with 
the concept in terms of rapid and concise interpretation of 
information. 

The qualitative monitor technique can be a useful one for 
the DSN, particularly for monitor and control. The monitor 
and control operator must observe many data values from 
many different pieces of equipment and usually in different 
formats and ranges of values. The qualitative monitor offers a 
concise and consistent means to observe station operation. 
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Fig. 2. Qualitative Monitor Display for MDA 
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This article describes the test results obtained for the Master Timing System (MTS) 
manufactured for Goddard Space Flight Center and deployed at the Tracking and Data 
Relay Satellite Systems Station in White Sands, N,M. 


I. Introduction 

The tests described in this article were performed on the 
Master Timing System (MTS) presently in use at the Tracking 
and Data Relay Satellite Systems (TDRSS) Ground Station at 
White Sands, NM. This system was manufactured for Goddard 
Space Flight Center (GSFC) by TRAK Systems, and is similar 
in all respects to the master clocks to be utilized in the MK 
IVA Frequency and Timing Subsystem. Three of these clocks 
are being purchased by GSFC to be transferred to JPL upon 
consolidation of the DSN and STDN. The fourth clock is being 
procured by JPL for future implementation into the DSN. 
Operating parameters and performance criteria were necessary 
to ensure subsystem performance and institute design changes 
for the new clocks. 

This report is a composite of two test runs performed on 
the MTS. The primary tests were conducted on 9 December 
1980 through 11 December 1980. These tests were inconclu- 
sive due to test equipment malfunctions. Follow-up tests were 
conducted on 3 February 1981 through 5 February 1981, 
with all results being verified. 


II. Test Descriptions 

During the two test runs, the following tests were per- 
formed and the results of these tests explained. 


(1) Port-to-port isolation. The output ports were tested for 
isolation by observing the output while subjecting the 
adjacent ports to varying loads, i.e., open load to short 
circuit. Both amphtude and phase of the test ports 
were monitored. The first trial of this test produced a 
complete failure of all outputs; however, subsequent 
exhaustive tests failed to produce any perturbations. 
No significant change in amplitude or phase could be 
measured at any output port. 

(2) Sine wave outputs. All sine wave outputs were mea- 
sured for isolation, amplitude and distortion. Distor- 
tion was not measured by a qualitative means; how- 
ever, no distortion could be detected on a 400-MHz 
storage oscilloscope. All outputs were within specified 
limits for isolation (<0.1% change in amplitude with 
any output open or shorted) and amplitude (adjust- 
able from 2 to 5 V(p-p)), and the distortion could be 
stated to be below five percent. 

(3) Serial time codes. The serial time codes were clean and 
within specifications. The amplitude was measured to 
be 3 V(p-p), and the mark-space ratio was adjustable 
from 2:1 to 4:1. 

(4) Timing pulses. Pulse amplitude, wave shape, rise time, 
and jitter were measured, and the results of each cate- 
gory are presented below. All measurements were made 
with the output terminated into 50 ohms. 
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(a) Amplitude. The positive or “on” portion of each 
pulse was measured to be +3.5 ±0.2 V and the neg- 
ative or “off’ portion was +0.2 + 0.1 V. 

(b) Wave shape. Each set of pulses displayed a differ- 
ent shape, which can best be described by draw- 
ings of those waveforms. These drawings are pre- 
sented in Fig. 1. 

(c) Rise time. Rise time was measured utilizing a 
400-MHz storage oscilloscope. All pulses displayed 
a rise time of less than 15 nsec between the 10 and 
90% points. 

(d) Jitter. Pulse jitter was measured by redundant sys- 
tems utilizing an HP 5370 time interval counter 
and a standard clock with an HP 9835 computer 
doing the controlling and data reduction for each 
of the two systems. The in-house frequency stan- 
dard was utilized for external synchronization of 
the counters and clocks to insure maximum coher- 
ence of each test channel. The standard clocks 
were measured to have pulse jitter of less than 50 
picosec. The pulse output of the MTS was com- 
pared against these standard clocks for this test. 
The data were taken at 100-sample increments and 
an accumulation of the drift, rms (la deviation), 
positive and negative peak deviation was compiled 
by the computers. Plots of the reduced data are 
presented for each of the pulse rates in Figs. 2 
through 5. The maximum accumulated values for 
each pulse rate are given in fractional seconds 
below: 


Pulse 

rate, 

pps 

RMS jitter 

Positive 
peak jitter 

Negative 
peak jitter 

1 

4.73 X 10-10 

1.25 X 10-9 

-1.44X 10-9 

10 

4.75 X 10-9 

1.47 X 10-9 

-9.34 X 10-9 

100 

8.03 X 10-10 

2.24 X 10-9 

-3.05 X 10-9 

1000 

4.41 X 10-9 

1.08 X 10-8 

-9.24 X 10-9 


(e) Duty cycle. The pulses displayed at 80% on-time 
and a 25% off-time duty cycle as specified by God- 
dard. This can be modified to provide a 10-/isec 
pulse width to meet JPL requirements. 

(f) Slew rates. All pulse slew rates were surprisingly 
accurate; even the 50-nsec/sec rate was exactly 50 
nsec per step. 

(g) Pulse coincidence. The timing coincidence be- 
tween different output ports was less than 5 nsec 
between the slowest and fastest pulse of a set (e.g., 
all 1 pps ports). 

(h) Zero crossing coincidence. The worst case offset of 
the zero crossings of the 5 -MHz reference signal 
and the leading or rising edge of the 1 pps was less 
than 50 nsec. 


III. Conclusion 

The test results indicate that the system performs within 
the specifications established by Goddard Space Flight Center. 
The overall system appears to be well designed and built to a 
good commercial grade. The Triple Redundant Time Code 
Generators, whose outputs are majority voted, performed ade- 
quately as far as could be determined on an operational sys- 
tem. The automatic fault isolation circuitry could not be 
tested, due to the system being in an operational situation. 

This timing system has been in operation at the TDRSS sta- 
tion for more than one year with no reported failures. 

The failure noted in paragraph II (1) apparently was caused 
by a disturbed reference cable or loose connection external to 
the clock. 

The inconclusive results obtained during the first test run 
were caused by a perturbation in the data, for which the cause 
could not be determined. The first test run did not utilize 
redundant test systems. 
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Fig. 1. Timing pulse waveforms 
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Fig. 2. 1-pps timing pulse 
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Fig. 3. 10-pps timing pulse 


73 



PULSE STABILITY vs TIME 


POSITIVE PEAK 
DRIFT 

NEGATIVE PEAK 
RMS JIHER 


POSITIVE PEAK 
DRIFT 

NEGATIVE PEAK 
RMS JITTER 


SECONDS 




Fig. 4. 100-pps timing pulse 
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This is the second article on the subject of Optimum Equipment Maintenance and j or 
Replacement Policy which employs the optimization technique called Markov Decision 
Process. In the first article, dynamic programming was utilized as an alternative 
optimization technique to determine an optimal policy over a given time period. 
According to a joint effect of the probabilistic transition of states and the sequence of 
decision making, the optimal policy is sought such that a set of decisions optimizes the 
long-run expected average cost (or profit) per unit time. Provision of an alternative 
measure for the expected long-run total discounted cost is also considered. A computer 
program based on the concept of the Markov Decision Process was developed and tested. 
The program code listing, the statement of a sample problem, and the computed results 
are presented in this report. 


I. Introduction 

In the first article {TDA Progress Report 42-66: September 
and October 1981), an optimal decision-making policy utiliz- 
ing the dynamic programming technique was presented. The 
aim was to make optimal decisions, over a finite number of 
time periods, regarding the equipment maintenance and/or 
replacement for a given system. 

When a system is required to be existing indefinitely, a best 
operation policy which gives the optimal long-run cost (or 
profit) may be estimated by successive approximations with 
dynamic programming techniques, providing a very large 
number of time periods is assumed. However, there is no way 
of knowing when to terminate the successive approximations; 
there is no procedure for deciding how large a number of time 
periods is sufficient. 


This report presents another decision-making technique to 
obtain a long-run optimal policy utilizing the Markov Decision 
Process Concept. The optimal policy is evolved over time 
periods according to the joint effect of the probabilistic 
transition of the condition of the system and the sequence of 
decision making. It is assumed that for a system under 
consideration, there exists a policy at any time period. The 
system changes to a new state at the next time period accord- 
ing to a known probability after a decision is made at the 
present time period and at the present state. 

In such decision making, different transition matrixes 
result, corresponding to the decisions made at each observed 
time period. Accordingly, a change of state as well as the 
associated value of cost (or profit) of the system in the next 
time period is governed by the transition matrix. The optimal 
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is sought such that the set of decisions made will optimize the 
long-run expected average cost (or profit) per unit time period. 

However, in systems involved with a long time-horizon, the 
changing time value of money is of importance, and the* 
expected long-run total discounted cost (or profit) should be 
determined with respect to a specific discount factor. The 
discount factor a is the present value of one unit money in one 
time period, expressed as: 


where i is the rate of return on the money for one time period. 
After M periods, a unit of money will be worth 

Derivation of mathematical equations, discussion of the 
policy-improvement computational algorithm, and a sample 
problem are presented in the following sections. A computer 
program employing the discussed algorithm is given in Appen- 
dix A. 


il. Theoretical Model 

Brief descriptions of the Markov Decision Process and the 
necessary equations are discussed as follows. More detailed 
derivations of the equations may be found in Refs. 1, 2. 

Consider a system which at a particular time period (^ = 1, 
2, 3, . . .) is in one state i out of A/ states. The system changes 
from one of these admissible states, /, to another state, /, ruled 
by a transition matrix, P = (p^y). The elements, p,y, are defined 
as the probabilities the system is in state / at f, given that it 
was in state i at (t -1). Further, it is assumed that the 
transition matrix/^ is not time-dependent. 

Let be the expected cost (or profit) incurred when 
the system, which originally is in state /, changed after a 
decision k is made to a state / at the next observed time 
period. Then, 

M 

/=! 

where is the cost incurred at the first observed time period 
as a result of the current state i and the decision D^{R) = k 
when operating under policy/?. 

By introducing F/^(/?) as the total expected cost of a 
system starting in state i (at the first observed time period) and 


evolving for N time periods following a policy/?, the recursive 
equation can be written as 

M 

Vf {R) = +53 p,! (k) (R) (3) 

/=! 

The second term on the right-hand side of Eq. (3) is the total 
expected cost of the system evolving over the remaining 
(A^- 1) time periods. 

Let g(R) be the long-run expected average cost per unit 
time following a policy/?. As one of the Markovian properties, 
it can be shown that the g(R) is independent of the starting 
state i as the number of time periods N approaches infinity. 
Hence, Vf^(R) may be approximated by 

V^(R) ^ Ng(R) + V^(R) (4) 

where V^(R) can be interpreted as the effect on the total 
expected cost due to starting in state /. Thus, from Eq. (4) 

Vj^(R) ^ V^(R)~ V.(R) (5) 

the term [ V^{R) - ^y(/?)] is a measure of the effect of starting 
in state / rather than state /. Substituting the linear, approxi- 
mate relations of Eqs. (4, 5) into Eq. (3) leads to the recursive 
equation 

M 

g(R) + (R) = + X) P/y W (^) (6) 

/=1 

Equation (6) represents one of the M equations corresponding 
to the state / for / = 1, 2, . . . , M. 

When a system operates according to the Markov chain, 
there are needed (M + 1) values of^(/?), Kj(/?), V^iR ), . . . , 
Vj^(R) which satisfy the set of M equations of the form 
Eq. (6). Note that there are M equations and (M + 1) un- 
knowns; one of the unknowns, say F^(/?), can be arbitrarily 
set to equal zero. Following a given policy/?, the correspond- 
ing values of g(RX V^(R), V^(R ), . . . , k"(Af_i)(/^) can then be 
obtained by solving the set of Af simultaneous linear equations. 

In principle, all policies can be enumerated to find the 
policy which optimizes the g(R)- However, even for a 
moderate number of states and decisions, this enumeration 
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technique is cumbersome. A different approach called Policy- 
Improvement can be used to evaluate policies and find the 
optimal set of decisions without a complete enumeration. The 
mechanism of this algorithm is presented in the next section. 


III. Computational Algorithm 

The Policy-Improvement algorithm consists of two steps: 
the Value-Determination step and the Policy-Improvement 
step. These steps are described as follows. 

(1) Value-Determination Step: For an arbitrary policy 
(with decisions D^(R^)= and the corresponding 
values of Pijik^), and P^(^i) = 0), this step 

solves the set of M equations of Eq. (6), or 

M 

«(«,)+ (',(«,) = C„ r,(R,) (7) 

/=! 


for /= 1, 2, . . . ,M. Hence, values of the 

V^(Rj) are obtained under 

policy R 1 . 

(2) Policy-Improvement Step: Using the above calculated 
values of the F’s, find the alternative policy 
that for each state /, ~ ^2 decision which 

optimizes^(i?2)> 


M 


emj = C, 


^ Ep„(y;(«.)- n<«i) w 


ik /= 1 




Af 


/=! 


for /= 1, 2, . . . ,M. Here, the is the expected long-run 
total discounted cost of the system starting in state i and 
continuing indefinitely. The VfR) can be evaluated in a 
similar fashion as computing the average cost. 

A computer program is written to incorporate both the 
average cost of Eq. (6) and the discount cost of Eq. (10). The 
optimal policy is determined utilizing the Policy-Improvement 
algorithm. This BASIC program code is presented in Appen- 
dix A. A sample problem, adapted from Ref. 2 and presented 
in the next section, is used for testing purposes. The calculated 
results of both the averaged cost and the discounted cost are 
presented in Figs. 1 and 2, respectively. 


IV. Sample Problem 

For the purpose of testing the computet program, a sample 
problem is taken from Ref. 2 and summarized as follows. 

The condition of a given system is inspected and classified 
into one of four possible states as shown in Table 1. It is also 
assumed that the state of the system evolves according to some 
known probabilistic transition matrix given in Table 2. After 
each periodic inspection of the system, a decision must be 
made as to which action to take: Decision 1 is doing nothing; 
Decision 2 is overhauling the system; Decision 3 is replacing 
the system. 

In addition, the following assumptions are made: 


That is, for each state i=l, 2, find the 

appropriate value of k 2 such that 

OPTIMUM [g(R^)] (9) 

= 1,2,... ,JC 

In turn, let optimal value of k 2 , 

which defines a new policy ^ 2 - 

Using the new policy R 2 , the Value-Determination step is 
repeated. This iterative procedure continues until two succes- 
sive iterations lead to identical policies, which signifies that the 
optimal policy has been obtained. 

If the expected long-run total discounted cost is of interest, 
the above algorithm can be used with a modification to the 
recursive equation Eq. (6). With a specified discount factor a, 
the recursive relation of Eq. (6) can be modified as 


(1) When the system becomes inoperable (State 4) and 
replaced (Decision 3), the system is found to be in 
State 1 at the time of regular inspection. It is assumed 
that the total cost incurred when the system is in 
State 4 is the sum of a replacement cost of $4000 and a 
cost of lost production of $2000. 

(2) When the system is overhauled, the system is returned 
to State 2 (operable with minor deterioration) at the 
time of regular inspection at the end of next time 
period. The cost of the overhaul process is taken as 
$2000 and requires one time period to complete. 

(3) When the system is in States 2 or 3, defective items 
may be produced during the following operating 
period. The expected costs due to producing defective 
items are $1000 when the system is operable with 
minor deterioration and $3000 when the system is 
■operable with major deterioration. 
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(4) The total expected cost incurred per one time period 
depends on the state the system is in and the decision 
made. The total expected costs (the maintenance cost, 
the cost due to producing defective items, and the cost 
from lost production) are tabulated in Table 3. 

The above information completed the necessary inputs to 
the computer program. Figure 1 presents the optimal policies 
and the expected average cost of the sample problem. As the 
result of the Markov Decision Process, an average cost of 
$1667 can be expected when the policy is to do nothing when 
the system is found to be in States 1 and 2, to overhauling the 
system when it is in State 3, and to replace the system when it 
is in State 4. 

In the second case, as presented in Fig. 2, an interest rate of 
11% (or discount factor of 0.9) was assigned. However, with 
the same policy as in Case 1, a discounted cost of $14,950 can 


be expected if the system started in State 1, $16,260 if it 
started in State 2, etc. 


V. Summary 

The Policy-Improvement algorithm using a Markov Decision 
Process is incorporated in a computer program and tested with 
a sample problem on a Hewlett-Packard 2647A terminal. This 
computer program is capable of finding the best maintenance 
policy with respect to an optimal long-run average cost or the 
long-run discounted cost for a system with known transition 
probabilities. 

From the standpoint of management and operation, the 
algorithm provides a useful tool in obtaining an optimal 
maintenance schedule which gives the best return on capital 
invested. 
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Table 1 . States of the system 


State 

Condition 

1 

Good as new 

2 

Operable with minor deterioration 

3 

Operable with major deterioration 

4 

Inoperable 


Table 2. Transition matrix of the system 


State 

1 

2 

3 

4 

1 

0 

7/8 

1/16 

1/16 

2 

0 

3/4 

1/8 

1/8 

3 

0 

0 

1/2 

1/2 

4 

0 

0 

0 

1 


Table 3. Total expected cost per one time period 


State 


Decision 

* 

1 

2 

3 

1 

0 

$4000 

$6000 

2 

$1000 

$4000 

$6000 

3 

$3000 

$4000 

$6000 

4 

- 

- 

$6000 


* * JK # )K INPUT INFORMATION #)!(### 


MSTATE 4 

NDECISION 3 
MAXTRIAL iO 
MAXIMUM 0 

DISCOUNT 0 


STATE 

I 


DECISION 

D(I) 


TIE-BREAKER 

TIED(I) 


1 i 0 

2 10 

3 10 

4 3 0 


# * # FINAL RESULTS 


LONG-TERM AVERAGE COST/RETURN= 




1666.67 


STATE POLICY VALUE 


1 

2 

3 

4 


1 1666.67 

1 1666.67 

2 1666.67 

3 1666.67 





JK * * # )f( INPUT INFORMATION )«c * 3«c * * 


MSTATE 4 

NDECISION 3 
MAXTRIAL iO 
MAXIMUM 0 

DISCOUNT ii.iiU 


STATE 

I 


DECISION TIE-BREAKER 
D(I) TIED(I) 


1 

2 

3 

4 


i 

i 

i 

3 


0 

0 

0 

0 


* * * * # * )K FINAL RESULTS * * jK * * ♦ * 


STATE POLICY UALUE 


1 

2 

3 

4 


1 14948.6 

1 16261.6 

2 18635.5 

3 19453.7 




Appendix A 

Computer Program Listing 


i REM 
a REM 


3 REM EQUIPMENT MAINTENANCE POLICY 

A REM MARKOVIAN DECISION ALGORITHM 

5 REM 

6 REM 


10 DIM Irow<iO),Jcol(iO),Jord<iO),Y(iO) 

11 DIM Dd(10>,D(iO),P(S,5,5),R(5,5,5),SuM(5,5),TiGd(iO),Q(5,5) 

15 LONG A<U,ii),X(iO),Eps,SiMul 

16 INTEGER D, Dd , Tr ial ,Tied 


REM 

REM 

REM 

REM 

REM 

REM 


20 REM 




100 

REM. 


BEGIN OF DATA LIST 

.REM 

102 

DATA 

0;6 

1 ASSIGN I/O DEVICES 


105 

DATA 

0 

! PRINTOUT OPTION 


110 

DATA 

4 

1 NUMBER OF STATES IN CONSIDERATION 


120 

DATA 

3 

! NUMBER OF DECISION ALTERNATIVES 


130 

DATA 

0 

! MAXIMIZE COST/RETURN IF >=1 


140 

DATA 

10 

! MAXIMUM TRIALS ALLOWED 


150 

DATA 

11 .1111 

! DISCOUNT RATE IN WHOLE NUMBER 


160 

REM. 


PROBABILITY MATRIX, P(I,J,K) 

.REM 

170 

DATA 

0, .875, . 

0625, . 0625 


180 

DATA 

0,1, 0,0 



190 

DATA 

1,0, 0,0 



200 

DATA 

0,.75,.l 

25, . 125 


210 

DATA 

0,1, 0,0 



220 

DATA 

1,0, 0,0 



230 

DATA 

0,0, .5, . 

5 


240 

DATA 

0,1, 0,0 



250 

DATA 

1,0, 0,0 



260 

DATA 

0,0, 0,1 



262 

DATA 

0,1, 0,0 



264 

DATA 

1,0, 0,0 



270 

REM. 


COST/RETURN MATRIX ( R(I,J,K) > 

.REM 

280 

REM 

DATA 

0,0,0,00 


290 

REM 

DATA 

0,4000,0,0 


300 

REM 

DATA 

6000,0,0,0 


310 

REM 

DATA 

0,1000,1000,1000 


320 

REM 

DATA 

0,4000,0,0 


322 

REM 

DATA 

6000,0,0,0 


330 

REM 

DATA 

0,0,3000,3000 


340 

REM 

DATA 

0,4000,0,0 


350 

REM 

DATA 

6000,0,0,0 


360 

REM 

DATA 

0,0,0,1.E30 


370 

REM 

DATA 

0,1.E30,0,0 


380 

REM 

DATA 

6000,0,0,0 


381 

REM. 


COST/RETURN MATRIX, Q( I ,K )=(P ( I , J ,K ) *R ( I , J ,K ) 

.REM 

382 

DATA 

0,4000,6000 


383 

DATA 

1000,4000,6000 


384 

DATA 

3000,4000,6000 


385 

DATA 

1E30,1E30,6000 


390 

REM. 


TIE-BREAKER, TIED(I) 

.REM 

395 

DATA 

0,0,0, 0 



400 

REM. . 


INITIAL POLICY, D(I) 

.REM 

410 

DATA 




500 

REM. . 


END OF DATA LIST 

.REM 

998 

REM 



REM 


GENERAL INPUT 


REM 


1000 REM 

1005 READ Ki,Ko 

1010 REM ASSIGN READ/PRINT FILES REM 

1015 ASSIGN “OUTPUT" TO #Ko 
1020 IF Ki<=0 THEN 1190 
1025 ASSIGN -INPUT*' TO *Ki 
1030 READ *Ki>Iprint 
1035 READ *Ki>M 5 tate 
1040 READ ♦Ki>Ndecision 
1045 READ *l<i>Maxi 
1050 READ *Ki>Maxtrial 
1055 READ *Ki>Discount 

1060 REM INPUT TRANSITION MATRIX REM 

1065 FOR 1=1 TO Mstate 
1070 FOR K=1 TO Ndecision 
1075 FOR J=1 TO Mstate 
1080 READ ♦Ki>P(I,J^K) 

1085 NEXT J 
1090 NEXT K 
1095 NEXT I 

1100 REM INPUT COST/PROFIT MATRIX REM 

1105 FOR 1=1 TO Mstate 
1110 FOR K=1 TO Ndecision 
1115 REM FOR J=1 TO Mstate 
1120 REM READ ♦5iR<I,J;K) 

1125 REM NEXT J 
1130 READ #Ki>Q(I,K) 

1135 NEXT K 
1140 NEXT I 

1145 REM INPUT TIE-BREAKER FLAG REM 

1150 FOR 1=1 TO Mstate 
1155 READ *Ki>Tied(I) 

1160 NEXT I 

1165 REM INPUT INITIAL POLICY REM 

1170 FOR 1=1 TO Mstate 
1175 READ *5>D<I) 

1180 NEXT I 

1185 GOTO 1350 

1186 REM REM 

1190 READ Iprint 

1195 READ Mstate 
1200 READ Ndecision 
1205 READ Maxi 
1210 READ Maxtrial 
1215 READ Discount 

1220 REM INPUT TRANSITION MATRIX REM 

1225 FOR 1=1 TO Mstate 
1230 FOR K=1 TO Ndecision 
1235 FOR J=i TO Mstate 
1240 READ P(I,J,K) 

1245 NEXT J 
1250 NEXT K 
1255 NEXT I 


1260 REM INPUT COST/PROFIT MATRIX REM 

1265 FOR I=i TO Mstate 
1270 FOR K=i TO Ndecision 
1275 REM FOR J=i TO Mstate 

1280 REM READ R(I,J,K) 

1285 REM NEXT J 

1290 READ Q<I,K) 

1295 NEXT K 
1300 NEXT I 

1305 REM... INPUT TIE-BREAKER FLAG REM 

1310 FOR I«1 TO Mstate 
1315 READ Tied<I) 

1320 NEXT I 

1325 REM INPUT INITIAL POLICY REM 

1330 FOR 1=1 TO Mstate 
1335 READ D(I) 

1340 NEXT I 

1345 REM REM 

1350 FOR 1 = 1 TO Mstate 
1355 FOR K=1 TO Ndecision 
1360 IF Q<I,K>O0 THEN 1425 
1365 NEXT K 
1370 NEXT I 

1375 FOR 1=1 TO Mstate 

1380 FOR K=i TO Ndecision 

1385 Q<I,K)=0 

1390 FOR J=i TO Mstate 

1395 Q<I,K)=Q(I,K)+P(I,J,K)*R<I,J,K) 

1400 NEXT J 
1405 NEXT K 
1410 NEXT I 

1415 REM rEM 

1420 REM 

1425 GOSUB 10205 ! PRINT INPUT INFORMATION 

1430 REM 

1435 REM rEM 

1440 Discount=i/(i+Discount/i00 ) 

2000 REM ‘ rEM 

2032 IF lprint>0 THEN GOSUB 10410 

2034 Trial=0 

2035 Trial=Trial+l 

2036 IF lprint>0 THEN GOSUB 10510 

2050 REM VALUE DETERMINATION REM 

2060 FOR 1=1 TO Mstate 
2070 K=D<I) 

2090 FOR J=i TO Mstate 

2100 A(I, J)=-P(I, J,K))|cDiscount 

2110 IF I=J THEN A(I,J>=i+A(I,J) 

2122 NEXT J 

2124 IF Discount=i THEN A( I ,Mstate)=i 

2125 A<I,Mstate+i):=Q(I,K) 

2130 NEXT 'I 

2140 N=Mstate 

2141 Indic=i 

2142 Eps=lE-20 
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2iS0 IF Iprint>i THEN GOSUB i06i0 
2170 GOSUB 5000 
2180 IF Discount<l THEN 2300 
2190 G=X(hstate) 

2220 X<Mstate)=0 

2300 IF lprint>0 THEN GOSUB 10710 

3010 REM POLICY IMPROVEMENT REM 

3020 FOR 1=1 TO Mstate 

3025 FOR K=1 TO Ndecision 

3026 Sum(I,K>=0 

3030 FOR J=1 TO Mstate 

3040 Sum<I,K)=Su«(I,K)+P(I,J,K)*X( J) 

3041 NEXT J 

3042 SuM(I,K)=SuM(I,K)*Discount+Q(I>K) 

3054 IF Discount=l THEN Sum(I ,K )=Sum(I , K )-X( I ) 

3058 NEXT K 
3060 E=Sum(I,1> 

3070 Dd<I)=l 

3106 FOR K=2 TO Ndecision 
3110 IF Maxi>0 THEN 3150 
3120 IF E<=Sum<I,K) then 3170 
3130 E«Su«<I^K) 

3140 Dd(I)=K 
3142 GOTO 3200 

3150 IF E>=Sum(I,K) then 3170 
3160 GOTO 3130 

3170 IF Sum(I,K)OE THEN 3200 
3190 IF Tied(I)>0 THEN Dd<I>=K 
3200 NEXT K 

3202 IF lprint>0 THEN GOSUB 10810 
3310 NEXT I 

3320 FOR 1=1 TO Mstate 
3330 IF D<I>ODd(I) THEN 3360 
3340 NEXT I 
3350 GOTO 10905 

3360 IF lprint>0 THEN GOSUB ilOiO 
3362 FOR 1=1 TO Mstate 
3370 D<I)=Dd(I) 

3380 NEXT I 

3390 IF TriaKMaxtrial THEN 2035 
3410 GOTO lliiO 
5000 REM 

5005 REM 

5006 REM FUNCTION SIMUL ( N, A,X,EPS, INDIC.NRC ) 

5007 REM 

5008 REM INDIC=~1, COMPUTE THE INVERSE OF THE N X N MATRIX 

5009 REM INDIC= 0, THE SET OF EQUATIONS A(N, N)»X < N)=A< N+i > N+i) IS 

5010 REM SOLVED AND THE INVERSE IS COMPUTED 

5011 REM INDIC=+1, THE SET OF EQUATIONS A<N,N)*X<N)«A<N+i , N+i) IS 

5012 REM SOLVED BUT THE INVERSE IS NOT COMPUTED 

5013 REM 

5014 REM EPS =MINIMUM ALLOWABLE VAULE FOR A PIVOT ELEMENT 

5015 REM A =AUGHENTED MATRIX OF COEFFICIENT, A=A<I,J) 

5016 REM 
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5017 RCM N ^NUMBER OF ROUS IN A 

5018 REM 

5019 REM X =S0LUT10N VECTOR, X=X(I) 

5020 REM 

5021 Max=N 

5022 IF lndic>=0 THEN Max=N+l 


5023 REM 

5024 REM IS N LARGER THAN 50 

5025 IF N<=50 THEN 5031 

5026 PRINT *6>" N IS GREATER THAN 50" 

5027 Si«ul=0 

5028 RETURN 

5029 REM 

5030 REM BEGIN ELIMINATION PROCEDURE 

5031 Deter^l 

5032 FOR K=1 TO N 

5033 Km1=K-1 

5034 REM 

5035 REM SEARCH FOR THE PIVOT ELEMENT 

5036 Piuot=0 

5037 FOR 1=1 TO N 

5038 FOR J=1 TO N 

5039 REM 

5040 REM SCAN IROW AND JCOL ARRARYS FOR INVALID PIVOT SUBSCRIPTS 


5041 IF K=1 THEN 5048 

5042 FOR Iscan=l TO KmI 

5043 FOR Jscan=l TO KmI 

5044 IF I=Irow(Iscan) THEN 5052 

5045 IF J=Jcol(Jscan) THEN 5052 

5046 NEXT Jscan 

5047 NEXT Iscan 

5048 IF ABS(A(I,J) )<=ABS(Pivot) THEN 5052 

5049 Pivot=A(I,J) 

5050 Irow<K)=I 

5051 Jcol(K>=J 

5052 NEXT J 

5053 NEXT I 

5054 REM 


5055 REM INSURE THAT SELECTED PIVOT IS LARGER THAN EPS 

5056 IF ABS<Pivot)>Eps THEN 5062 

5057 PRINT *6>" ABS(PIVOT)=" > ABS(Pivot > > " IS LESS THEN ">Eps 

5058 S1mu1=0 

5059 RETURN 

5060 REM 

5061 REM UPDATE THE DETERMINANT VALUE 

5062 Irowk=Irow<K) 

5063 Jcolk=Jcol(K) 

5064 Deter=Deter)KPi VO t 

5065 REM 

5066 REM NORMALIZE PIVOT ROW ELEMENTS 


5067 FDR J=1 TO Max 

5068 A ( Ir owk , J >=A< Irowk , J) /Pi VO t 

5069 NEXT J 

5070 REM 


5071 REM CARRY OUT ELIMINATION AND DEVELOP INVERSE 

5072 Adrowk, Jcolk) = i/Pivo t 

5073 FOR I=i TO N 

5074 Ai jck=A<I, Jcolk) 

5075 IF I=Irowk THEN 5080 

5076 A<I;Jcolk)=-Ai jck/Pivot 

5077 FOR J=i TO Max 

5078 IF JOJcolk THEN A< I , J)=A( I , J)-Ai jck*A< Irowk , J) 

5079 NEXT J 

5080 NEXT I 

5081 NEXT K 

5082 REM 

5083 REM ORDER SOLUTION VALUES (IF ANY) AND GREAT JORD ARRAY 

5084 FOR I«1 TO N 

5085 Ir owi=Ir ow< I ) 

5086 Jcoli=Jcol(I> 

5087 Jord<Irowi)=Jcoli 

5088 IF lndic> = 0 THEN X < Jcol i)=:A( Ir owi ,Ma x ) 

5089 NEXT I 

5090 REM 

5091 REM ADJUST SIGN OF DETERMINANT 

5092 lch=0 

5093 Nm1=N-1 

5094 FOR I«1 TO NmI 

5095 Ipl=I+l 

5096 FOR J=Ipl TO N 

5097 IF Jord<J)>=Jord<I) THEN 5102 

5098 Jte«p=Jord ( J) 

5099 Jord<J)=Jord(I) 

5100 Jord(I)«TeMp 

5101 Ich=Ich+l 

5102 NEXT J 

5103 NEXT I 

5104 IF Ich/2*20Ich THEN Deter=-Deter 

5105 Si«ul=Deter 

5106 REM 

5107 REM END OF SUBROUTINE 

5108 RETURN 

5109 REM 
10100 
lOliO 
10120 
10130 
10140 
10200 
10205 
10210 
10220 
10225 
10230 
10235 
10240 
10250 
10255 


REM REM 

PRINT *Ko 
PRINT *Ko 
PRINT *Ko 
RETURN 

REM REM 

COMMAND "M F H Hp-Ib#l» 

PRINT ♦Ko>LIN(2)>TAB<10)>“5»c # 5K * INPUT INFORMATION 

PRINT ♦Ko>LIN(2)>TAB(10)>"MSTATE ">Mstate 

PRINT ♦KoiTAB<10))”NDECISION ”>Ndecision 

PRINT ♦Ko>TAB(10)>"MAXTRIAL ‘'^Maxtrial 

PRINT ♦Ko>TAB(10) > "MAXIMUM ">Maxi 

PRINT #Ko;TAB(10)> "DISCOUNT "^Discount 

PRINT #Ko>LIN(l) >TAB(10) ) "STATE"," "> "DECISION" /"TIE-BREAKER " 

PRINT ♦Ko>TAB(10)>" I D(I> TIED(I) " 


87 


;D(I),Tied<I> 


PRINT *Ko>TAB(iO> i " % " 

FOR 1=1 TO Mstate 
PRINT ♦Ko>TAB<iO))I/' 

NEXT I 

IF Iprint<2 THEN RETURN 

PRINT ♦Ko>LIN(l)>TAB(10)>" STATE 


DECISION 


PRINT ♦Ko>TAB(10)> " STATE % ” DE 
PRINT #KoiTAB(10))”AT CURRENT%“AT 
PRINT ♦Ko>TAB(10)>" STAGE 
PRINT *KoiTAB(iO)>” ( I ) 

PRINT *Ko>TAB(iO)>" 

FOR I=i TO Mstate 
PRINT *Ko 

FOR K=1 TO Ndecision 

PRINT #Ko 

FOR J=i TO Mstate 

PRINT ♦Ko>TAB(10) >I,“%K,J,P(I,J,K) 
NEXT J 
NEXT K 
NEXT I 

PRINT #Ko>LIN<2)>TAB(iO)>» STATE 


DECISION STATE » 
AT CURRENT”/' AT NEXT 
STAGE “ , ” STAGE 
< K ) < J ) 


STATE 

”, "PROBABILITY” 
P(I,J,I<) ” 
P(I,J,K) ” 


DECISION " 


PRINT ♦Ko>TAB(iO)>"AT CURRENT”, "AT CURRENT” 
PRINT ♦Ko)TAB(lO) /* STAGE STAGE 

PRINT #Ko>TAB(iO)/' ( I ) ",” ( K > "j 

PRINT #KOiTAB(i0) >” ”j 

FOR 1 = 1 TO Mstate 

FOR K=1 TO Ndecision 

PRINT ♦Ko>TAB(iO)>I,””,K ,Q<I,I<) 

NEXT K 
PRINT *Ko 
NEXT I 
RETURN 

REM 


”, "COST/PROFIT” 
Q(I,K) ” 


PRINT ♦Ko>LIN(2)iTAB(iO)>” ***** INTERMEDIATE RESULTS *****" 
RETURN 

REM REM 


PRINT *Ko 

PRINT #Ko>" TRIAL =”>Trial 

RETURN 

REM 


PRINT *Ko 

FOR 1=1 TO Mstate 

FOR J=1 TO Mstate+1 STEP 4 

PRINT ♦Ko>TAB(10)/'A<”>I>”,”>J/') = "iA(I,J)> ”A<"iI>”,”>J+l/*) = ”>A(I,J+l); 
PRINT ♦Ko/*A(">I/',”)J+2> ”) = ">A<I,J+2)>"A< ">!/•, ">J+3/') = ” >A( I, J+3) 

NEXT J 
PRINT #Ko 
NEXT I 
RETURN 

REM REM 


PRINT *Ko 

IF Discount=l THEN PRINT *Ko >TAB( 10 ) > ”G 
FOR 1=1 TO Mstate 


10740 PRINT #Ko>TAB(10)>”U(”>I>”)=”>X(I) 
10750 NEXT I 
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10760 PRINT *Ko 
10770 RETURN 

10800 REM REM 

10810 PRINT ♦Ko>LIN(l)>TAB(10)>“STATE%"DECISION“>”VALUE“ 

10820 PRINT ♦Ko>TAB(iO) >” % “ " 

10830 FOR K=1 TO Ndecision 

10840 PRINT ♦Ko>TAB(10);I,K,Su«(I,K) 

10850 NEXT K 
10860 RETURN 

10900 REM REM 

10905 PRINT #Ko>LIN(3)> 

10910 PRINT ♦Ko>TAB(i0)i"* )K )|c * * * )K FINAL RESULTS # * * “ 

10915 PRINT ♦Ko>LIN(3) 

10920 IF Discount«l THEN PRINT #Ko >TAB( 10 )> "LONG-TERM AVERAGE COST/RETURN=" > G 
10925 PRINT ♦KoiLIN(2>>TAB(i0)> “STATE",” "> "POLICY" , “VALUE” 

10930 PRINT ♦Ko>TAB(10) >” ” 

10935 FOR 1=1 TO Mstate 
10940 K=D<I) 

10945 PRINT tKo^TABdO))!," “ >D< I ) ,Sum( I ,K ) 

10950 NEXT I 
10955 GOTO 11210 

11000 REM REM 

ilOiO PRINT ♦Ko>LIN<2) 

11012 FOR L=1 TO Mstate 

11014 PRINT ♦Ko>TAB(10>>”OLD-D(”>L>“)=”>D(L),”NEW-D<”>L>”)=”>Dd(L) 

11016 NEXT L 
11020 RETURN 

11100 REM REM 

11110 PRINT *Ko >LIN(3) >TAB(20) > ” 

11120 PRINT ♦Ko>TAB(20) /'EXCEEDED THE MAXIMUM NUMBER OF TRIALS =”>Maxtrial 
11130 PRINT #Ko>TAB(20) > " 

11200 REM REM 

11210 PRINT *Ko>LIN(5)>TAB(20)>”* * 5K * )K END OF TASK * * * *" 

11220 COMMAND "M f h hp-Ib*l” 

60000 END 
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The Global Positioning System (GPSj is a worldwide navigation system using a con- 
stellation of earth satellites with onboard clocks. The GPS is also usable to transfer time 
and frequency. JPL has provided part of the funding to the National Bureau of Standards 
(NBS) to build GPS timing receivers. Preliminary tests with the breadboard receiver at 
NBS has produced remarkable precision in time and frequency transfer between the 
United States Naval Observatory (USNO) and NBS. JPL plans to install receivers in the 
DSN to demonstrate their ability to transfer time and frequency within the DSN and 
between the DSN and outside agencies. 


I. Introduction 

JPL is procuring two Global Positioning System (GPS) satel- 
lite timing receivers from the National Bureau of Standards 
(NBS). JPL has funded a large part of the development of the 
receiver design at NBS. The GPS timing receivers are to be evalu- 
ated for use in the DSN as a method of time and frequency 
transfer among the stations of the DSN, and between the DSN 
and world time and frequency standards. 

JPL plans to demonstrate the use of a GPS timing receiver 
system as a method of significantly reducing the cost of meeting 
the functional requirements for time and frequency transfer. 

II. The Global Positioning System 

The GPS is a Department of Defense worldwide navigation 
system that uses a constellation of orbiting Earth satellites, 
each of which carries an onboard clock (Ref. 1). Stored on each 


space vehicle (SV) is information on its position with respect 
to an Earth-based rectilinear coordinate system. If this time 
and position data were transmitted to an Earth-based station 
with its own clock, then the difference in time, neglecting 
relativistic effects, would be the propagation delay of the 
signal. Hence the observer can determine his position as being 
somewhere on a sphere which has its center at the spacecraft. 
By viewing three SVs and by having a reasonably good clock, 
an observer can determine his position exactly. By viewing 
four SVs, the observer can determine his position using a poor 
clock. 

The GPS satellite configuration calls for 24 SVs. This num- 
ber has been reduced to 18 due to budgetary considerations. 
The SVs are distributed in three orbital planes inclined 63° 
with respect to the Earth’s equatorial plane, and offset from 
each other by 120° longitude. Each SV is in 12-hour prograde 
orbit. This distribution is meant to maximize the number of 
SVs visible to an observer anywhere on the earth at any time. 
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The GPS SV signal is at two frequencies: link 1 at 1575.42 
MHz and link 2 at 1227.6 MHz. These frequencies are multi- 
ples of 10.23 MHz, namely: 

Hnkl = 1575.42 MHz = 10.23 MHz X 154 

link 2 = 1227.6 MHz = 10.23 MHz X 120. 

This is a particularly good frequency region for a couple of 
reasons. First, there is little use at these frequencies, so it is 
relatively easy to obtain an allocation. Second, the effects of 
the ionosphere on the propagation delay are smaller than at 
lower frequencies. 

The GPS SV transmits two PN codes on link 1 ; they are 
named the P code and the C/A code. The GPS timing receiver 
uses only the C/A code. The C/A code is 1023 bits long and 
is transmitted at a 1023-Mbps rate; hence its duration is 1 ms. 
Each SV’s code is unique so that all the SVs can use the same 
RF carrier frequency and yet a receiver can pick out a particu- 
lar SV’s signal. The data are transmitted at a 50-bps rate, so 
there are 20 C/A code sequences for each data bit. The data 
consist of 4 blocks. Data Block 1 contains the clock correc- 
tions for the SV clock with respect to GPS time. Data Block 2 
contains ephemeris information for that particular SV. There 
is a message block that contains any messages that might be 
needed. Finally Data Block 3 contains an almanac of the loca- 
tions of all of the SVs in the GPS constellation. The complete 
message takes 6 seconds to transmit. 


III. Description of the Receiver 

The NBS GPS timing receiver system consists of three 
chassis which are rack-mountable and a separate low-noise 
amplifier and down-converter which are mounted at the 
antenna. The rack-mountable chassis take a combined vertical 
rack space of about 18 inches, a keyboard which takes about 
4 inches more and a printer which takes about 6 inches. 

A 100-MHz local oscillator frequency is sent to the antenna- 
mounted down-converter. This signal is subsequently multi- 
plied (X15) and mixed with a 1.57542 GHz out of the wide 
band front end amplifier to produce a 75.42 MHz IF frequency. 

The balance of the receiver is what is known as a Tau Dither 
Delay Lock Receiver (Ref. 1). The receiver has two loops. The 
inner loop is a delay lock loop. Its function is to lock to the 
incoming PN code, thereby producing a signal the outer phase 
lock loop can lock to. 

The receiver must perform a cross correlation operation to 
extract a signal from a particular SV. To do this, the PN code 
of the particular SV is replicated in the receiver. Then the 


replicated code is shifted (delayed) until a correlation is 
found. With correlation, the receiver locks to the code with 
the delay lock loop and to the signal with a phase lock loop. 

IV. Timing and Frequency Transfer 

A. Synchronization 

One fundamental aspect of time is that of simultaneity.^ A 
definition given by NBS for simultaneity is: “Two events are 
simultaneous if equivalent signals, propagating in a given media, 
arrive coincidentally at a common point in space which is 
geometrically an equal distance from the source of each 
event. In practice a much broader definition is often used for 
clock synchronization; i.e., two clocks have the same reading 
in a specific reference frame” (Ref. 2). 

If two clocks are beating simultaneously, in a specific refer- 
ence frame, they are said to be synchronized. All of us have 
synchronized one clock to another by setting our wrist watch 
to the mantel clock. We make the time reading (date) of one 
dock agree with the date of another clock. 

If clocks are located some distance from one another, the 
problems of accurate synchronization become difficult. Of 
course we use radios and telephones to get the date, but these 
methods are not usable to transfer the date with an accuracy 
of microseconds or nanoseconds. 

Early in the 17th century Galileo devised a method of date 
transfer using observations of the moons of Jupiter. The Jovian 
satellite system is a clock that is observable from any place on 
the Earth. All one needed was an ephemeris table and a reason- 
able telescope and he could synchronize any clock on Earth. 
This is an application of the first definition of simultaneity. 
Note that every place on Earth is approximately the same dis- 
tance to Jupiter. This technique was used extensively for over 
100 years for map-making; however, it was not usable on ships 
at sea because of the difficulty of looking through a telescope 
from a pitching deck. 

With good portable cesium clocks the definition of sya- 
chronization of widely separated clocks is a “clock trip” from 
one clock to the other, then back to the first clock. This is 
called a clock closure. This method invokes the second defini- 
tion of synchronization. The most careful clock trips take into 
account the path of the traveling clock to correct for relativistic 
effects and refer the clocks to the same reference frame 


^ There are three aspects of time: time interval, date (the clock reading), 
and simultaneity. 
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(Ref. 3). This is one of the methods presently used in the 
DSN.2 

Over the last several years people have been getting back to 
Galileo’s method of using satellites for time transfer. For 
example NBS has been using the GOES weather satellites for 
time transfer. The DSN has used a technique of bouncing 
timing pulses off the Moon. Now both NBS and the United 
States Naval Observatory (USNO) are using the GPS satellites 
for time transfer. 

The GPS can be used in several ways to accurately transfer 
time between two or more known locations on the Earth. I 
shall describe two methods which are planned during the JPL/ 
NBS demonstration. 

The first of these two methods is a simultaneous common 
view (mutual view) of a single GPS satellite. With this tech- 
nique the common mode ephemeris errors cancel, and the satel- 
lite clock error contributes nothing. If the two or more clocks 
are within 3000 km of each other so that the observation angles 
are high, one can expect a measurement accuracy of 10 nano- 
seconds time difference between the observers’ clocks. If the 
propagation and ephemeris errors can be reduced, the accuracy 
can be reduced to 1 nanosecond. 

The specific procedure for a mutual view observation is 
quite simple. The observers agree to observe a given satellite at 
a given time. Alernatively, if an observer knows another’s 
observation schedule, he only has to make his schedule agree 
with the other’s. The data of a complete observation of a satel- 
lite would include the date of the observation to the second 
and the time offset from the observer’s clock to the GPS 
dock. At some time later the observers can exchange data. 

Ideally an observatory’s daily observation schedule would 
be on sidereal time.^ This would mean that the SV would 
appear in the same position in the sky every day at observa- 
tion time. The advantage of this method is the elimination of 
variables associated with direction. As the two clocks are sepa- 
rated, the angles of mutual view of the satellites become more 
shallow. As this happens, the ephemeris and propagation errors 
become larger, so the measuring accuracy is decreased. 


^In addition to clock closures, VLBI methods are now being used 
between the DSN stations. In Spain the Loran-C navigation system is 
used and in Australia television synchronization signals are used to 
transfer time. 

^Sidereal time is the time used by astronomers. It uses observations on 
the fixed stars rather than the Sun. A sidereal day is about 4 minutes 
shorter than a mean solar day; hence the time of observation would 
advance about 4 minutes per day mean solar time. 


The second method is to make successive observations of 
the same SV clock (sequential view). The first observer syn- 
chronizes his clock with the SV clock when the SV is at a good 
viewing angle. A second observer then synchronizes his clock 
to the same SV at some time which is less than 12 hours later 
than the first observation, and when the SV is at a good view- 
ing angle for him. If the SV that is observed has a cesium dock, 
then the date uncertainty will be about 5 ns over a 12-hour 
period. The ephemeris errors will tend to cancel because the 
same SV is used. 

As in the case of portable clock trips, the relativisitic effect 
needs to be taken into account for a precision time transfer 
(Ref. 3). 

B. Syntonization 

Closely related to synchronization is the concept of syntoni- 
zation. Syntonization is the act of putting two oscillators on 
the same frequency (tone). It is a musical term and is found in 
some dictionaries so defined. The term syntonization accuracy 
refers to a number (/-/o)//o = A///, which is a measure of 
how close the frequency of one oscillator is to that of another, 
where is the frequency of the reference oscillator and / is 
the frequency of the measured oscillator. 

As presently conceived, the operation of the GPS receiver 
would be used to take a reading of time offset once a day 
between two docks. If clock A’s oscillator is running fast 
with respect to clock B’s oscillator, then clock A’s value of 
date will advance with respect to clock B’s value of date. 
Therefore, by taking daily clock offset readings one can deter- 
mine the difference in rates between the two clocks oscillators, 
given A///. 


V. Present Performance of the Receivers 

At present NBS and USNO are conducting daily mutual 
view observations of SVs. This means that NBS acquires and 
receives data from the same SV at the same time USNO is 
doing so. NBS is using a breadboard version of a receiver of 
NBS design. USNO is using a receiver purchased from Stanford 
Telecommunications Incorporated (STI). 

USNO observes five SVs daily. The offsets between GPS 
time and USNO Master Clock (MC) are published weekly in a 
USNO Bulletin (Fig. 1). In addition, the data are available 
daily by telephone using a MODEM and a terminal. The data 
consist of the offset in nanoseconds between GPS time calcu- 
lated by each satellite and the time of the MC. Also the data 
includes the date to the second that the data were taken. 
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The USNO measures the offset between GPS time and MC. 
The MC is the physical realization of UTC (USNO)"^ and is 
designated UTC (USNO, MC). Because of the present method 
of steering the MC, the time scale UTC (USNO, MC) differs 
from UTC (USNO) by some small amount (several nanosec- 
onds). This offset is available on the telephone on a daily basis 
in the same way as the GPS data. 

NBS observes four SVs daily, scheduled so that it has a 
mutual view of USNO. The SVs are SV-5 and SV-9, which 
have cesium clocks on board, and SV-6 and SV-8 which have 
rubidium clocks. NBS measures the offset between GPS time 
by each SV and NBS clock 9 [(UTC (NBS, CL9))] . 

In the same way, UTC (USNO, MC) has some offset from 
UTC (USNO), UTC (NBS, C19) as some small offset from 
UTC (NBS). These data are not published but are available 
from NBS. 

There is a bit of a scheduling problem in the mutual view 
exercise between NBS and USNO. The NBS receiver corrects 
its viewing time 4 minutes a day to make up the difference 
between mean solar time and sidereal time. The USNO receiver 
on the other hand is corrected 28 minutes once a week. This 
means there is a slight departure from sidereal time as the 
week wears on. I understand that USNO will eventually operate 
the receiver on sidereal time. At the present NBS has to change 
the schedule every day. 

The big advantage of viewing the SV on a sidereal schedule 
is that it will appear in the same position in the sky every day. 
This is important if there are any multipath problems and with 
respect to syntonization. To measure frequency offset, the 
important thing is change in time. If we look at the SV in the 
same place in the sky every day, a lot of variables cancel. 

Present measurements between USNO and NBS have shown 
that the short-term characteristics of the GPS signal have white 
noise phase modulation (white PM) as a limiting noise process. 
The method of determining this was a new technique called 
the Modified Allan Variance (mod o^y(r)). Mod as its 

name implies is a modification of the traditional method of 
characterizing frequency stability, the Allan Variance. Using 
mod a2(r), one can easily distinguish between white PM and 
flicker phase noise (Flicker PM). Flicker PM is a type of power 
spectral density which is inversely proportional to the spectral 
frequency o?/27r. If Flicker PM noise is present then there is a 


^UTC universal coordinated time. The new UTC was adopted by the 
International Radio Consultive Committee (CCIR) in 1971 and put 
into effect 1 Jan. 1972. In this system all clocks run at the same rate, 
based on an atomic cesium clock. 


logarithmic divergence of the standard deviation which increases 
with the number of samples taken (Ref. 4). So it is important 
to determine if this noise process is present. 

Knowing the limiting noise process to be white PM, one can 
safely proceed using the statistics of the normal distribution. 
The results of this analysis have shown that one can measure 
time differences to less than one nanosecond for averaging 
times of 4 minutes (Ref. 5). The data from the two receivers 
are being analyzed by both NBS and JPL. Results have been 
published by NBS (Ref. 7). Figure 2 is a graph of UTC (USNO) - 
UTC (NBS) versus days. The graph covers about a four-month 
period. During this period several clock trips were made 
between USNO and NBS. The agreement is shown on the 
graph.^ 

The syntonization between the rates of UTC (USNO) and 
UTC (NBS) can be measured using the GPS timing receivers. 
If a set of time difference measurements are made using a 
single SV, then a slope can be found of time difference versus 
time, by finding a least square linear fit to the data. This 
assumes that the noise causing the variations in readings is 
white. This procedure is repeated for several SVs which pro- 
duces a set of slopes (Fig. 3a,b,c,d). The slopes are a measure 
in the difference in rates of the two clocks as measured by 
using different SVs. 

A second assumption is now made; that is, the measures of 
these slopes are independent. Next, the mean and standard 
deviation and the standard deviation of the mean are found 
for the set of slopes. 

For the period 5 Sept 81 to 16 Sept 81^ the slopes of 
[UTC (USNO)] — [UTC (NBS)] by the several SVs are: 

by SV 5 -11 .32 nanosec/day 
bySV9 -11.35 nanosec/day 
by SV 6 -10.73 nanosec/day 
by SV 8 -1 1.87 nanosec/day 

The mean is -1 1 .07 nsec/day 

The standard deviation is -0.67 nsec/day 

The standard deviation of the mean is -0.33 nsec/day 


^The measurement of the time offset via GPS has an offset of +445 
nanosec with respect to the clock closure. The reason is unknown at 
this time; it is probably due to design differences in the USNO and 
NBS receivers. The offset was removed from the GPS data. 

^These dates are arbitrarily taken. The period happens to be a set of 
continuous data. Other sets have produced better results. 
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The standard deviation of the mean is found by 
- -0.33 nsec/day 

where SD is the standard deviation and N is the number of 
degrees of freedom. Since it was assumed that the measure- 
ments via the different SVs were independent, this number 
is 4, the number of SVs. 

The value of 0.33 nanosecond per day corresponds to 


0.33 X 10*"^ sec/day 
86,400 sec/day 


-3.82 X 10'^^ 


M. 

f 


which is a measure of syntonization. 


VI. The GPS Timing Receivers in the DSN 

The GPS timing receivers will be deployed in the DSN in 
several steps. The first step will be to establish time and fre- 
quency transfer between NBS and the Time Standard Ensemble 
located at Goldstone. The final step will be to establish time 
and frequency transfer throughout the DSN. 

The first step of the demonstration will be to install the 
two GPS receivers at the Goldstone clock ensemble at Gold- 
stone. The receivers will be scheduled to perform mutual view 
observations with NBS and USNO. 

This configuration and procedure will repeat the demon- 
stration that is presently being performed between NBS and 
USNO. There will be a receiver at each of these two clock 
ensembles. The ensembles are located some distance from each 
other but not so great a distance that the mutual viewing 
angles are shallow. Finally there is a regular traveling clock trip 
between these clock ensembles. This demonstration will not 
only confirm the NBS-USNO experiment, but it will establish 
a time and frequency transfer between the DSN and NBS. The 
immediate results of this configuration and procedure will be 
to reduce the necessity of clock trips between the Goldstone 
clock ensemble and NBS from every 120 days to once a year 
or even less often. It is planned to permanently have one of 
the receivers at the Goldstone clock ensemble. Therefore this 
saving can be continually realized. 

Another saving will be realized when the Goldstone clock 
ensemble is moved to DSS 14. This move is plaimed for some- 
time in 1983. The presence of a GPS receiver will greatly 
help maintain time and frequency during the move. 


The reason for placing two GPS receivers at the same place 
is to check the receivers themselves and act as a control for the 
Goldstone-NBS experiment. The two receivers should give 
identical results if their antennas are in the same location. In 
addition this will give us a chance to check the receiver before 
it is sent to another station in the DSN. 


I plan to repeat the procedures that are now being done 
between NBS and USNO. The difference is time and frequency 
transfer from NBS to the Goldstone clock ensemble rather than 
USNO. In addition time comparisons will be made between 
the two GPS receivers which are located next to one another 
at the Goldstone clock ensemble. These measurements should 
be identical except for any errors in the receivers themselves. 


Beyond repeating the NBS to USNO demonstration and the 
establishment of time and frequency transfer from NBS to the 
Goldstone clock ensemble, an even more informative measure- 
ment will be available. If the NBS, USNO and the Goldstone 
dock ensemble GPS receivers are all on the same schedule, for 
mutual view, then a three-station measurement closure will be 
available. This will be somewhat akin to the three baseline 
VLBI closures reported in Ref. 8 except that all of the stations 
will be located in the continental U.S. 

This will allow a simultaneous measurement of the time 
offset [UTC (USNO) - UTC (DSN)] , [UTC (USNO) - UTC 
(NBS)] and [UTC (USNO) - UTC (DSN)] . There wiU be two 
checks of the accuracy of time transfer using the GPS receivers. 
First the offsets [UTC (USNO) - UTC (NBS)] and [UTC 
(NBS) — UTC (DSN)] are well-known and kept with traveling 
clocks.'^ Second, by the closure measurement it should be pos- 
sible to get a measure of the absolute accuracy of the GPS 
receiver system independent of clock trips. 


The second step will be to move one of the GPS timing 
receivers to DSS 63 in Spain. Spain was chosen as the second 
installation because it is closer to Goldstone GTS than DSS 43 
in Australia. The great circle distance between GTS and Madrid 
is approximately 85° whereas the great circle distance between 
Goldstone and Canberra is approximately 110°. 

It is also possible to complete a closure of measurements 
which includes DSS 63. USNO is approximately 55° great 


'^AU times have been subtracted from UTC (USNO). This corresponds 
to using USNO as the Standard. This is done for convenience. JPL 
traditionally uses NBS as the Standard for the DSN because of its 
proximity. 
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circle distance from Madrid. There will be mutual view mea- 
surements (once per day, it is hoped) among the four locations 
USNO, NBS, Goldstone and DSS 63. This daily mutual view 
might involve some low viewing angles from DSS 63 because 
of nonconformance to USNO’s schedule. Daily mutual or 
sequential view between Goldstone and DSS 63, which is 
independent of outside agencies, is being planned. 

There are plans by NBS to install a receiver in Paris, France, 
at the Bureau International de THeure (BIH), sometime in 
1982. Paris is only 10° great circle distance from Madrid; 
therefore one could expect the same performance as between 
USNO and NBS. This assumes a daily high angle observation 
from BIH. 

The third step will be to install a receiver at DSS 43. This 
is the most difficult step because of the distance. The great 
circle distance from GTS to DSS 43 is approximately 110° and 
from DSS 43 to DSS 63 approximately 160°. Mutual views at 
these distances are right above the horizon; therefore a schedule 
of sequential views will probably be the more useful. Because 
JPL will have only two receivers we are planning to use one of 
NBS’s receivers for the DSS 43 demonstration. 

An exciting experiment will be a three-station closure of 
DSS 43, DSS 63 and GTS. The results of this will be com- 
pared to the VLBI three station closure in 1979 (Ref. 7). 
The DSN receivers v^ill be scheduled to make a three-station 
closure every day if possible. Certainly we will make an off- 
set measurement, [UTC (DSN) — UTC (DSS 43)] and [UTC 
(DSN) - UTC (DSS 63)] , every day. 


VII. Conclusion 

It is no mean task to determine the accuracy of time trans- 
fer over intercontinental distances using the GPS receivers. 
Suppose a sequential view of a SV is taken by two receivers 
anyplace on the Earth. Furthermore, suppose the SV is equipped 
with a cesium clock. Then time transfer accuracies of 10 nsec 
to 50 nsec are anticipated (Ref. 8). This means for a three- 
station closure the expected error is \/J times the individual 
experiment error or, 17 - 86 nsec. By using multireceiver 
closures, and repeating the short distance (<4000 km), trans- 
fers which were demonstrated in the continental U.S., the 
accuracies can certainly be confirmed. 

The obvious way to determine the accuracy of GPS receiver 
time transfer is to make an independent measurement using 
established methods. Data on the regular clock closures give 
errors of 300 nsec, which is too large to measure GPS receiver 
performance. A careful clock closure could yield an error of 
35 nsec for a trip of <40 hours each way, with three such trips 
2 months apart. Frequency could be measured with an error of 
1 X 10-14 

The other method of time transfer used in the DSN is VLBI . 
Present VLBI errors are approximately 100 nsec from DSS 14 
to DSS 63, and 200 nsec from DSS 14 to DSS 43. The DSS 43 
to DSS 63 baseline measurement is not done often but the 
accuracy of that baseline is approximately 300 nsec. Within 
the next few months, there are plans to implement new tech- 
niques that are expected to decrease the errors from one to 
two orders of magnitude (Refs. 9, 10). This would mean that 
VLBI could be used as an independent measurement of the 
time offsets between the clocks in the DSN. 
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GLOBAL POSITIONING SYSTEM (GPS) 





VALUES PRESENTED BELOW FOR NAVSTAR GPS SATELLITES ARE THE RESULT OF A LINEAR FIT THROUGH 

APPROXIMATELY 100 DATA POINTS, REFERRED TO THE BEGINNING OF THE TRACKING PERIOD (TIME 

OF MEASUREMENT). 

SATELLITES ARE TRACKED FOR APPROXIMATELY TEN MINUTES. 




NAVSTAR 1 

NAVSTAR 3 

NAVSTAR 4 

NAVSTAR 5 

NAVSTAR 6 



SV#4 

SV#6 

syo 8 

SV*5 

SV'9 


MJD 

MC-GPS UT 

MC-GPS UT 

MC-GPS UT 

MC-GPS UT 

MC-GPS UT 

NOV. 9 

44917* 

_ 


_ 


_ 

10 

44918 

-49.805(180448) 

-49.877(175330) 

-49.879(174224) 

-49.866(191530) 

-49.867(190530) 

11 

44919 

-49.697(180430) 

-49.805(175330) 

-49.811(174230) 

-49.798(191548) 

-49.792(190430) 

12 

44920 

-49.705(180442) 

-49.743(175354) 

-49.755(174230) 

-49.737(191548) 

-49.732(190430) 

13 

44921 

-49.621(180448) 

-49.677(175330) 

-49.688(174230) 

-49.673(191554) 

-49.676(190448) 

14 

44922 

-49.538(180430) 

-49.592(175330) 

-49.607(174230) 

-49.599(191530) 

-49.595(190430) 

15 

44923 

-49.419(180424) 

-49.544(175400) 

-49.545(174218) 

-49.534(191548) 

-49.546(190424) 

16 

44924 

-49.434(180454) 

-49.475(175348) 

-49.480(174218) 

-49.472(184830) 

-49.458(183818) 

17 

44925 

-49.306(173730) 

-49.404(172630) 

-49.412(171518) 

-49.381(184830) 

-49.398(183730) 

18 

44926 

-49.253(173730) 

-49.352(172654) 

-49.353(171530) 

-49.775(184830) 

-49.012(183730) 

*NO DATA DUE TO RECEIVER MALFUNCTION 





Fig. 1. Time differences between USNO master clock minus GPS time 
(from a U.S. Naval Observatory publication) 
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Noise Temperature and Noise Figure 
Concepts: DC to Light 

C. T. Steizried 
TDA Mission Support Office 


Deep Space communication systems require extremely sensitive receiving systems. The 
NASA Deep Space Network is investigating the use of higher operational frequencies for 
improved performance. Noise temperature and noise figure concepts are used to describe 
the noise performance of these receiving systems. It is proposed to modify present poise 
temperature definitions for linear amplifiers so they will be valid over the range 
(hf/kT) « 1 « (hf/kT). This is important for systems operating at high frequencies and 
low noise temperatures, or systems requiring very accurate calibrations. The suggested 
definitions are such that for an *Hdear amplifier, T^= (hf!k)= T^ and F= 1. These 
definitions revert to the present definition for (hf/kT) «L Noise temperature 
calibrations are illustrated with a detailed example. These concepts are applied to system 
signal- to-noise analysis. The fundamental limit to a receiving system sensitivity is 
determined by the thermal noise of the source and the quantum noise limit of the 
receiver. The sensitivity of a receiving system consisting of an **ideaV" linear amplifier 
with a 2.7 K source, (-194.3 dBm/Hz assuming (hf/kT )«1) degrades significantly at 
higher frequencies. 


i. Introduction 


II. Theory 


Deep space communication systems require extremely 
sensitive receiving systems (Ref. 1). The Deep Space Network 
is investigating the use of higher operational frequencies for 
improved performance. Noise temperature and noise figure 
concepts are used to describe the noise performance of these 
receiving systems (Refs. 2, 3, 4). These concepts are reviewed 
for application to higher frequencies by removal of the low- 
frequency restriction, (JifjkT)« 1. 


The available thermal (Refs. 5, 6, 7) noise power from a 
source at the amplifier output (G » 1) is given by 


P = kTBG 

n 


hfIkT 

e'’r/kT_i 


) 


( 1 ) 


100 


where 


T = source temperature, K 
h = Planck’s constant = 6.6256 X 10”^^ J-S 
k = Boltzman’s constant = 1.3805 X 10“^^ J/K 
f - operating frequency, Hz 

B - ^ I = noise bandwidth, Hz 

^ Jo 

G(f) = available power gain, ratio 

G = maximum available power gain, ratio 

The amplifier output is approximately (disregarding the 
contribution of the amplifier: hf « kty 


These correction terms are shown plotted in Fig. 2. For 
example, at 32 GHz, the cosmic background temperature 
(Refs. 9, 10) of = 2.7 K, correctly defined for use with Eq. (1) 
is “corrected” to ^ 2.0 K for use with Eq. (2). 

The quantum noise limit (Ref. 7, 11) of a linear amplifier 
(where ‘both phase and amplitude information is retained), a 
manifestation of the quantum mechanics uncertainty princi- 
ple, is given by^ 

= hfBG (5) 

Fundamental limits of an ideal receiving system sensitiv- 
ity are determined by the sum of the source thermal noise and 
the quantum noise limit of an ideal amplifier (since these noise 
sources are uncorrelated (Ref. 7)). 


= kTBG (2) 


Eqs. (1) and (2) are shown plotted in Fig. 1 for a large 
range of hfIkT values. Most microwave applications with 
QifIkT) « 1 are restricted to the region near the origin. 

or 

It is computationally convenient to define a temperature T* 
such that {PjkBG) = + 


P = kTBG +hfBG = kiT'+T )BG 

n \^hflkT_^f <l' 


( 6 ) 


p 


n 


= ktBG = kTBG 


( 


hflkT \ 


T* can be conveniently found by subtracting a correction 
(Ref. 8) from T 


This is plotted if Figs. 1, 3 and 4 as functions of QiflkT)yf and 
r. In Fig. 4, the value of PjkBG for {hf!kT» 1) is given by 
hfik (or TJ and for (hflkT« 1), by T+ /i//2/: (or T+ Tjl), 
It is easily shown (Fig. 1) that the quantum noise limit and 
thermal noise are equal when (hfIkT) = In 2 = 0.69. 


(PjkBG) ^ r = 


where 


Tc = 



hfIkT \ 

ehfIkT.J 


III. Noise Temperature and Noise Figure: 
(hf/kT) « 1 

A receiving system noise performance is characterized by 
the operating noise temperature (Refs. 2, 3, 4, 13) defined by 
(linear amplifier, single channel, matched source, G » 1) 


or conveniently 

T^ s 0.024/(GHz) - 0.000192 |/(GHz)]2/r + • • • (4) 


‘Note that (.hfIkT) « 0.048 f (GHz)/r(Ar) « 0.00048 at 1 GHz and 
100 K and » 1 at 208 GHz and 10 K indicates that Eq. (2) is an 
extremely good approximation for most microwave applications. 


IkBG (7) 

^ o 


^Eq. (5) is appropriate only for linear amplification. At optical 
frequencies, using discrete photons, techniques may exist (Ref. 7) to 
circumvent this limitation. Equating kTB to quantum noise hfB results 
in an equivalent quantum noise temperature, = (hf/k). Although 
is a fictitious temperature, it is useful for computational analysis 
and can be used with Eq. 2 to compute 
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where 


Substituting T\ for in Eq. (8)^ 



receiver system total output noise power, within 
the frequency band 5, excluding output load noise, 


W 


op 


= 

I e 


( 12 ) 


Also 

Top = Ti + T^ ( 8 ) 

where 

= input source noise temperature, K 
= effective input noise temperature of the receiver, K 

The definition (Refs. 2, 3, 4) of an amplifier noise figure is 
F = N^ (T. = TjIkTFG (9) 

O 

where 

(T. = T^) = receiver total output noise power with 
^ input termination at temperature W 

T = 290K 

o 


where 



Following the spirit (F = 1 for ah “ideal” amplifier) of Eq. 
(9), define 

F = N^ (r, = rj/k (r; + TJ BG (13) 

o ^ 

where 

(T. = = receiver total output noise power with 

^ input termination at temperature W 



T = 290K 

o 


From Eqs. (7), (8) and (9) 


From Eqs. (1 1), (12) and (13) 


F = \+{TjT^) 


( 10 ) 


F = (i+r/r)/(i+VO 


(14) 


or 


or 


r = (F- \)T 


r = (F- i)r' +FT 

e ^ ' o q 


For an “ideal” amplifier, = 0, F= 1, and 


For an “ideal” amplifier = 1, and 


^ op 


r„„ = r.'+r 

op i q 


IV. Noise Temperature and Noise Figure: 
0 < (hf/kT) < 00 

As before. 


These equations all revert to present definitions. Section III, 
when (hfjkT) « 1 ; i.e., = 0, 7/= 7}, = 7^. 




(11) 


3 

is obtained by analysis or measurement as shown in Appendix A. 
For an ideal amplifier, = hfjk = 
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V. Noise Temperature and Noise Figure: 
(hf/kT)»-i 

As before, 


where 


S. = input signal power, W 


o 


( 15 ) FoT(0<(hflkT)<co), 


At these very high frequencies, thermal noise is negligible 
compared to so that T/= 0, and from Eq. (12) 


(W = ■ 


= s,ik(r;t T^)B 


op 


= T, 


(16) 


= Sjk (t ; + (f - 1) r; + Frp B (20) 


and from Eq. (13), 


F = Nj, IkT^BG 
= IhfBG 

O 


andfor(/i//A:70»l, 


(17) 


(SjNo) = SjikT^^B) 


= SjkT B 

I e 


From Eqs. (15), (16), and (17), or directly from Eq. (14), 

F = TJT^ (18) 


= SjFkT^B 

= SjFhfB (21) 


or 


T = FT 

e q 

For an “ideal” amplifier, F = 1 , and 


op 



The performance of a receiving system composed of an 
“ideal” amplifier (T^ = T^) and a source at the cosmic 
background temperature {T^ = 2.7 K), is shown in Table 1 and 
Fig. 5. is the input signal power required for (SJN^) = 1. 
Table 1 and Fig. 5 demonstrate the loss in sensitivity at very 
high frequencies relative to low frequencies for a conventional 
receiving system with an “ideal” linear receiver. 


VI. System Performance 

The output signal-to-noise power ratio for a receiving 
system is given by (QtfIkT) « 1), 

(W = 

= Sjk(Tj + T^)B 

= Sjk(T, + (F-l)TjB (19) 


VII. Conclusion 

The equations developed in the previous sections are 
tabulated in Table 2. The use of noise temperature and noise 
figure concepts require special consideration when 
QifIkT) « 1. This is usually important for systems operat- 
ing at high frequencies with low noise temperatures, or sys- 
tems requiring precise calibrations. 

Fundamental limits to a receiving system are determined by 
the thermal noise of the source and the quantum noise limit of 
the amplifier. The sensitivity of a receiving system consisting 
of an “ideal” linear amplifier with a 2.7 K source (-194.3 
dBm/Hz assuming QifjkT) « 1) degrades significantly at 
higher frequencies. Of course, there are tremendous advantages 
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at higher frequencies, such as greater potential bandwidth, less 
spectrum crowding, link security, smaller antennas, photon 
counting schemes, etc. 

A frequent end goal of using the more complicated 
equations of Section IV is to obtain the most accurate value of 
for use with system performance evaluation such as in 
Section VI. consists of the sum of r/and 7^. The decrease 


of 7/ and increase of 7^ with frequency tends to cancel, 
minimizing the error in 7^^. Therefore, the techniques of 
Section III can be used to very good accuracy to extremely 
high frequencies (see Table A-1 and Fig. A-2) for calculation 


However, for precise calibrations, the proposed equations 
of Section IV should be used. 
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Table 1. Tabulationof the sensitivity of an “ideal'’ receiver (Tg = Tq) with an input source 
temperature of 2.7 K as a function of frequency 


Equations 

Equations 

Equations 

(8). (19) 

(12), (20) 

(16), (21) 

Frequency, GHz 


Parameter 

0.0 

32 

200 

400 

2000 

20,000 
(\ = 15 Mm ) 

T\,Yi 

2.7 

2.0 

0.3 

0.0 

0.0 

0.0 

7’e = 7’,.K 

0.0 

1.5 

9.6 

19.2 

96.0 

960.0 


2.7 

3,5 

9.9 

19.2 

96.0 

960.0 

dBm/Hz 

-194.3 

-193,1 

-189.0 

-185.8 

-178.8 

-168.8 


Table 2. Summary of noise temperature and noise figure concepts: dc to light (Linear amplifier, single channel, matched source, G > > 1) 

(hflkT) = TqlT 

Parameter 0 < (hflkD < ~ QiflkT) « 1 QiflkT) » 1 


T'op=N^JkBG 

= T. + T 
1 e 

II 

= *[?’;+ T^] 

= 

II 

F = Nj, (r = rj/k (r; + tj bg 

= (1 + Tjr^Kx + r^/r;) 

= (T^ = T^)lkT^BG 

o 

= i + (TglT^) 

= N IhfBG =Nj, IkT BG 
‘o '' 

= i+(W 

= *[i] 

= *[1] 

= •[!] 

= S/k{T'. + T^)B 

= S.lk{T. + T^)B 

= S.lkT^B 

= S.lk(Tl + (F - 1) + FT^) B 

= SJk(T. + {F-\)T^)B 

= S.lFhfB=SjFkT^B 

= *[S^lk(T‘+T^)B] 

= *[SjkT^B] 

= *[S^JhfB=S.!kT^B] 

*[ ] “Ideal” ampUfier (i.e., 7^ = 

' 'Lw-J 


S, = Input signal power, W 

. 


Nj, = Total output noise power within 
o bandwidth W 

n ■ rj 'I 



T = 290 K 
o 



= W/k) s 0.048/ (GHz) 
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(P ATBG) 



Fig. 1. (Pp IkTBG) vs (hf/kT), showing the thermal and 
quantum noise contributions 



Fig. 2. Thermal noise temperature correction Tq as a function of 
frequency and temperature 



Fig. 3. (PnIkTBG) vs frequency for various temperatures, showing 
thermal and quantum noise contributions 
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Appendix A 

Amplifier Noise Temperature Caiibration 


Consider the Y factor equations for the amplifier in the test 
configuration shown in Fig. A-1 using two matched calibration 
terminations at physical temperatures Tj and Tj. 

Casel:(fif/kT)«l 


Y = 




(A-1) 


Solving 


T 

e 


YT. 


Y- 1 


(A-2) 


where 





Consider the following parameters applicable to Fig. A-1. 
Tj = 100 K = 10 K 


Case2:0<(;f//A:r)<«^ 


= 400K Y = 3.7 


Y = 




Solving 


+ YT[ 
Y- 1 


(A-3) These values are used with Eqs. (A-2) and (A4) for a wide 
range of frequencies. The results are tabulated in Table A-1 
and plotted in Fig. A-2. T'. and can be estimated from 
Fig. 3. For these parameters, decreases with frequency, 
while increases. Since is the sum of T/and there is 
a compensatory effect; however, there is a resultant net 
. . increase in with frequency. For these parameters, the 

^ ' amplifier is nearly “ideal” at 400 GHz. 
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Table A-1, Tabulation of receiver and system parameters, using 
test parameters: Ti = 100 K, T 2 = 400 K, Tj = 10K and Y = 3.7 



Equations 

(A-l).(8) 


Equations 

(A-2),(12) 




Frequency GHz 


Parameter, K 

0.00 

32 

200 

400 

7 ; 

100 

99.23 

95.28 

90.71 

r- 

400 

399.23 

395.22 

390.48 

r. 

10 

9.25 

5.96 

3,30 

= wik) 

0.0 

1.54 

9.60 

19.20 

T'e 

11.11 

11.88 

15.81 

20.32 

T 

op 

21.11 

21.13 

21.77 

23.62 
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TIGER Reliability Analysis in the DSN 

J. M. Gunn 

Deep Space Network Support Section 


TIGER is a computer program designed to simulate a system over a period of time to 
evaluate system reliability and availability. Results can be used in the Deep Space Net- 
work for initial spares provisioning and system evaluation. This article describes the 
TIGER algorithm, the inputs to the program and the output. 


I. Introduction 

Efficient spares provisioning is a continuing problem in the 
Deep Space Network (DSN). Correct determination of equip- 
ment to be spared, spares stock and optimal network alloca- 
tion becomes increasingly critical as budgetary restrictions 
reduce the margin of error. The Deep Space Network Support 
Section (377) has modified a computer program developed by 
the Department of the Navy to evaluate system reliability and 
availability which can aid in this process. Known as TIGER, 
this program’s applications in the DSN include initial spares 
provisioning, system reliability and availability evaluation and 
maintenance planning. A major asset is the program’s ability 
to model a wide range of systems. Systems as complex as an 
entire tracking station to the simplest station subsystem can be 
accommodated. Models can be simulated on the network, 
complex, station, system and subsystem level. 

Presently the DSN relies on the engineering judgment of 
the Cognizant Design Engineer and Cognizant Operations Engi- 
neer, aided with results from the DSN Efficient Sparing Pro- 
gram, to provision spares for the Deep Space Stations in its 
tracking network. By taking into account factors such as failure 
rates, repair rates and shipping time, TIGER can provide the 
engineers with an analytic tool that takes into consideration 
all major system parameters. 


II. Simulation Theory 

TIGER uses Monte Carlo simulation methods to model the 
system. This procedure assumes a constant failure rate over a 
period of time for the system equipment. The probability of a 
failure occurring before time t can therefore be determined by 
the exponential distribution function F{() = 1 - where/, 
is the failure rate and t is the time to failure (TTF). Substitut- 
ing the reciprocal of the mean time between failure (MTBF) 
for L and solving for f, the equation used in the TIGER algo- 
rithm is derived: TTF = -{\jMTBF) In [1 - F’(O) * To obtain 
an equipment TTF a random number between zero and one is 
generated and substituted into the equation for [1 - F(r)] 

(Fig- !)• 

Similarly, equipment time to repair (TTR) is simulated by 
drav^ng from an exponential distribution with the mean equal 
to the mean time to repair (MTTR). Here the random num- 
ber is substituted into the equation TTR=-{\IMTTR) In 
[1 - F(/)] to generate the time of equipment repair. 

The system is simulated over a meaningful period of operat- 
ing time, referred to as a mission. System up and down times 
are determined by generating a TTF for each lowest replace- 
able element (LRE) in the system. Initially each LRE is in an 
upstate. To detect system failure before completion of the 
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operational period simulated, the system is tested as each LRE 
fails. With each failure a TTR is generated for the LRE and a 
spare, if available, is put in its place. A new TTF is generated 
for each LRE experiencing repair to allow further failure after 
repair. This procedure continues until the mission aborts or 
the specified mission period is exceeded. The mission is aborted 
if the specified allowable downtime is exceeded due to lack of 
a critical equipment or subsystem. A new set of randomly gen- 
erated values is used to determine equipment TTF and the pro- 
gram is rerun. At the end of each set of 50 mission simulations, 
reliability and availability figures are computed. Running a 
large number of missions (typically 500-1000) and averaging 
the results gives a high degree of accuracy. 

TIGER is designed to simulate systems with variations in 
configuration. The operational sequence simulated is comprised 
of one or more operational phases. Each phase represents a 
unique configuration for the system. A maximum of 95 phases 
of 6 phase types are allowed. The configuration is based on a 
reliability diagram of the system. The system is divided into as 
many levels of subsystems as appropriate vrith the lowest level 
being that for whfch the MTBFs andMTTRs are supplied. From 
a maintenance standpoint this is the LRE. 


III. Program Features 

Much flexibility is provided by the program. Standby equip- 
ment may be designated to come up in the event of equipment 
failure. Operating rules may also be specified to cause a desig- 
nated string of equipment to go down if one of the string 
equipment fails and no spares are available. These rules may 
vary by phase or remain constant throughout the mission. 

Duty cycle, or percent of time used, may also be specified 
for each LRE by phase. This option is useful in systems where 
some equipment is used less rigorously than other system 
equipment. In the DSN this might apply to subsystems that are 
used less heavily during nontracking periods. Variable MTTRs 
may also be specified to provide for situations in which an 
equipment cannot be repaired during a particular mode of 
operation. 

A sophisticated, lo^sitics system is incorporated into 
TIGER. Numbers and types of spares may be specified as well 
as at what level the spares will be stocked. Three levels of 
sparing and associated administrative delay times are provided 
for. In the DSN this generally corresponds to the station, com- 
plex and network levels. 

System and subsystem maximum allowable downtimes are 
user-specified to define mission abort conditions. Both sus- 
tained and total downtimes are considered. 


User inputs to TIGER are equipment MTBF,MTTRand sys- 
tem configuration. Steps to using the program are as follows: 

Define operational period to be simulated and allowable 
downtime. 

Determine subsystem reliability configuration, noting 
parallel and series equipment. 

Determine LREs. 

Break system down into subsystems as warranted. 

Gather MTBF and MTTR data for LREs 

Determine sparing locations and associated logistics delay 
times. 

Determine program options desired 
Code input data in TIGER acceptable form. 

Run program. 


IV. Program Report 

In the final report estimators for total system reliability, 
instantaneous availability, average availability and system 
readiness are given. In addition, mean uptime, mean downtime 
and mean time between mission failures are tabulated. 

TIGER uses the generally accepted definition of relia- 
bility: the probability that the system will perform its in- 
tended function for a specified interval under the stated condi- 
tions. Mathematically this is calculated by dividing the number 
of mission aborts by the total number of simulated missions. 

Two availability estimators are calculated. The average 
availability is the probability the system will be in satisfactory 
operating condition at a random point in time. In the DSN 
this is functional availability, or the fraction of scheduled ser- 
vice time the system performs its intended function. This is 
calculated mathematically as : 

Average AvaOabmty 

The instantaneous availability is the same probability for a 
specific point in time (i.e., the beginning and the end of the 
phase.) 

System readiness is the probability the system will be in 
satisfactory operating condition when there is neither a mis- 
sion abort nor a system failure. The readiness estimator is 
generally considered the lower bound on the availability 
estimator. 
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Most useful from a maintenance standpoint are the indi- 
vidual equipment unreliability and unavailability calculations. 
This part of the report ranks the equipment in order of contri- 
bution to system unreliability and unavailability. This isolates 
the major contributors to system downtime and calculates the 
individual equipment percent contribution to unreliability. 

Spares usage is reported by LRE. Average spares used per 
mission and the maximum amount of spares used in a single 
mission are given. For initial spares provisioning, an unlimited 
spares stock may be specified. In this case the results would 
reflect the optimal spares stock for 100% system reliability if 
it is obtainable with the specified system components. 

A more detailed report, if desired, gives a complete listing 
of system equipment events. Up and down times of all critical 
LRE are listed in chronological order. 


V. Example 

To demonstrate the usefulness of the TIGER analysis the 
program results for the command modulator assembly (CMA) 
are included. The system has 34 LREs of 13 types. The system 
consists of two CM As, one of which must be operational to 
prevent system downtime (Fig. 2). CMA 801 will come up only 
upon the failure of CMA 800. 

The first part of the report (Fig. 3a,b) lists pertinent sys- 
tem parameters. Two simulation report options are available. 
The management summary prints a message each time a mis- 
sion abort occurs (Fig. 4). The phase type, phase sequence, 
mission number and time of mission abort are listed. The engi- 
neering summary is a more complete report that indicates each 
time an equipment changes status (up to down or down to up) 
or a phase boundary is crossed (Fig. 5). 

After each set of 50 system simulations a “figures of merit 
summary” is printed. Reliability and availability figures are 
calculated for each phase of the mission and for the entire mis- 
sion timeline. The system report (Fig. 6) is for 250 missions. 


The equipment failures and corrective maintenance summary 
(Fig. 7) tabulates for each LRE the number of failures for all 
missions. The average number of failures per mission and the 
average number of corrective maintenance hours expended are 
listed. 

The spares report lists the average number of spares used in 
any single mission. The report for the CMA (Fig. 8) indicates 
that the spares stock could be more optimally allocated since 
the entire stock is not used. 

The critical equipment list ranks each LRE in order of 
greatest contribution to system unavailability and unreliability 
(Figs. 9 and 10). The results for the CMA indicate the subcar- 
rier synthesizer (SC SYN) is a major contributor to both sys- 
tem unavailability and unreliability. Attention should be given 
to improving the design of the SC SYN to increase the MTBF 
and to procuring SC SYN spares. 

VI. Conclusions 

Used in combination with the DSN Efficient Sparing Pro- 
gram, TIGER can potentially aid the cognizant engineers in 
both initial and ongoing spares provisioning and improve 
operational availability of system equipment. 

The simulation techniques used in TIGER are in wide use 
today for system reliability and availability evaluation. A sec- 
ond version of TIGER called TIGER/MANNING will soon be 
completed. This version will include the simulation process of 
the original program augmented with calculations of optimal 
worker allocation in system maintenance and repair duties. 

The only significant disadvantage, common to all simula- 
tions of this type, is the considerable amount of computer time 
and storage required when a large complex system is modeled. 
However, most foreseeable applications of this program within 
the DSN are not on such a scale. Another minor problem is the 
inevitable random errors found in all system simulations. In 
actual use the major problem has been a lack of reliable failure 
data for system equipment. The results, of course, are limited 
by the accuracy of the parameters supplied. 


114 



























SIMULATION CONTROL SPECIFICATIONS 

STATISTICAL PARAMETERS USED FOR THIS SIMULATION: 

250 - MAXIMUM NUMBER OF MISSIONS TO BE RUN 
250 - OPTIMUM NUMBER OF MISSIONS TO BE RUN 
.998 - SPECIFICATION REQUIREMENT FOR RELIABILITY 
1 .280 - STANDARD DEVIATION TO BE USED TO CALCULATE LOWER CONFIDENCE LIMIT 


985 RANDOM NUMBERS WERE REJECTED PRIOR TO SIMULATION. 
SIMULATION BEGINS WITH RANDOM SEED = .94062208 

1 PHASE TYPES ARE USED. 

1 TIMELINE IS SIMULATED. 


PHASE SEQUENCE TYPE DURATION CUMULATIVE TIME 
1 1 2184.00 2184.00 

PRINTOUT OPTION: MANAGEMENT SUMMARY 
REPAIR OPTIONS; 

PHASE TYPE 1: REPAIR ALLOWED 
100 PERCENT OF REPAIRS ARE PERFORMED ON STATION. 

MISSION ALLOWABLE DOWNTIME IS 273.00 HOURS. 

MTBF MULTIPLIER = 1.00 
MTTR MULTIPLIER = 1.00 
BLOCK III CMA 


EQUIPMENT CHARACTERISTICS 


TYPE 

NOMENCLATURE 

MTBF 

1 

SYCONT 


25000.0 

2 

SC SYN 


12000.0 

3 

DC SYN 


12000.0 

4 

COMREG 


25000.0 

5 

SC PFC 



25000.0 

6 

DC PFC 



25000.0 

7 

DFG 



25000.0 

8 

SFG 



25000.0 

9 

MOD 



25000.0 

10 

D/A 



25000.0 

n 

POWER SUPPLY 

10000.0 

12 

DELAY 



25000.0 

13 

D/A DELAYED 

25000.0 

TYPE 

EQUIPMENT 


1 

1 

18 



2 

2 

19 



3 

3 

20 



4 

4 

12 

21 29 


5 

5 

22 



6 

6 

23 



7 

7 

13 

24 30 


8 

8 

14 

25 31 


9 

9 

15 

26 32 


10 

10 

27 



11 

11 

28 



12 

16 

33 



13 

17 

34 



SPARES COMPLEMENT 


TYPE 

STATION 

COMPLEX 

NETWORK 

1 

0 


1 

0 

2 

0 


1 

0 

3 

0 


1 

0 

4 

0 


2 

0 

5 

0 


1 

0 

6 

0 


1 

0 

7 

3 


2 

0 

8 

3 


2 

2 

9 

0 


2 

0 

10 

0 


1 

0 

11 

0 


1 

0 

12 

0 


1 

0 

13 

0 


3 

0 


BLOCK III CMA 



DUTY 

ADMIN DELAY TIME 

MTTF 

CYCLE 

COMPLEX 

NETWORK 

2.00 

1.00 

24.0 

336.0 

2.00 

1.00 

24.0 

336.0 

2.00 

1.00 

24.0 

336.0 

2.00 

1.00 

24.0 

336.0 

2.00 

1.00 

24.0 

336.0 

2.00 

1.00 

24.0 

336.0 

2.00 

1.00 

24.0 

336.0 

2.00 

1.00 

24.0 

336.0 

2.00 

1.00 

24.0 

336.0 

4.00 

1.00 

24.0 

336.0 

2.00 

1.00 

24.0 

336.0 

2.00 

1.00 

24.0 

336,0 

4.00 

1.00 

24.0 

336.0 


SPARES MULTIPLIER = 1.00 


Fig. 3. TIGER Program output: simulation control specifications and equipment characteristics 




SIMULATION 


MISSION 1 ABORTED IN PHASE SEQ 1 (TYPE 1) AT TIME 

BECAUSE CMA EXCEEDED ALLOWABLE DOWNTIME ( 273.0 HRS.) IN MISSION. 

POWER SUPPLY (EQUIPMENT 11) IS DOWN. 

IT IS NOT REPAIRABLE DUE TO LACK OF SPARES. 

MISSION 5 ABORTED IN PHASE SEQ 1 (TYPE 1) AT TIME 

BECAUSE CMA EXCEEDED ALLOWABLE DOWNTIME ( 273.0 HRS.) IN MISSION. 

DC PFC (EQUIPMENT 6) IS DOWN. 

IT IS NOT REPAIRABLE DUE TO LACK OF SPARES. 

MISSION 17 ABORTED IN PHASE SEQ 1 (TYPE 1) AT TIME 

BECAUSE CMA EXCEEDED ALLOWABLE DOWNTIME ( 273.0 HRS.) IN MISSION. 

SC SYN (EQUIPMENT 2) IS DOWN. 

IT IS NOT REPAIRABLE DUE TO LACK OF SPARES. 

MISSION 24 ABORTED IN PHASE SEQ 1 (TYPE 1) AT TIME 

BECAUSE CMA EXCEEDED ALLOWABLE DOWNTIME ( 273.0 HRS.) IN MISSION. 

SYCONT (EQUIPMENT 1) IS DOWN. 

IT IS NOT REPAIRABLE DUE TO LACK OF SPARES. 

MISSION 31 ABORTED IN PHASE SEQ 1 (TYPE 1) AT TIME 

BECAUSE CMA EXCEEDED ALLOWABLE DOWNTIME ( 273.0 HRS.) IN MISSION. 

DC PFC (EQUIPMENT 6) IS DOWN. 

IT IS NOT REPAIRABLE DUE TO LACK OF SPARES. 

MISSION 40 ABORTED PHASE SEQ 1 (TYPE 1) AT TIME 

BECAUSE CMA EXCEEDED ALLOWABLE DOWNTIME ( 273.0 HRS.) IN MISSION. 

SC SYN (EQUIPMENT 2) IS DOWN. 

IT IS NOT REPAIRABLE DUE TO LACK OF SPARES. 

MISSION 49 ABORTED IN PHASE SEQ 1 (TYPE 1) AT TIME 

BECAUSE CMA EXCEEDED ALLOWABLE DOWNTIME ( 273.0 HRS.) IN MISSION. 

DC SYN (EQUIPMENT 3) IS DOWN. 

IT IS NOT REPAIRABLE DUE TO LACK OF SPARES. 

SC PFC (EQUIPMENT 22) IS DOWN. 

IT WILL COME UP AT TIME 968.1 


1816.6 

1273.9 

1170.6 

1550.1 

759.9 

940.4 

967.4 


Fig. 4. TIGER Program output: management report 




SIMULATION 



START OF MISSION 

] ******************** 



SFG 

(EQUIPMENT 14) FAILED AT TIME 

232.1 


SFG 

(EQUIPMENT 14) CAME UP AT TIME 

232.9 


SYSTEM WENT DOWN AT TIME 232.1 DOWN TIME IS 

.9 HRS. 


POWER SUPPLY 

(EQUIPMENT 28) FAILED AT TIME 548.5 



POWER SUPPLY 

(EQUIPMENT TYPE It) CONSUMED ALL SPARES AT TIME 

548.5 

POWER SUPPLY 

(EQUIPMENT 28) CAME UP AT TIME 

574.7 


MOD 

(EQUIPMENT 32) FAILED AT TIME 

968.9 


MOD 

(EQUIPMENT 32) CAME UP AT TIME 

995.0 


POWER SUPPLY 

(EQUIPMENT 11) FAILED AT TIME 

1543.6 


DC SYN 

(EQUIPMENT 20) FAILED AT TIME 

1906.4 


MISSION ABORTED AT TIME 1816.6 



BECAUSE CMA EXCEEDED ALLOWABLE DOWNTIME (273.0 HRS.) IN MISSION. 


POWER SUPPLY 

(EQUIPMENT 11) IS DOWN. 



IT IS NOT REPAIRABLE DUE TO LACK OF SPARES. 



DC SYN 

(EQUIPMENT TYPE 3) CONSUMED ALL SPARES AT TIME 

1906.4 

DC SYN 

(EQUIPMENT 20) CAME UP AT TIME 

1930.6 


SYSTEM WENT DOWN AT TIME 1543.6 DOWN TIME IS 

387.0 HRS. 


POWER SUPPLY 

(EQUIPMENT 28) WILL FAIL AT TIME 

2512.2 


2184.0 END OF SEQ 

1 (TYPE 1). 



SYSTEM WAS DOWN FOF 

253.4 HRS. IN PHASE SEQ 1 



AND REMAINED DOWN AT END OF PHASE. 



START OF MISSION 

2 ******************** 



DC SYN 

(EQUIPMENT 3) FAILED AT TIME 

428.7 


DC SYN 

(EQUIPMENT TYPE 3) CONSUMED ALL SPARES AT TIME 

428.7 

DC SYN 

(EQUIPMENT 3) CAME UP AT TIME 

454.4 


SYSTEM WENT DOWN AT TIME 428.7 DOWN TIME IS 

25.8 HRS. 


D/A DELAYED 

(EQUIPMENT 17) FAILED AT TIME 

1157.5 


D/A DELAYED 

(EQUIPMENT 17) CAME UP AT TIME 

1185.9 


SYSTEM WENT DOWN AT TIME 1157.5 DOWN TIME IS 

28.4 HRS. 


SC SYN 

(EQUIPMENT 19) WILL FAIL AT TIME 

2629.6 


2184.0 END OF SEQ 

1 (TYPE 1) . 



START OF MISSION 

3 ******************** 



D/A 

(EQUIPMENT 27) FAILED AT TIME 

1270.2 


D/A 

(EQUIPMENT TYPE 10) CONSUMED ALL SPARES AT TIME 

1270.2 

D/A 

(EQUIPMENT 27) CAME UP AT TIME 

1297.4 


D/A 

(EQUIPMENT 27) FAILED AT TIME 

2001.9 


COMREG 

(EQUIPMENT 4) FAILED AT TIME 

2015.4 


COMREG 

(EQUIPMENT 4) CAME UP AT TIME 

2040.9 


SYSTEM WENT DOWN AT TIME 2015.4 DOWN TIME IS 

25.6 HRS. 


MOD 

(EQUIPMENT 15) WILL FAIL AT TIME 

2206.0 


2184.0 END OF SEQ 

1 (TYPE 1). 




Fig. 5. TIGER Program output: engineering report 


FIGURES OF MERIT FOR PHASE SEQUENCE 1 (PHASE TYPE 1): 


BEGINNING; 
IN PHASE: 


END OF PHASE: 


AVAILABILITY (INSTANTANEOUS) = 1.0000 

RELIABILITY = .9240 

AVAILABILITY (AVERAGE) = ,9A65 

READINESS = .9545 

AVAILABILITY (INSTANTANEOUS) = .8960 


FIGURE OF MERIT SUMMARY FOR A GRAND TOTAL OF 250 MISSIONS: 

RELIABILITY = .92400 SPECIFICATION = .9980 

LOWER CONFIDENCE LIMIT (FOR STANDARD DEVIATION = 1.2800) = .9025 
AVAILABILITY (AVERAGE) = .94648 
(INSTANT) = .89600 
READINESS = .95448 

MEAN TIME BETWEEN MISSION FAILURES (MTBMF) = 27910.6 

LOWER CONFIDENCE LIMIT (90 PERCENT) = -5101 .4 (STD DEV = 1 .28) 

VARIANCE = 665157560.0 

SYSTEM MEAN UP TIME, MUT (MTBF ESTIMATOR) = 1288.7 

SYSTEM MEAN DOWN TIME, MDT (MTTR ESTIMATOR) = 21.2 

SIMULATION COMPLETE - SPECIFIED OPTIMUM NUMBER OF MISSIONS WERE RUN. 
SYSTEM FAILS REQUIREMENT. 

BLOCK III CMA 


Fig. 6. TIGER Program output: system reliability, maintainability and 
readiness report 
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EQUIPMENT FAILURES AND CORRECTIVE MAINTENANCE (CM) SUMMARY 



EQUIPMENT 

TOTAL 

MISSION AVERAGES 

NUMBER 

TYPE 

NOMENCLATURE 

FAILURES 

FAILURES 

CM MANHOURS 

1 

1 

SYCONT 

20 

.080 

.160 

2 

2 

SC SYN 

41 

.164 

.328 

3 

3 

DC SYN 

49 

.196 

,392 

4 

4 

COMREG 

25 

.100 

,200 

5 

5 

SC PFC 

22 

.088 

,176 

6 

6 

DC PFC 

21 

.084 

.168 

7 

7 

DFG 

19 

.076 

.152 

8 

8 

SFG 

20 

.080 

.160 

9 

9 

MOD 

15 

.060 

.120 

10 

10 

D/A 

20 

.080 

,320 

11 

11 

POWER SUPPLY 

51 

.204 

,408 

12 

4 

COMREG 

19 

.076 

.152 

13 

7 

DFG 

27 

.108 

.216 

14 

8 

SFG 

17 

.068 

.136 

15 

9 

MOD 

25 

.100 

.200 

16 

12 

DELAY 

17 

.068 

.136 

17 

13 

D/A DELAYED 

12 

.048 

.192 

18 

1 

SYCONT 

23 

.092 

.184 

19 

2 

SC SYN 

50 

.200 

.400 

20 

3 

DC SYN 

40 

.160 

.320 

21 

4 

COMREG 

20 

.080 

.160 

22 

5 

SC PFC 

24 

.096 

.192 

23 

6 

DC PFC 

18 

.072 

.144 

24 

7 

DFG 

13 

.052 

.104 

25 

8 

SFG 

24 

.096 

.192 

26 

9 

MOD 

14 

.056 

.112 

27 

10 

D/A 

20 

.080 

.320 

28 

11 

POWER SUPPLY 

45 

.180 

.360 

29 

4 

COMREG 

22 

.088 

.176 

30 

7 

DFG 

23 

.092 

.184 

31 

8 

SFG 

21 

.084 

.168 

32 

9 

MOD 

24 

.096 

.192 

33 

12 

DELAY 

16 

.064 

.128 

34 

13 

D/A DELAYED 

14 

.056 

.224 




831 

3.324 

7.176 


BLOCK III CM A 


Fig. 7. TIGER Program output: equipment failures and corrective 
maintenance summary 


NUMBER OF SPARES PER MISSION 


EQUIPMENT 


STATION 
STOCK USAGE 


COMPLEX 
STOCK USAGE 


NETWORK 
STOCK USAGE 


TYPE 

NOMENCLATURE 


MAX 

AVG 


MAX 

AVG 


MAX 

AVG 

1 

SYCONT 

0 

0 

.00 

1 

1 

.16 

0 

0 

.00 

2 

SC SYN 

0 

0 

.00 

1 

1 

.32 

0 

0 

.00 

3 

DC SYN 

0 

0 

.00 

1 

1 

.32 

0 

0 

.00 

4 

COMREG 

0 

0 

.00 

2 

2 

.34 

0 

0 

.00 

5 

SC PFC 

0 

0 

.00 

1 

1 

.17 

0 

0 

.00 

6 

DC PFC 

0 

0 

.00 

1 

1 

.15 

0 

0 

.00 

7 

DFG 

3 

3 

.33 

2 

0 

.00 

0 

0 

.00 

8 

SFG 

3 

3 

.33 

2 

0 

.00 

2 

0 

.00 

9 

MOD 

0 

0 

.00 

2 

2 

.30 

0 

0 

.00 

10 

D/A 

0 

0 

.00 

1 

1 

.14 

0 

0 

.00 

11 

POWER SUPPLY 

0 

0 

.00 

1 

1 

.33 

0 

0 

.00 

12 

DELAY 

0 

0 

.00 

1 

1 

.13 

0 

0 

.00 

13 

D/A DELAYED 

0 

0 

.00 

3 

1 

.10 

0 

0 

.00 


BLOCK III CMA 


Fig. 8. TIGER Program output: number of spares per mission 




EQUIPMENT CONTRIBUTIONS TO SYSTEM UNAVAILABILITY 


EQUIPMENT 

DOWNTIME 

UNAVAILABILITY 

PERCENT 

NUMBER 

TYPE 

NOMENCLATURE 

HOURS 

FRACTION 

CONTRIBUTION 

2 

2 

SC SYN 

9181.3 

.01682 

31.4 

11 

11 

POWER SUPPLY 

5521.9 

.01011 

18.9 

6 

6 

DC PFC 

3346.0 

.00613 

11.5 

3 

3 

DC SYN 

2995.1 

.00549 

10.2 

10 

10 

D/A 

2450.6 

.00449 

8.4 

1 

1 

SYCONT 

1714.5 

.00314 

5.9 

5 

5 

SC PFC 

859.3 

.00157 

2.9 

15 

9 

MOD 

720.3 

.00132 

2.5 

4 

4 

COM REG 

636.1 

.00116 

2.2 

12 

4 

COMREG 

501.2 

.00092 

1.7 

16 

12 

DELAY 

397.4 

.00073 

1.4 

9 

9 

MOD 

369.6 

.00068 

1.3 

17 

13 

D/A DELAYED 

334.1 

.00061 

1.1 

13 

7 

DFG 

63.5 

.00012 

.2 

7 

7 

DFG 

47.4 

.00009 

.2 

8 

8 

SFG 

44.8 

.00008 

.2 

14 

8 

SFG 

38.9 

.00007 

.1 


BLOCK Hi CMA 


Fig. 9. TIGER Program output: equipment contributions to system unavailability 



EQUIPMENT CONTRIBUTIONS TO SYSTEM UNRELIABILITY 


EQUIPMENT 

NUMBER OF 

UNRELIABILITY 

PERCENT 

NUMBER 

TYPE 

NOMENCLATURE 

FAILURES 

FRACTION CONTRIBUTION 

2 

2 

SC SYN 

6 

.02400 

31.6 

11 

11 

POWER SUPPLY 

4 

.01600 

21.1 

1 

1 

SYCONT 

2 

.00800 

10.5 

6 

6 

DC PFC 

2 

.00800 

10.5 

10 

10 

D/A 

2 

.00800 

10.5 

3 

3 

DC SYN 

2 

.00800 

10.5 

5 

5 

SC PFC 

1 

.00400 

5.3 


Fig. 10. TIGER Program output: equipment contributions to system unreliability 
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Special Activity Utilization of GDSCC Antennas During 1981 

E. B. Jackson 
Goldstone Operations Section 


In addition to direct spacecraft support, the GDSCC antennas also support '*speciaV* 
activities. These activities can be categorized as Advanced Systems Program, Radio 
Astronomy Program, Crustal Dynamics Program, and Operations and Support Activities. 
This article briefly discusses the goals and categorizes each of the activities that received 
tracking support at Goldstone during 1981. All Goldstone stations (DSSs 11, 12, 13 and 
14) provided a total of 2353.55 hours of support to **special” activities during the period 
January through December 1981. 


I. Introduction 

The antennas at the Goldstone Deep Space Communica- 
tions Complex (GDSCC), in addition to providing telemetry, 
ranging, and commanding support to spacecraft projects, also 
support a number of other tasks which require “tracking.” 
These special activities are categorized in this report into 
Advanced Systems Program, Radio Astronomy Program, 
Crustal Dynamics Program, and Operations and Support 
Activities. This report discusses the special activities which 
were supported with tracking by the GDSCC antennas during 
the period January through December, 1981. 

II. Advanced Systems Program 

Activities in this category are funded by NASA through the 
Office of Space Tracking and Data Systems (OSTDS). This is 
the development work which culminates in provision, to the 
Deep Space Network (DSN) stations, of new capability with 
which each spacecraft project’s needs can be met. Related 
activities are associated and organized into Research and 


Technology Operations Plans (RTOPs). GDSCC activities in 
this category are tabulated in Table 1 with hours and stations 
involved. 

A. RTOP 60: Radio Metric Technology Development 

The goal of this effort is the development and demonstra- 
tion of advanced radio metric systems for navigation and radio 
science, with a specific angular accuracy goal of 50 nanoradians 
for delta VLBI measurements. 

In general, using VLBI, this RTOP generates a catalog of 
suitable (strength and location) extragalactic radio sources 
(EGRSs). Using this catalog of EGRSs, the ephemerides of the 
planets, particularly Mars, are tied into the “radio reference 
frame.” Then techniques are developed for differential mea- 
surements between spacecraft and EGRSs which locate the 
spacecraft within the radio reference frame and hence the 
spacecraft with reference to the target planet(s). Table 1 lists 
the various activities within the RTOP for which tracking 
support was provided at GDSCC during 1981. 
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B. RTOP 65: Antenna Systems Development 

In recognition of the necessity of improving antenna 
performance to provide increased planetary communications 
capability, this activity’s goal is to enhance antenna tech- 
nology. Specifically, this RTOP looks at possible improve- 
ments in electronic-microwave capabilities and mechanical- 
structural design. Application of technology to provide higher 
frequency operation of existing antennas, along with more 
sophisticated techniques of antenna pointing, is also a part of 
this activity. 

As part of the analysis performed in this area, techniques 
for measurement of antenna gain were refined and used to 
make extensive measurements at DSS 12. These antenna gain 
measurements were used to quantify and evaluate the changes 
in microwave performance resulting from main reflector 
surface alignment and changes in subreflector position as a 
function of antenna attitude. 

The particular activity supported during 1981 at DSS 13 
was “Analytical Techniques and Procedures,” which has as 
its objective the development of analysis software for applica- 
tion to antenna structural and mechanical performance. (An 
innovative application used a cooled infrared camera to obtain 
thermal measurements of the DSS 13 26-m antenna backup 
structure. These measurements supplied information on 
possible thermally induced structural deformations which was 
not heretofore available.) This activity also aims to extend the 
capability for automatic antenna structural design optimiza- 
tion and organize existing software to simplify maintenance 
and execution of improvements. Again, Table 1 depicts track- 
ing support during 1981. 

C. RTOP 68: Station Monitor and Control Technology 

Recognizing that a substantial portion of the operating 
cost for DSN stations is allocated to personnel, this RTOP 
seeks to develop automation technology that will allow entire 
station operation from a remote point and prove the reliabil- 
ity, safety, and efficiency of unattended operation. Addi- 
tionally, this RTOP plans to develop the data base with which 
reliability, costs, and productivity can be monitored, using 
DSS 13 as a demonstration unattended station. 

Table 1 lists only tracking support to this activity; substan- 
tial other work was also ongoing at DSS 13 during the year as 
the equipment used to make DSS 13 operate in the unattended 
mode evolved and additional systems were added to those 
controlled remotely. 

D. RTOP 70: High-Speed Signal Processing Research 

The objective of this RTOP is to develop high-speed digital 
signal processing techniques for use in the DSN. A major part 


of this task is the development of a test bed as a research tool 
used to explore high-speed techniques. The test bed is to be 
used for demonstrations of various development capabilities, 
such as high-performance array processors and wideband 
correlation subsystems. This work also demonstrates commer- 
cial LSI signal processing building blocks and utilizes promis- 
ing commercial design systems. 

The particular activity supported at GDSCC during 1981 
was Radar Data Acquisition, which uses various planets and 
asteroids as targets for demonstration of signal processing 
techniques. 

In Table 1 are listed the various targets which were explored 
in 1981 as the demonstration test bed was used to develop sig- 
nal processing techniques. 


III. Radio Astronomy Program 

In general, activities in this category are either sponsored 
by the Office of Space Science (OSS) or authorized by the 
Radio Astronomy Experiment Selection (RAES) panel. 

Experimenters within NASA who desire OSS support of 
their activities submit an observing plan in the form of an 
RTOP. If approved, OSS requests OSTDS to provide support, 
and time is made available on suitable facilities of the DSN. 
Experimenters from universities who have observing plans 
that require NASA facilities submit their observing plan to 
the RAES panel. If approved, time is made available through 
an agreement whereby NASA makes available a percentage 
of the operational time on DSN facilities. Table 2 depicts 
the support at GDSCC in this area during 1981. 

A. OSS-Sponsored Activities 

1. Pulsar rotation constancy. This experiment seeks to 
monitor short-term variations in the period of the relatively 
young VELA pulsar (PSR 0833-45). Additionally, by observ- 
ing 23 older, more stable pulsars, this experiment hopes to 
obtain data that will provide further tests of the hypothesis 
that pulsars are neutron stars resulting from supernova explo- 
sions, which impart high velocities to the resulting pulsar. 

2. Planetary radio astronomy. This experiment has the dual 
objectives of studying the properties of the planet Jupiter’s 
radio emission and measurement of the thermal emission from 
the atmosphere of the outer planets. A new objective is to 
search for a link between variations in the lo plasma Torus and 
Jupiter’s radio emissions. These observations will aid construc- 
tion of improved models of Jupiter’s radiation belt environ- 
ment as well as atmospheric models of the outer planets, par- 
ticularly Uranus. 


122 


3. SETI Receiving System Stability Measurements. In 
anticipation of a SETI observing program (which never mater- 
ialized), the phase and amplitude stability of the various 
receiving systems at DSS 14 was measured. 

B. JPL Director’s Discretionary Fund Sponsored 
Activities 

1. Science utilization of the Radio Frequency Interference 
Surveillance Subsystem (SURFISS). As part of the task of 
detection and characterization of radio frequency interference, 
the DSN has developed a unique capability to perform high- 
resolution spectrum analysis. This van-mounted subsystem can 
perform real-time spectral analysis of a 20-MHz band of fre- 
quencies, using 2^^ channels, resulting in a resolution of 
305 Hz, The objective of this task was to develop computer 
programming and observing techniques which would best 
make use of this high-resolution, wide bandwidth capability 
to perform radio astronomical observations. In particular, 
the RFI subsystem, with altered software, would be used to 
observe fine structure and narrow spectral features in radio 
astronomical objects. 

The RFI subsystem, coupled to the 26-m antenna, was 
used at DSS 13, to observe, among others, the H142a line in 
Ml 7. This line was observed at approximately 2272.3 MHz. 

C. RAES Panel Sponsored Activities 

1. VLBI investigation of SS-433 (RA-175). SS-433 has a 
bizarre optical spectrum that exhibits three sets of emission 
lines. One of these three emission lines shows near-zero radial 
velocity, while the other two show large and variable shifts 
to the blue and the red. The variability is periodic, with a 
period of about 164 days. 

This experiment conducts regular VLBI observations of 
this source, using many simultaneous baselines. These multiple 
baselines should enable determination of the source angular 
size and structure, which wiU aid niodeling. These observa- 
tions show that the compact radio structure is also variable, 
with a period of 164 days. The involved stations at GDSCC 
include DSS 13 and DSS 14, along with a number of other 
radio observatories in the United States and Europe. 

2. VLBI investigation of “twin” quasistellar objects (QSOs) 
0957 + 561A, B (RA 176). These QSOs, designated 0957 + 
561A, B, are separated by only 6 arc seconds in angle and have 
equal optical redshifts and remarkable similarity in their 
optical spectra. This similarity is inconsistent with chance 
alone, so the objects are in some manner physically asso- 
ciated. One theory proposes that there is only one object, 
whose radiation is gravitationally bent about an intervening 
massive lens, a “gravitational lens.” 


This experiment uses VLBI observations with many simul- 
taneous baselines to provide data with which detailed structure 
maps can be prepared. Repeated observations could monitor 
changes in the radio brightness distribution. At GDSCC, the 
involved antenna is DSS 14. 

3. Milliarcsecond nuclei in quasars and galaxies (RA-177). 
The 64-m antennas at Goldstone, California, and Madrid, 
Spain, connected as a Very Long Baseline Interferometer 
(VLBI), have a detection sensitivity of 1 mJy at S-band and 
2 mJy at X-band. This task uses this sensitive interferometer 
to determine the frequency of occurrence of milliarcsecond 
radio nuclei in radio quasars and radio galaxies. These two 
antennas, used in the VLBI configuration, generate detailed 
maps of the nuclei of these extragalactic objects. A knowledge 
of the statistical occurrence and properties of these nuclei 
will aid in a more detailed understanding of the role these 
nuclei play in the energetics and evolution of extragalactic 
objects. 


IV. Crustal Dynamics Program 

A Crustal Dynamics Program goal is the demonstration of 
the capability of VLBI systems to make highly accurate 
geodetic measurements. The particular activity supported at 
GDSCC during 1981 is Mobile VLBI Field Operations, The 
thrust of the Mobile VLBI activity is to demonstrate the 
geodetic performance of highly mobile (antenna size approxi- 
mately 4 meters) VLBI stations while providing accurate data 
of significant geophysical interest. By moving around California 
and making VLBI measurements with base stations located at 
GDSCC, Owens Valley Radio Observatory, and other fixed 
points, the Mobile VLBI activity is able to accumulate data on 
regional deformation and strain accumulation, particularly as 
associated with the San Andreas Fault. 

These measurements are usually scheduled for periods of 
7 to 10 days, in excess of 24 hours observing at a time, and a 
number of suitable radio sources are observed simultaneously 
and repeatedly by all the stations involved. During the period 
January 1 through DecemberSl, 1981, DSS 13 (Venus 
Station) provided 418.75 hours of tracking support to these 
measurements of cmstal deformation along the San Andreas 
Fault. 


V. Operations and Support Program 

The maintenance of a capability to observe at various 
frequencies, on a routine basis, involves a number of activities 
which do not easily fit any of the specific programs supported 
at GDSCC, but these activities are essential to continued sup- 
port of any or all programs. For the purposes of this progress 
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report, these activities are grouped under the heading “Opera- 
tions and Support.” These activities are tabulated in Table 3, 
along with hours and stations involved. 

A. Antenna Boresight Offset Measurements 

At DSS 13, where developmental Cassegrain feedcones 
are routinely used for program support, continued knowledge 
of the offset between the antenna mechanical axis and the 
antenna electromagnetic axis, as a function of antenna attitude, 
is essential. By use of suitable radio calibration sources, on 
which the antenna is “boresighted” (usually at 15-minute 
intervals of time), a table of axis offsets as a function of 
antenna attitude can be generated. This table can then be used 
to provide assurance of accurate “blind pointing” on objects 
such as weak radio sources where no real-time indication of 
correct pointing can be obtained. 

B. High Resolution Analysis, Transmitter(s) 

Radiated Spectrum 

As part of the overall electromagnetic compatibility pro- 
gram at DSCC, the actual radiated spectrum of the various 
transmitters at DSS 14 was measured. These measurements 
utilized the Department of Commerce radio frequency mea- 
suring van, coupled with the JPL developed high resolution 
spectrum analysis subsystem, usually referred to as the “RFI 
Van.” Measurements were made on the 20- and 100-kW 
S-band transmitters and the 400-kW X-band transmitter. 
Partial measurements were also made on the 400-kW S-band 
transmitter. These measurements characterized the radiated 
spectrum of each of the transmitters, including carriers and 
associated spurious signals, with a maximum resolution of 
305 Hz. 

C. Precision Signal Strength Measurements, ISEE-3 

Signal strength measurements of ISEE-3, made in the 
conventional fashion (using automatic gain control voltage 
as an analog of signal strength) were not of appropriate preci- 
sion or resolution. Using equipment normally used for radio 
astronomy activities, DSS 12 and 13 made direct measurement 
of the power ratio in a harrow bandwidth (10 kHz nominal) 
with and without the spacecraft signal present. With knowl- 
edge of the receiving system temperature, the spacecraft 
signal strength can be accurately computed from these power 
ratios. 


D. Electromagnetic Compatibility Testing (Sodium 
Lamp) 

Nighttime illumination around the antennas is essential, 
but should be of minimum energy consumption. Sodium vapor 
lighting is efficient but the potential for radio frequency inter- 
ference to the DSN antennas was unknown. At DSS 13, the 
sodium lamp, in the fired and unfired state, was positioned 
over the antenna feedhorn and, using specialized equipment, 
the possible level of interference was' determined. It was con- 
cluded that lamps of this type are suitable for ground installa- 
tion around DSN antennas. 

E. Precision Signal Strength Measurements, 
ALSEP-16 

Antenna gain measurements can be, and are, made using 
radio calibration sources. However, from time to time, it is 
desirable to have available a far-field source which radiates a 
strong, narrrow-band, stable signal. The radio transmitters left 
on the Moon by the Apollo astronauts are quite suitable for 
this purpose, but occasional characterization of the signals is 
necessary. These measurements were to determine the con- 
tinued suitability of ALSEP-16 as a stable far-field signal 
source. On the DSS 13 26-m antenna, ALSEP-16 was mea- 
sured at -1 10 dbm signal level. 

F. VLBI Software Development 

When performing VLBI measurements, which depend upon 
correlation of two (or more) signals, real-time information 
about proper functioning of the station is desirable. Utilizing 
RF signals generated at the DSS 14 antenna, modifications 
were made to in-use software to enhance its utility for real- 
time performance evaluation. 

G. Antenna Location Determination, National 
Geodetic Survey 

VLBI measurements of vector distances between partici- 
pating stations require accurate knowledge of the geodetic 
location of a reference point on the antennas used to make 
these observations. On the DSS 13 antenna, this point is the 
intersection of the azimuth and elevation axes. Personnel 
from the National Geodetic Survey agency required access 
to the DSS 13 26-m antenna in order to tie this axis intersec- 
tion into the continental grid locations. To aid these measure- 
ments, the 26-m antenna was intermittently positioned as 
required over a period of a few weeks. 
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Table 1. Advanced Systems Program utilization of GDSCC antennas during 1981 


Station and tracking hours provided 

Activity 11 12 13 14 


RTOP 60 radio metric technology development 


Mark IV radio metric system development and 


101.75 

5.25 

demonstration 

VLBI radio source identification 


102.00 

16.50 

5 nanoradian VLBI system analysis and modeling 


52.00 


Microwave phase delay calibration 


44.75 

14.25 

Clock synchronization VLBI 


14.25 


RTOP 65 antenna systems development 




Analytical techniques and procedures 

140.75“ 

33.25 


RTOP 68 station monitor and control technology 




DSS 13 S/X unattended system development 


22.50 


RTOP 70 high-speed signal processing research 




Radar data acquisition 




Asteroid Melpomene 



6.00 

Planet Mercury 


4.25 

58.00 

Planet Saturn and Rings of Saturn 



20.25 

Callisto (satellite of Jupiter) 



29.25 

Totals 

140.75 

374.75 

149.50 


^Antenna gain measurements. 



Table 2. Radio astronomy utilization of GDSCC antennas during 1981 



Station and tracking hours provided 

Activity 

11 

12 

13 

14 

OSS sponsored activities 





Pulsar rotation constancy 

117.00 

3.25 

184.75 

81.00 

Planetary radio astronomy 



207.25 

30.50 

SETI receiving system stability measurements 




26.50 

JPL Director’s Discretionary Fund sponsored activities 





Science utilization of the Radio Frequency 



45.25 


Interference Surveillance Subsystem (SURFISS) 





RAES panel sponsored activities 





VLBI investigation of SS-433-RA 175 



57.25 

26.00 

VLBI investigation of “twin” quasistellar objects 



- 

52.25 

(QSOs) 0957 + 561 A, B and 1038 + 528A, B- 
RA 176 





Milliarcsecond nuclei in quasars and galaxies-RA 177 




15.50 

Totals 

117.00 

3.25 

494.50 

231.50 


Table 3. Station operations and associated activity utilization of GDSCC antennas during 1981 


Activity 

Antenna boresight offset measurements 

High resolution analysis, transmitter(s) 
radiated spectrum 

Precision signal strength measurements, ISEE-3 
Electromagnetic compatability testing (sodium lamp) 
Precision signal strength measurements, ALSEP-16 
VLBI software development 

Antenna location determination, National Geodetic 
Survey 


Station and tracking hours provided 
U 12 13 m” 

80.50 

17.00 

5.00 4.25 

5.00 
2.75 

2.00 

171.25 


5.00 263.75 19.00 


Totals 
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Maintenance and Operations Tasks Accomplished 
at DSS 12 During the Antenna Panel 
Replacement Downtime 

L Butcher, T. Jonas, and W. Wood 
Goldstone Operations Section 


The heavy schedule of tracking activities at the Echo Deep Space Station (DSS 12) 
prevents some time-consuming maintenance tasks from being performed. Careful coor- 
dination prior to and during a mandatory task (antenna panel replacement) made it 
possible to do a large number of unrelated tasks that ordinarily would have to be de- 
ferred, It is not the intent to relate the activities of the antenna panel replacement but to 
describe the maintenance and operations tasks accomplished during the downtime. 


I. Introduction 

An eight-week period, from 22 February to 19 April 
1981, was set aside to correct antenna gain deficiencies in the 
34-meter antenna at Goldstone, California. To take advantage 
of the station availability during the antenna downtime, the 
Goldstone Maintenance and Integration Unit worked closely 
with DSN Engineering (Section 355). This coordination made 
time available to do things that would not normally be done 
because of the heavy tracking commitments. 

II. Maintenance Activities 

The following maintenance and operations tasks were 
accomplished during the antenna panel replacement down- 
time: 

(1) Seven Engineering Change Orders were implemented: 


Number 

Title 

Subsystem 

79.006 

Symbol Synchronizer Assembly 
Upgrade 

40.3 

79.152 

Host Processor Installation 

Various 

79.159 

Radio Frequency Power Supply 
Monitor 

2.1 

79.269 

High Voltage Cable Isolation 

2.0 

79.272 

Angle Encoder Translator 
Damage 

6.7 

79.292 

Add Radiometer Experiment 

Various 

80.205 

Generator Control Logic 
Modification 

39.2 
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(2) Coriective maintenance had to be carefully scheduled 
and accomplished on a not-to-interfere with the 
antenna panel replacement work basis. 

(a) All Frequency and Timing System cables were 
checked and rerouted to eliminate a long standing 
problem of intermittent failures. Included in this 
work was the retermination of system ground 
cables to reduce ground loop problems that had 
existed for years. 

(b) The coolant lines at the 20-kW S-band transmitter 
were disassembled and resealed to prevent leaks. 

(c) The hour angle and declination axes cable wraps 
were reworked to eliminate existing and future 
cable damage. 

(d) Defective wiring in the backup S*band maser 
control cable was repaired by using spare wires 
within the existing cable. 

(e) A large amount of corrective maintenance was 
accomplished on the Digital Instrumentation 
Assembly. 

(3) Workmanship Assurance Discrepancy Report work. 
During the antenna downtime, 26 discrepancy reports, 
containing a total of 59 separate items, were corrected 
and cleared. Deep Space Station 12 is now at an alltime 
low in the number of outstanding Workmanship 
Assurance Discrepancy Reports. 


(4) Systems Testing: 

(a) To prepare for expected DSN support of the 
International Sun-Earth Explorer III (ISEE-3) 
satellite ranging compatibility tracks, the power 
klystron from the 20-kW transmitter was removed 
and tested at the Goldstone Microwave Test 
Facility. Data for retuning the klystron were 
obtained to permit rapid frequency changing 
during later operational testing. 

(b) Following the completion of the antenna panel 
work, seven radio source star tracks were carried 
out to determine antenna gain and to realign the 
antenna subreflector for optimum performance. 

(c) To assure full and complete performance of all 
station systems, a complete set of System Perfor- 
mance Tests (SPTs) was performed for telemetry, 
command and radio metric data. 

III. Summary 

With careful planning and coordination with the Engineer- 
ing Task Manager, time-consuming maintenance tasks were 
accomplished during the major downtime activity. These 
tasks could not have been accomplished during the normally 
scheduled short maintenance periods because of tracking 
commitments. All scheduled maintenance and operations 
activities were successfully accomplished on schedule. 
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